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^, ' Abstract. Half a century has passed since the foundation of nuclear astrophysics. 

Since then, this discipline has reached its maturity. Today, nuclear astrophysics 

constitutes a multidisciplinary crucible of knowledge that combines the achievements in 

I . theoretical astrophysics, observational astronomy, cosmochemistry and nuclear physics. 

fT^ I New tools and developments have revolutionized our understanding of the origin 

OQ ' of the elements: supercomputers have provided astrophysicists with the required 

^ J , computational capabilities to study the evolution of stars in a multidimensional 

r~^ ' framework; the emergence of high-energy astrophysics with space-borne observatories 

^"^ ■ has opened new windows to observe the Universe, from a novel panchromatic 

perspective; cosmochemists have isolated tiny pieces of Stardust embedded in primitive 

meteorites, giving clues on the processes operating in stars as well as on the way matter 

condenses to form solids; and nuclear physicists have measured reactions near stellar 

energies, through the combined efforts using stable and radioactive ion beam facilities. 

j^ ' This review provides comprehensive insight into the nuclear history of the Universe 

and related topics: starting from the Big Bang, when the ashes from the primordial 

explosion were transformed to hydrogen, helium, and few trace elements, to the rich 

variety of nucleosynthesis mechanisms and sites in the Universe. Particular attention 

is paid to the hydrostatic processes governing the evolution of low-mass stars, red 

giants and asymptotic giant-branch stars, as well as to the explosive nucleosynthesis 

occurring in core-collapse and thermonuclear supernovae, 7-ray bursts, classical novae. 

X-ray bursts, superbursts, and stellar mergers. 



(Some figures in tliis article are in color only in the electronic version) 
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1. Introduction 

Most of the ordinary (visible) matter in the Universe, from a tiny terrestrial pebble 
to a giant star, is composed of protons and neutrons (nucleons). These nucleons can 
assemble in a suite of different nuclear configurations. The masses of such bound nuclei 
are actually smaller than the sum of the masses of the free individual nucleons. The 
difference is a measure of the nuclear binding energy that holds the nucleus together, 
according to Einstein's famous equation E = Amc^, where c is the speed of light and 
Am denotes the mass difference between the bound nuclei and the free nucleons. For 
example, the nuclide ^^Fe is among the most tightly bound species, with a binding 
energy near 8.8 MeV per nucleon. There are 82 elements that have stable isotopes, all 
the way from H to Pb (except for Tc and Pm). Several dozen elements have only unstable 
isotopes, naturally abundant or artificially produced in nuclear physics laboratories (the 
super-heaviest nucleus discovered to date has 118 protons and a half-hfe of 0.89 ms [1]). 
The distribution of the nuclide abundances in the Solar System, plotted as a function of 
the number of nucleons (mass number). A, is shown in Fig. [H It displays a complicated 
pattern, with H and He being by far the most abundant species. The abundances 
generally decrease for increasing mass number. Exceptions to this behavior include the 
light elements Li, Be and B, which are extremely underabundant with respect to their 
neighbors, a pronounced peak near ^^Fe (the so-called iron-peak) , and several noticeable 
maxima in the region oi A> 100, for the distributions of both odd- and even-A species. 

The reasons for favoring the presence of some species and for the specific pattern 
inferred for the Solar System abundances have spawned large debates. Following earlier 
attempts, Hoyle proposed in 1946 that all heavy elements are synthesized in stars [3j. A 
major breakthrough was achieved with the detection of technetium, an element with no 
stable isotopes, in the spectra of several S stars by Merrill [1] in the early 1950s. Since the 
technetium isotopes have rather short half-lives on a Galactic time scale, the obvious 
conclusion was that technetium must have been produced in the observed stars and 
thus nucleosynthesis must still be going on in the Universe. Shortly after, two seminal 
papers that provided the theoretical framework for the origin of the chemical species 
elements, were published (almost exactly a century after Darwin's treatise on the origin 
of biological species) by Burbidge, Burbidge, Fowler & Hoyle [5], and independently, by 
Cameron [6]. These works identified stars as Galactic factories that synthesize atomic 
nuclei via nuclear reactions. The stars return this processed material to the interstellar 
medium through various hydrostatic or explosive means during their evolution and 
thereby enrich the interstellar medium in metals, that is, any species heavier than He 
according to the astronomical jargon. Out of this matter new generations of stars are 
born and a new cycle of nuclear transmutations will be initiated. In fact, without 
concourse of the many nuclear processes that take place in the stars and that are 
discussed in this review, the Universe would have remained a chemically poor place, 
and probably no life form would have ever emerged. 

This review summarizes more than 50 years of progress in our understanding of the 
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Figure 1. Solar System isotopic abundances versus mass number, A. Abundances are 
based on meteoritic samples as well as on values inferred from the solar photosphere 
(Sect. 2). Values taken from Lodders |2]- 



nuclear processes that shaped the current chemical abundance pattern of the Universe. 
The structure of the paper is as follows. In Section 2, we briefly address the main 
aspects associated with the determination of the Solar System abundances. This is 
followed by an overview on Big Bang nucleosynthesis (Sect. 3), and by an account of 
the nuclear processes governing low-mass stars (Sect. 4). Particular attention is paid 
to the Standard Solar model and to lessons learned from neutrino physics. In Sect. 
5, we outline the basic properties of red giant and AGB stars, with emphasis on the 
s-process. The series of explosive processes that characterize core-collapse supernovae 
and 7-ray bursts (e.g., p-, v-, up-, a-, and r-processes) are reviewed in Sect. 6. This is 
supplemented with an overview of different astrophysical scenarios that involve binary 
stellar systems, such as classical and recurrent novae (Sect. 7), type la supernovae 
(Sect. 8), X-ray bursts and superbursts (Sect. 9), and stellar mergers (Sect. 10). We 
briefly address as well non-stellar processes, such as spallation reactions (Sect. 11), and 
conclude with a number of key puzzles and future perspectives in this field (Sect. 12). 
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2. An insight into Solar System abundances 

The determination of the chemical composition of a distant star is not a trivial task. 
It basically relies on the analysis of the electromagnetic radiation emitted at different 
wavelengths (so-called spectra). For the case of the Solar System abundances, several 
sources of information, that combine radiation and matter samples, are used to this end: 
the Sun (mainly photospheric absorption lines), meteorites, planets and moons (mainly 
Earth and Moon). Recently, several space probes have for the first time carried material 
from the interplanetary medium to Earth. NASA's Genesis mission collected about 0.4 
mg of solar wind particles for two years, using 250 hexagonal (10 cm wide) wafers made 
of silicon, sapphire, diamond, and gold. During descent in 2007 the parachutes never 
opened and the capsule crashed, but fortunately the samples could be retrieved and 
analyzed. NASA's Stardust mission, in turn, flew through the tail of comet Wild-2 in 
2004 and captured dust particles at a speed of 6.1 km s^^ both in aerogel and in Al 
foils. Those samples were brought to Earth in 2006. 

Inferring the elemental composition of the Solar System from this variety of sources 
becomes a complex enterprise. First, Earth does not show a homogeneous chemical 
composition. Second, meteorites do not constitute a chemically unique class of objects, 
but different types exhibit large abundance differences. The usual approach to infer the 
bulk of Solar System abundances relies on a very rare and primitive type of meteorites: 
CI carbonaceous chondrites, representing the least modified and chemically unaltered 
class of meteorites. Named after the Ivuna meteorite, the CI carbonaceous chondrite 
class includes only four other members: Orgueil, Alais, Tonk and Revelstoke. Because 
of the extreme similarities of the abundances in CI chondrites (determined through 
mass spectroscopy, wet chemistry. X-ray fiuorescence, and neutron activation) and those 
inferred from the solar photosphere (for non- volatile elements), it is believed that this 
class of meteorites represents pristine material from which the Solar System formed 
about 4.6 Gyr ago. Unfortunately, noble gases and other very volatile elements do not 
easily condense into solids. Hence, values inferred from solar photospheric spectra are 
used instead to account for the C, N or O abundances, while other specific techniques 
are adopted for some special elements, such as Ar, Kr, Xe, or Hg. Moreover, the 
meteoritic abundances have to be expressed relative to an element other than H. Silicon 
is the species traditionally chosen, with a defined number abundance of N(Si) = 10^. 
This requires renormalization of the meteoritic values, such that the photospheric 
and meteoritic Si abundances agree. The procedure guarantees that both abundance 
determinations rely on the same absolute scale. 

It is worth noting that photospheric abundances face two major drawbacks (besides 
the fact that Li is depleted in the solar atmosphere). On the one hand, the precision 
achieved in photospheric abundance determinations is much lower than what can be 
obtained using mass spectroscopy measurements of meteoritic samples in laboratories. 
On the other hand, photospheric abundances are not directly measured but inferred 
from observations. This requires state-of-the-art models of the solar atmosphere as 
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well as an accurate account of spectrum-formation processes. As an example, the 
recent improvements in atomic and molecular transition probabilities, coupled with the 
use of three-dimensional (3-D) hydrodynamic solar atmosphere models that relax the 
assumption of a local thermodynamic equilibrium in spectral line formation, translated 
into a severe reduction of the abundances of C, N and O compared to previous estimates 
[3 E]. Furthermore, the existing 3-D models proved successful in reproducing a suite of 
characteristic solar features that 1-D models cannot account for, including the observed 
solar granulation topology, typical length- and time-scales, convective velocities, and 
specific profiles of metallic lines. One of the major sources of uncertainty still remaining 
is the lack of accurate cross sections for excitation and ionization. 

Unfortunately, photospheric abundance determinations of a few elements are 
not possible through spectroscopy of the quiet Sun, and other techniques must be 
implemented. For instance, this is the case for the F and CI abundances, inferred 
from infrared spectroscopy of sunspots (through molecular lines of HF and HCl, 
respectively). With regard to noble gases, no photospheric lines are available because 
of their high excitation potentials. Indirect determinations, using helioseismology (He), 
X-ray and ultraviolet spectroscopy of the solar corona and solar flares (Ne), or solar 
wind measurements (Ne, Ar), among others, are used to this end. 

It is remarkable that while variations in the elemental composition have been found 
within the different existing sources, that can be attributed to a variety of physical, 
chemical and geological processes, the bulk isotopic composition of the Solar System 
is extremely homogeneous. This circumstance turns out to be extremely useful since 
isotopic ratios of the Sun cannot be derived easily. In fact, direct isotopic measurements 
in the Sun can only be performed through the analysis of CO vibration-rotation lines in 
the infrared spectrum. Therefore, except for H and the noble gases, isotopic abundances 
of the Solar System are directly obtained from terrestrial samples. Other sources are 
used for estimating isotopic ratios of the most volatile elements (e.g., ^He/'^He), such as 
the solar wind or lunar samples. In other cases, such as deuterium or ^^N/^^N, estimates 
rely on measurements of the Jovian atmosphere. 

It is worth noting that independent estimates of the chemical composition of the 
Sun can also be obtained by comparison with solar-like stars, OB stars (hot, short-lived, 
massive main sequence stars), H II regions (large clouds of gas and plasma, characterized 
by the presence of ionized H, with recent star forming episodes, such as the Orion 
nebula), or with the interstellar medium in the solar neighborhood. Once the effects 
of diffusion in the Sun and the chemical enrichment of the Galaxy during the past 4.6 
Gyr are taken into account, very good agreement is achieved. Many compilations of the 
Solar System abundances have been published and the interested reader is referred to 
Refs. [ilHElElIIO]. 

Starting in the 1960s, studies have revealed that a minor fraction of the material 
embedded in meteorites (some micron-sized and smaller dust particles) showed huge 
isotopic anomalies when compared with standard Solar System material (see Fig. 
[2]). This has been interpreted as heterogeneities present already in the Solar nebula 
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Figure 2. Oxygen isotopic (number) abundance ratios measured in presolar oxide and 
silicate grains. The ellipse in the inset covers the range of ^^'^*0/^^0 ratios, on a linear 
scale, measured in almost all Solar System samples (dashed lines indicate terrestrial 
ratios) . Presolar grains make up at most a few hundred parts per million (ppm) of the 
mass in the most primitive meteorites. Notice that their oxygen isotopic ratios deviate 
by orders of magnitude from what is found in bulk solar system matter. Grain data 
are based on Nittler et al. (2008) [11], Gyngard et al. (2010) [12], and Nguyen et al. 
(2010) [T^, and references therein. The terrestrial and Solar System values are based 
on Clayton (1993) [H] and Lodders (2003) ^. Figure courtesy of L. Nittler. 



before the first solids began to condense and hence have been coined presolar grains. 
Unfortunately, the tiny amount of radioactive material present in those grains does not 
allow us at present to directly confirm this conjecture through radioactive dating. These 
grains are presumably produced by ancient stars of previous generations. So far, silicon 
carbide (SiC), graphite (C), diamond (C), sihcon nitride (Si3N4), silicates, and oxides 
(such as corundum and spinel) have been identified as presolar grains. Ion microprobe 
analyses of single presolar grains reveal a variety of isotopic signatures that allow, 
via comparison with predictions from stellar models, an identification of their parent 
stars. Low-mass red giants and asymptotic giant branch (AGB) stars (Sec. 5) have 
been identified as stellar sources of the majority of SiC and oxide grains. Supernovae 
(most likely of the core-collapse variety; see Sect. 6) are likely the sources of a small 
fraction of the SiC and Si grains, of low-density graphite grains, and of a few corundum 
grains. Moreover, nanodiamonds carry a supernova signature in their Xe and Te isotopic 
abundances. Finally, classical novae have been claimed as the likely stellar parents of a 
small number of SiC, graphite and oxide grains. There are many excellent reviews on 
presolar meteoritic grains, and the interested reader will find more information in Refs. 

dsi [mini [TBI [19]. 
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3. Big Bang nucleosynthesis 

3.1. Early universe. Microwave background radiation 

Initial ideas for the Big Bang model were proposed by Gamow, Alpher, Bethe, Herman 
and collaborators in the 1940s in order to explain the origin of the chemical elements. 
They proposed that the Universe was initially very dense and hot, and that it expanded 
and cooled to its present state. The elements would then be synthesized during an 
early time when the temperature and density conditions were appropriate for nuclear 
reactions to occur |20l [21] . Such a model would also predict a relic background radiation 
of photons with a present temperature of a few kelvin [22] . While the bulk of the elements 
are produced in stars rather than the early Universe, the idea proved correct for the 
origin of the light species ^H, ^He, ^He and ^Li. In addition, the observation of the 
cosmic microwave background radiation |23], corresponding to a blackbody spectrum 
with a temperature of 2.73 K, was of paramount importance in this regard. It singled 
out the Big Bang theory as the prime candidate for the model of our Universe. 

Modern theories of cosmology are based on the assumption of a homogeneous 
and isotropic Universe, as implied by the Cosmological Principle. The geometry and 
evolution of the Universe can then be predicted by the theory of general relativity, 
with a number of cosmological parameters describing the spatial curvature and overall 
expansion. The Hubble parameter, H, provides a measure for the expansion rate of the 
Universe and its present-day value is called Hubble's constant, Hq. The mass density of 
baryons (more precisely, of nucleons), pb, can be expressed relative to the critical density, 
Pc, by introducing the cosmological baryon density parameter, flh = Pb/Pc- The critical 
density describes the borderline between a closed and an open Universe, that is, the total 
density at which the Universe is spatially fiat, and is defined by the Friedmann equation 
as Pc = 3Hq/87tG, with G the gravitational constant. The number density of photons 
is expected to have been constant since the epoch of electron-positron annihilation, 
which occurred during a time of ~4-200 s after expansion began. The photon number 
density after that epoch can be found from the precise value for the present-day 
temperature of the cosmic microwave background radiation, T = 2.725 ± 0.001 K |24j . 
and is given by n^ = 410.5 cm~^. Using the definitions r] = Ub/n^, rjiQ = rjlO^^ and 
h = i^o/lOO km s~^ Mpc~^, one finds for the relationship of the parameters introduced 
above the expression fib h? = rjio/273.9. 

The microwave background radiation carries a record of the conditions in the early 
Universe at a time of last scattering, when hydrogen and helium nuclei recombined 
with electrons to form neutral atoms. As a result, photons decoupled from baryons 
and the Universe became transparent to radiation. Any oscillations in the photon- 
baryon fluid around that time (~400,000 years after the beginning of expansion) would 
give rise to tiny variations of temperature (on the 10~^ level) in different parts of the 
present microwave sky. The precise mapping of these anisotropics by NASA's Wilkinson 
Microwave Anisotropy Probe (WMAP) represented a spectacular success for cosmology. 
The observed anisotropics are usually decomposed in terms of spherical harmonics. 
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where each term describes the magnitude of the anisotropy on a particular angular 
scale. The observed features in the resulting angular power spectrum are closely related 
to specific cosmological parameters. The analysis of the 5-year WMAP data [251 EH] , 
when combined with other constraints [28| 129] . yields values oi h = 0.705 ± 0.013 for 
the dimensionless Hubble parameter and Qb h"^ = 0.0227 ± 0.0006 for the physical 
baryon density. From these results one can derive a number of other parameters: 
to = 13.72 ±0.12 Gyr for the age of the Universe, p^ = (9.3 ±0.2) x 10-^° g cm'^ for the 
critical density, ^b = 0.044 ±0.003 for the baryon density, and rj = (6.2 ±0.2) x 10"^° for 
the baryon-to-photon ratio. The WMAP data also provide estimates of Qc = 0.23 ±0.01 
for the dark matter density and ^a = 0.73 ± 0.02 for the dark energy density. These 
imply, respectively, that ordinary (baryonic) matter makes up only ~20% of all matter 
and that the expansion of the Universe is presently accelerating. 

3.2. Nucleosynthesis. The lithium problem 

Besides the cosmic microwave background, the other relic of the Big Bang is the 
abundance distribution of the light nuclides ^H, ^He, ^He and ^Li. Obviously, it is of 
utmost importance for any useful cosmological model that these two signatures provide 
a consistent description. When the Universe was less than 0.5 s old, at temperatures of 
T > 15 GK, the energy density was dominated by radiation (photons and neutrinos). 
The process e~ + e~^ i-^ u + i? equilibrated the electron and neutrino gases, while the 
weak interactions e~ + p ^ Ue + n, e+ ± n -H- Ve + V and n -^ p + e^ + v^ coupled 
the electron and neutrino gases to the baryon gas. In equilibrium the neutron-to- 
proton number ratio is determined by a Boltzmann distribution, n^/np = e"'^/'^^, 
where Q = 1293.3 keV [27] denotes the neutron-to-proton mass difference. During 
expansion and cooling, a temperature is eventually reached where the neutrino (weak) 
interaction processes become too slow to maintain equilibrium. The freeze-out of the 
weak interactions depends on their cross sections and occurs near T ~15 GK, at a time 
of ~0.5 s, when the neutron-to-proton number ratio is near 2/5. Beyond a time of ~10 s, 
when T ~3 GK, the decay of free neutrons to protons, with a half-life of T1/2 = 614 s 
|30j . becomes the dominant weak interaction. As will be seen below, further expansion 
and cooling gave rise to the onset of primordial nucleosynthesis. This stage is reached at 
a temperature and density of T ~0.9 GK and p^ ~ 2 x 10~^ g cm~^, respectively, near 
a time of ~200 s, when nuclear reactions can compete with the destruction of nuclei by 
photons from the high-energy tail of the Planck distribution. By this time, the decay 
of free neutrons gave rise to a neutron-to-proton number ratio of about 1/7. 

The subsequent nuclear reactions are relatively fast and, for reasons given below, 
nearly all neutrons are incorporated into the tightly bound species ^He, while only 
very small amounts of other nuclides are synthesized. Under such conditions, the 
primordial ^He abundance can be estimated by using a simple counting argument: 
for a ratio Un/up = 1/7, two out of each eight nucleons will end up bound in ^He. 
Consequently, we find for the predicted primordial helium mass fraction a value of 
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XP^f^d, ^ 2nn/{np + nn) = 2/8 = 0.25. From the above arguments it is apparent that the 
primordial ^He abundance is determined by the weak interaction cross sections (which 
are normahzed by the neutron half-hfe), the neutron-proton mass difference, and the 
expansion rate, but is rather insensitive to either the baryon density or to any nuclear 
reaction cross sections. Observations of ^He in metal-poor clouds of ionized hydrogen (H 
II) in dwarf galaxies reveal that there is a small contribution from stellar nucleosynthesis 
and that it is correlated with metallicity. Extrapolation to zero metal abundance yields 
an observed primordial ^He mass fraction of X°^'^ = 0.248 ± 0.003 [31]. The agreement 
with the predicted value can be regarded as a key piece of evidence for the standard 
cosmological model. 

For a time window of ~200-1000 s after expansion began, corresponding to 
temperatures and densities of T ~0.9-0.4 GK and p?, ~ 2 x 10~^ — 2 x 10~^ g cm~^, 
respectively, the early Universe passed through an epoch of nucleosynthesis. Note that 
the cosmic microwave background radiation and primordial nucleosynthesis probe very 
different eras of the cosmic expansion. Initially, at the higher temperature end, the 
strong and electromagnetic interactions are sufficiently fast to ensure (quasi-static) 
equilibrium among the abundances of the light nuclides ^H, ^H, ^He, ^He, ^Li and 
^Be. As the Universe expands and cools, the nuclear reactions slow down, both because 
there is a decrease in the density and because the Coulomb barriers become harder to 
overcome. As a result, individual reactions freeze out of equilibrium at characteristic 
values of temperature. The final abundance of a particular nuclide in primordial 
nucleosynthesis is then mainly given by the rate ratio of the largest production and 
destruction reactions at the freeze-out temperature [321 [33]. Nuclear reactions cease 
completely once temperature and density are sufficiently low. Among the synthesized 
light species, ^He becomes by far the most abundant one, since it has a much higher 
binding energy per nucleon {B/A = 7.074 MeV [31]) compared to all other nuclides in 
this mass region. 

Primordial nucleosynthesis starts with the p(n,7)d reactioii||, resulting in a rapid 
increase of the deuterium abundance followed by destruction via d(d,n)^He, d(d,p)t and 
t(d,n)^He. Tritium is produced by the d(d,p)t and ^He(n,p)t reactions, and mainly 
destroyed via t(d,n)^He. The most important production and destruction reactions for 
^He are t(d,n)^He and ^He(t,7)^Li, respectively. The species ^Li is mainly produced 
as ^Be via '^He(a;,7)^Be, with ^Be(n,p)^Li and ^Li(p,Q;)'^He as the most important 
destruction mechanisms. For most of the reactions, direct cross section measurements 
exist at relevant energies corresponding to the individual freeze-out temperatures |35j . 
Only for the p(n,7)d reaction are the rates largely based on theory [136], although the 
rate uncertainty has been estimated to be on the order of only ~1%. Precise primordial 
abundances are computed numerically using reaction network codes [37], incorporating 

I The nuclear physics notation for interactions is as follows: p(n,7)d stands for the nuclear reaction 
p + n ^)- J + d, where 7 denotes a 7-ray; ^^N{(3'^i/)^^C denotes the /3 (positron)-decay of ^^N to the 
daughter ^"^C, where i' denotes a neutrino. The symbols p, n, d, t and a stand for protons, neutrons, 
deuterons (^H), tritons (■^H) and a-particles (^He), respectively. 
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Figure 3. Big Bang nucleosyntliesis computed using the WMAP value of the baryon- 
to-photon ratio, ?7io*^'^^ — 6.2±0.2. (Left panel) Mass fractions of the most important 
light species versus time. The hydrogen abundance (between 0.75 and 0.83 by mass) 
is almost constant on this scale and is not displayed. The neutron abundance (dashed 
line) declines even at late times because of radioactive decay. Note that the species 
^Li is mainly produced as radioactive ^Be. (Right panel) Evolution of temperature 
and density. Significant nucleosynthesis starts at t ~200 s, when T ^0.9 GK and 
p ^ 2x 10^^ g cm^'^. For comparison, the density of air at room temperature amounts 
to Pair ^ 1-2 X 10^"^ g cm^'^. Data are courtesy of A. Coc. 



all important production and destruction reactions. The rates of nuclear reactions 
depend on the baryon density, or the baryon-to-photon ratio rjiQ, which is the only free 
parameter for Big Bang nucleosynthesis in the standard cosmological model. If one 
adopts the WMAP value, rj^^^^ = 6.2 ± 0.2, in the nucleosynthesis calculations, then 
standard primordial nucleosynthesis becomes a parameter-free model. A comparison 
of computed final abundances to observed values can then be used to investigate the 
galactic chemical evolution of the light species or to search for possible extensions of 
the standard cosmological model. The results of such a reaction network calculation are 
displayed in Fig. [31 

The primordial abundance of ^He, computed using rj^^^^ and the mean neutron 
lifetime from Ref. [30], amounts to X^'^'^ = 0.2486 ± 0.0002 [38]. The result agrees 
well with the observational value quoted above. For the number abundance ratio 
of deuterium to hydrogen, the predicted value using rj^^^^ has been estimated to 
[D/Hy^'^'^ = (2.5 ±0.2) X 10~^ [38]. Deuterium is observed via isotope-shifted (Lyman- 
a) absorption lines arising from low-metallicity gas clouds that fall on the line of 
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sight between the observer and a high-redshift quasar. The low metalhcity imphes 
that the processing of deuterium by previous-generation stars, which always deplete 
deuterium during hydrogen burning p9], is negligible. Averaging the existing seven 
high-quahty measurements yields a value of (D/H)"''^ = (2.8 ± 0.2) x 10~^ |iO] . 
Considering that the computed deuterium abundance depends rather strongly on rjio, 
the agreement between predicted and observed values of D/H can be regarded as 
another key piece of evidence in favor of the standard cosmological model. Finally, 
the predicted number abundance ratio of ^Li to hydrogen at rj^'^^^^ amounts to 
{^Li/Hy^'^ = (5.2 ± 0.7) X 10"^° [38j. It has been known for some time that 
unevolved metal-poor dwarf stars exhibit a remarkably constant lithium abundance 
[H], independent of metalhcity. The "Spite plateau" has been interpreted as being 
representative for the primordial lithium abundance. Recent work has demonstrated 
that, at the lowest observed metallicities, there is a sizable scatter in lithium abundance 
and that the lithium abundance is correlated with metalhcity [Hj. However, the 
measurement of the primordial ^Li abundance remains a challenging problem. One 
source of systematic uncertainty is the determination of the effective temperature of the 
stellar atmosphere in which the lithium absorption line is formed. Another problem is 
the possible depletion of lithium, either before the currently observed stars were formed 
or within the stars we currently observe. Most recent studies of metal-poor field and 
globular cluster stars yield estimates for the primordial ^Li abundance in the range of 
CLi/Hy^ ~ 1.1 - 2.3 X 10-1° [321 SSI m], where the spread in values indicates the 
current systematic uncertainty. Despite these observational difficulties, the predicted 
and observed values are clearly in disagreement (at the 4 — 5a level). This "lithium 
problem" represents the central unresolved issue in primordial nucleosynthesis. 

At present it seems unlikely that this problem is caused by either erroneous reaction 
cross sections or by wrong effective temperatures. Recent work has focused on two 
possibilities: (i) depletion of ^Li in stellar interiors, and (ii) physics beyond the standard 
cosmological model. The first possibility has been discussed in the literature for a 
long time and many different mechanisms have been proposed, such as atomic and 
turbulent diffusion processes, meridional circulation, gravity waves, or rotational mixing 
|45j. The second possibility is more intriguing, and may involve a time variation of 
the fundamental coupling constants [Ml ST] , modification of the expansion rate during 
Big Bang nucleosynthesis [3H1 HI] , neutrino degeneracy [50] , or negatively charged relic 
particles [51]. Clearly, the future resolution of the lithium problem promises to be an 
exciting intellectual endeavor. 

4. The Sun and low-mass stars 

Stars form in large interstellar clouds of dim, cold gas by gravitational contraction, in 
timescales ranging between ~0.01 - 100 Myr [52]. Once a critical ignition temperature 
of ~ 10^ K is achieved, protostars begin to fuse deuterium. Stars reach equilibrium 
configurations in which the gravitational pull is balanced by the internal (thermal) 
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pressure. But unlike planets, stars achieve equilibrium when their central temperatures 
are high enough to induce thermonuclear fusion reactions, which become their main 
source of energy. The minimum mass required to form a star is about 0.08 M©. 
Protostars with masses below that threshold value are called brown dwarfs, and never 
reach sufficient temperatures to initiate hydrogen fusion. Nevertheless, brown dwarfs 
heavier than ~13 times the mass of Jupiter do fuse deuterium during a short period. 
The first thermonuclear fusion stage that all stars undergo is central H-burning (T ~ 10^ 
K), by far, the longest and most stable phase in the life of a star. The reason is twofold: 
first, hydrogen is the most abundant element in the Universe (and hence, in stars); and 
second, H has a small Coulomb barrier, which favors nuclear fusion even at moderately 
low temperatures. Astrophysicists refer to this initial stage as the main sequence, and its 
duration relies almost exclusively on the overall stellar mass. A simple estimate of the 
time spent by a star on the main sequence can be derived from a dimensional analysis: 
the energy emitted by a star, E, in a time t, can be expressed as L = E/t, where L 
is the star's luminosity. Since E is mainly provided by nuclear fusion reactions, it can 
be expressed as E = fMc^, with / being the fraction of the star converted into energy 
through H-burning. This yields Ims = fMc^/L oc M/ L. But the relationship between 
the star's luminosity and its mass is far from being simple. In contrast to what is usually 
found in many stellar astrophysics textbooks, no single power law for a mass-luminosity 
relationship fits the entire main sequence. Empirically, one finds L oc M", with a in the 
range of 1.9-4.8 |53]. Taking this expression into account, and using a mean value of 
a = 3.5, one finds Ims oc M/L oc M/M^-^ oc M~^-^. Consequently, the larger the mass 
of a star, the faster it evolves and eventually dies. Indeed, the duration of the main 
sequence phase decreases dramatically for massive stars. For instance, a star like the 
Sun spends about 10 Gyr at this stage, fusing ~600 million tons of hydrogen to helium 
per second. A star with only a tenth of a solar mass would require about 3 trillion years 
to consume its hydrogen fuel, a time that is much longer than the current age of the 
Universe. In sharp contrast, a star of 100 Mq would invest only 200,000 yr on the main 
sequence. Moreover, since radioactive dating yields a value of 4.6 Gyr for the age of 
the Solar System, it is expected that the Sun will spend another ~ 5 Gyr on the main 
sequence. 

One of the natural consequences of H-burning is the fact that the central layers 
of a star become progressively enriched with He. Actually, once the central hydrogen 
content is exhausted, the He-rich region can again gravitationally contract toward a new 
equilibrium configuration. Depending on its mass, a star can either become a stable He 
white dwarf or reach sufficiently high temperatures to fuse helium to carbon (Sect. 5). 
This requires temperatures near 10^ K because of the larger Coulomb barrier associated 
with the "^He nucleus. As a result, the star is again stabilized, winning its second match 
against gravity. It is also worth noting that the overall duration of the He-burning phase, 
about 100 Myr for a 1 Mq star [51], is shorter that the main sequence stage. This is 
because helium fusion produces less energy per gram of fuel than hydrogen fusion. 

Helium burning results in a carbon- and oxygen-rich core. Central He exhaustion 
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forces the stellar CO-ricli core to contract, rising once more its temperature. Stars with 
masses up to 8 Mq are not able to achieve sufficient temperatures to ignite another 
burning episode. For instance, carbon burning requires about 0.6 — 1 GK. Hence, such 
stars end their lives as CO white dwarfs, representing stellar remnants of planetary size 
that are supported by the pressure of degenerate electrons. 

4-1- Hydrogen burning: proton-proton chains vs. CNO-cycles 

During H-burning, four protons are transformed into a single "^He nucleus, with an 
energy release of 26 MeV per process. This is accompanied by the emission of two 
positrons and two neutrinos. Two major modes for H-burning are possible in principle. 
The ffist one is called proton-proton (or pp) chains and it is the dominant mode at 
low temperatures or in low metallicity environments. This burning regime is initiated 
by the reaction p + p— )-(i + e^ + z/, where e^ denotes a positron. Alternatively, two 
protons may fuse via the pep-reaction, p-|-p-|-e~— J-cZ-l-z/. The latter reaction depends 
on the electron density and thus can only compete with the former for densities of 
> 10^ g cm~^, which are usually not achieved during H-burning in low-mass stars. 
Subsequently, deuterium and hydrogen fuse via the reaction d(p,7)^He. This can 
be followed either by 3He(3He,2p)^He (ppl chain), 3He(^He,7)^Be(e-,z/)^Li(p,^He)^He 
(pp2 chain), or ^He(^He,7)'^Be(p,7)^B(/9+z/)^Be— )-^He-|-^He (pp3 chain). The relevance 
of each particular chain depends critically on the temperature as well as on the ^He 
content. In the Sun the central temperature amounts to 15.6 MK [5H] and about 90% 
of the nuclear energy is produced via the ppl chain. 

The second H-burning regime are the CNO-cycles. They require somewhat higher 
temperatures than the pp-chains, as well as the presence of suitable amounts of CNO- 
nuclei, which act as catalyzers. Because of the larger temperature dependence of the 
energy production rate through the CNO-cycle, as compared with that of the pp- 
chains, the CNO-cycle becomes the dominant mode above ~20 MK. Hence, the pp- 
chains dominate the energy production during H-burning in low-mass stars, whereas 
the CNO-cycles become the chief regime in massive stars, which are characterized by a 
higher central temperature. Since ^^C is the fourth most abundant nuclide and proton 
captures destroy ^^C much faster than ^^O (which is the third most abundant nuclide, 
see Fig. 1), the ^^C(p,7)^^N reaction often triggers the CNO-cycle. This is followed 
by ^3n(/3+i/)13C(p,7)^^N(p,7)150(/3+z/)15n(p^«)12c (CNOl or CN cycle). As was the 
case for the pp chains, four protons are transformed into a single helium nucleus, while 
two positrons and two neutrinos are emitted. The process ^^N(p,7)^^0 has by far the 
smallest reaction rate among the reactions that form the CNOl cycle and, therefore, 
the major fraction of the initial ^^C will be converted to ^^N at the end of H-burning. 
Alternative cycles are ^^N{p,-fy^O{/3+uy^N{p,-fy^O{p,-fy'^F{/3+uy'^0{p,ay^N (CN02 
or NO cycle), ^^n{p,-fy^O{p,-fy^F{P+iyy^O{p,-fy^F{P+uy^O{p,ay^N (CN03 cycle), 
and ^60(p,7)i^F(/3+z/)i70(p,7)^8F(/3+z/)i80(p,7)^9F(p,a)i60 (CN04 cycle). 

For a long time it was difficult to prove that energy production inside the Sun 
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Figure 4. Solar neutrino energy spectrum computed for a specific solar model |69) . 
Neutrino fluxes (in cm^^ s~^) correspond to the predicted values on Earth. In the 
figure, pep denotes the reaction p + p + e^ — ^ d + v, while hep corresponds to 
■^He(p,e"'"^)^He. The uncertainties associated with each neutrino signal are indicated. 
Figure from Bahcall, Serenelli & Basu (2005) [70], reproduced by permission of the 
American Astronomical Society. 



and other stars is achieved by nuclear fusion reactions. Direct evidence was obtained 
by detecting the neutrinos emitted during H-burning in the solar interior using an 
underground laboratory located in the Homestake gold mine, South Dakota (USA) 
[521 ED]- Interestingly, the number of detected solar neutrinos disagreed with the 
value expected from the standard solar model. A number of explanations were put 
forward to solve this solar neutrino problem. They included uncertainties in nuclear 
reaction rates, problems with the standard solar model itself, as well as new physics 
associated with neutrino properties [61] . A number of experiments were conducted to 
solve the controversy, including studies of the solar interior through helioseismology, 
improvements of the standard solar model, nuclear physics experiments, and extensive 
solar neutrino measurements at new facilities, such as GALLEX [62], SAGE [63], Super- 
Kamiokande [M], SNO [651 |66l [67] , and Borexino [68]. 

Models predict that a suite of neutrinos, with different fluxes and energies (Fig. H]), 
are produced in the solar interior [611 EO]- The dominant source is the p(p, e"*" z/)d 
reaction, generating neutrinos with energies below ~0.4 MeV. These can be detected 
by both the GALLEX and SAGE gallium detectors. Recent theoretical estimates [71] 
of the solar pp neutrino flux yield a value of 6 x 10^° cm~^ s~^ with an associated 
uncertainty of < 1%. A much smaller flux, 1.4 x 10® cm~^ s~^, is expected for neutrinos 
produced in the pep reaction with an energy of 1.4 MeV, which are detectable by liquid 
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scintillators (e.g., Borexino), or chlorine and water detectors. Neutrinos produced by 
^Be(e~,z/)^Li have discrete energies of 0.38 and 0.86 MeV, with a predicted flux of 
4.6 X 10^ cm~^ s~^ (±11%). These neutrinos should contribute substantially to the 
fluxes measured in gallium and chlorine detectors. The flux expected from ^B(/3+z/)^Be 
neutrinos is very small, about 4.9 x 10® cm~^ s~^. However, because of their high energy 
(up to ~ 15 MeV), the ^B neutrinos dominate the fluxes of many chlorine and water 
detectors, such as the pioneering Homestake mine experiment, or the more recent SNO 
and Super-Kamiokande experiments. Unfortunately, the flux uncertainty associated 
with these neutrinos (±16%) is rather large. Neutrinos produced in other reactions, such 
as 3He(p,e+z/)^He (hep), ^^N{l3+i^y^C, ^^0{l3+iyy^N, or i^F(/3+z/)i^0, have estimated 
fluxes of 8.2 X 10^, 2.1 x 10*, 1.5 x 10®, and 3.5 x 10® cm~2 g-^ respectively. So far, the 
solar *B neutrino flux has been directly measured by the SNO collaboration [66l [72| [73] 
and, more recently, by Borexino [68]. Moreover, the ^Be neutrino flux has also been 
measured by Borexino, with a 10% uncertainty [71]. Efforts to measure the neutrino 
fluxes associated with the /9-decays of ^^N and ^^O are currently ongoing at Borexino, 
Super-Kamiokande and the upcoming SNO±. 

One of the difficulties in solving the solar neutrino problem was related to the fact 
that the existing neutrino detectors were sensitive only to electron neutrinos. Indeed, 
the detected number of electron neutrinos was only a third of the predicted number. 
A clue came from Super-Kamiokande measurements in 1998, first suggesting a possible 
transformation of muon neutrinos, created in the upper atmosphere by energetic cosmic 
rays, to tau neutrinos. The first direct evidence of solar neutrino oscillations was found 
at SNO [67] in 2001. Solar neutrinos from ®B-decay were detected via charged current 
reactions on deuterium (sensitive to electron neutrinos only) and via elastic scattering 
on electrons (sensitive to all neutrino fiavors, i.e., electron, muon and tau neutrinos). 
The analysis of the data, after a thorough comparison with the elastic scattering results 
obtained at Super-Kamiokande, revealed that only ~ 35% of the solar neutrinos arriving 
at Earth are actually electron neutrinos (the number depends on the exact neutrino 
energy). With this in mind, the experimental neutrino fiux agreed with the theoretical 
predictions. In fact, the discovery of the oscillatory nature of neutrinos firmly established 
the need for new physics beyond the standard electroweak theory. 

Several nuclear physics experiments, aimed at improving estimates of solar neutrino 
fiuxes have been conducted in recent years, including measurements of ^Be(p,7)®B, 
^He(a,7)^Be, and ^^N(p,7)^^0. Indeed, the rate reduction by a factor of ~2 inferred for 
the ^^N(p,7)^^0 reaction after new measurements performed underground at LUNA 
[75| [76] and above ground at LENA [77], had important implications for several 
astrophysical scenarios [78j. First, it reduces the expected fiux of solar CNO neutrinos 
by a factor of 2. Second, it affects the time at which stars leave the main sequence 
after H-burning (so-called turnoff point). Indeed, the reduction of the ^'^N(p,7)^^0 
rate translates into a smaller age for a given turnoff luminosity. This results in an 
effective reduction of globular cluster ages by about 1 Gyr. Third, it also provides 
better agreement between stellar models and observations of AGB stars regarding the 
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chemical abundance pattern for species heavier than carbon. 

4-2. The standard solar model: current status 

The agreement between theoretical predictions and observations regarding the solar 
neutrino fluxes, achieved after the discovery that neutrinos change flavor, gives 
credibility to the existing solar models. This discovery spawned a frantic activity aimed 
at improving the input physics adopted in the solar models. Solar abundances received 
particular attention. Early measurements of the photospheric solar oxygen content [79] 
paved the road for a major revision of the solar abundances that resulted in a dramatic 
reduction of the overall solar metallicity (in particular, for CNO and Ne) with respect 
to previous estimates [7] . Interestingly, the new solar models computed with the revised 
solar abundances showed severe discrepancies for the properties of the internal solar 
structure when compared with helioseismology determinations. Ironically, shortly after 
the solar neutrino problem was solved, a new and yet unsolved solar abundance problem 
appeared. The reader is referred to Refs. [Ml [801 EB [82] for a thorough analysis of the 
impact of the new solar abundances on the properties of the solar interior. 

Standard solar models rely on the evolution of a 1 Mq star from the pre-main 
sequence up to the current Solar System age of 4.57 Gyr. Models are constrained by 
the observed values for the present solar radius and luminosity (6.96 x 10^° cm and 
3.84 X 10^^ erg s~^, respectively). In addition, they are required to match the observed 
metal-to-hydrogen ratio in the solar photosphere. Emphasis is placed on four important 
quantities that can be compared directly with helioseismology data: the surface helium 
mass fraction, Ys, the depth of the convective envelope, Rconv, and the average relative 
differences of the sound speed and density proflles. Whereas improved choices for the 
equation of state and for a number of nuclear reaction cross sections have shown little 
impact on the overall properties of the models, the adoption of different abundance sets 
translated into signiflcant differences [71]. For example, the abundances of Grevesse & 
Sauval [To], which were widely used in the past, yield values of Ys = 0.2423 and Rconv 
= 0.713 Kq, while models computed with the improved abundances from Asplund et al. 
[7| give values of Ys = 0.2292 and Rconv = 0.728 R©. For comparison, helioseismology 
yields Ys = 0.2485 ± 0.0035 [80], and Rconv = (0.713 ± 0.001) R© p]. This clearly 
stresses that recent revisions of the solar abundances result in larger disagreements with 
helioseismology data. Moreover, the most recent reanalysis of the solar abundances by 
Asplund et al. [8] , which resulted in slightly higher CNO and Ne abundances compared 
to their previous determinations [7], does still not restore the agreement between solar 
models {Ys = 0.2314 and Rconv = 0.724 Rq) and helioseismology. 

Similar discrepancies are obtained between theoretical predictions from solar 
models with improved abundances and the corresponding measurements of a suite of 
physical properties of the deepest layers of the Sun (probed through helioseismology 
measurements with high penetrative, low-degree modes), of neutrino fluxes, or even of 
the Solar System age. In all cases, models computed with higher solar metallicities (as 
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in Grevesse & Sauval [TU]) result in much better agreement with hehoseismology than 
those computed with lower metallicity values (as in Asplund et al. [HE]). In particular, 
the former yield an age consistent with current estimates for the Solar System (~ 4.6 
Gyr), whereas the latter suggest a value between 4.8-4.9 Gyr plj. 

Titanic efforts have been invested in studying the suitable changes in the physics 
adopted for solar models in order to achieve better agreement with results from 
hehoseismology measurements. Particular attention deserves the 10,000 Monte Carlo 
simulations of solar models with 21 input parameters that are randomly varied within 
expected uncertainties [82]. Many different solutions, but none satisfactory, have 
been proposed. The list includes changes in the solar composition (in particular, an 
enhancement of the neon abundances [HI [85], [86]), enhanced microscopic diffusion 
[80l [87l [88] , accretion of metal-poor material by the Sun jHHl [89] , and changes in the 
radiative opacities by 15% [HI [HSl [871 [90]. 

5. Red giants and AGB stars. The s-process 

5.1. Evolution of massive stars and AGB stars 

A typical massive star of solar metallicity (M > IOMq) spends about 10 My on the main 
sequence, fusing hydrogen to helium via the CNO cycles in its core. When the hydrogen 
in the center is exhausted, hydrogen burning continues in a shell. The core contracts and 
heats up until helium is ignited. This new source of nuclear energy heats the overlying 
hydrogen shell and the outer layers of the star expand greatly. The star shows up in 
the Hertzsprung-Russell diagram at the highest observed luminosities and becomes a 
supergiant. Examples are Rigel (blue supergiant) and Betelgeuse (red supergiant) in 
the constellation Orion. Core helium burning lasts for about 1 Myr. When helium 
is exhausted in the center, helium burning continues in a shell located beneath the 
hydrogen burning shell. Such massive stars are capable of initiating advanced burning 
stages in their cores, which are referred to as carbon burning, neon burning, oxygen 
burning, and silicon burning [39]. Each time a given fuel is exhausted in the core, the 
very same fuel continues to burn in a shell just outside the core. At the same time, the 
core contracts, thereby raising the temperature and pressure in order to stabilize the 
star. Eventually, the ashes from the previous core burning stage begin to fuse, providing 
a new source of nuclear energy and thus halting temporarily the contraction of the stellar 
core. The duration of each subsequent advanced burning stage decreases significantly. 
For example, while hydrogen burning may last several million years, core silicon burning 
may last only one day. The reasons are twofold. First, hydrogen burning releases far 
more energy per unit mass compared to subsequent burning stages and, hence, the 
hydrogen fuel is consumed much slower. Second, the manner by which the star radiates 
energy (or "cools" ) changes dramatically. For hydrogen and helium burning, the nuclear 
energy generated in the core is carried by photons, which scatter many times on their way 
to the surface before being emitted from the star. Beyond helium burning, most of the 
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stellar energy is radiated as neutrino-antineutrino pairs, produced via electron-positron 
pair annihilation (e^ + e^ — )■ z/ + z/) or the photo-neutrino process (7 + 6" — > e~ + z/ + z/). 
Since the neutrino energy losses increase dramatically during the advanced burning 
stages and the neutrinos escape freely from the star, the evolution of the massive star 
accelerates rapidly. At the end of silicon burning the star consists of several layers 
of different composition, separated by thin nuclear burning shells, surrounding a core 
mainly composed of ^^Fe (Sec. [6]). We do not describe the advanced burning stages in 
detail here since a modern discussion can be found in standard texts (see, for example, 
Ref. [39]). However, note that the nucleosynthesis in terms of nuclear processes and 
final abundances is very similar to explosive nucleosynthesis in massive stars, which is 
discussed in Sec. 6.3. As will become clear below, the explosive burning stages operate 
at higher temperatures and also occur in reversed order compared to the advanced 
burning stages in hydrostatic equilibrium. For more information on the evolution of 
massive stars, see Sec. |6]and Refs. [9T| [92]. 

Low-mass {OAMq < M < 2Mq) and intermediate-mass {2Mq < M < IOMq) 
stars of solar metallicity evolve very differently from their massive counterparts. After 
hydrogen is exhausted in the core, hydrogen burning continues via the CNO cycles in 
a shell and the star leaves the main sequence. The core contracts and the hydrogen 
burning shell becomes hotter. As a result, the surface expands dramatically and the 
star ascends the red giant branch. The envelope becomes convective, grows in size, and 
eventually dredges up matter to the surface that had experienced CNO processing (first 
dredge-up). When the temperature in the contracting core becomes sufficiently high, 
helium burning is initiated and the star now occupies the horizontal branch. Eventually 
the helium in the core is exhausted, the core contracts again, and the star ascends 
the asymptotic giant branch (AGB). At this time, stars with masses of M > 4:Mq 
experience a second dredge- up event, which carries products of shell hydrogen and core 
helium burning to the stellar surface. Eventually, helium burns in a shell. The star now 
burns nuclear fuel in two shells: helium in a region surrounding a degenerate carbon- 
oxygen core, and hydrogen in a region surrounding the helium burning shell. The two 
shells are separated by a helium-rich intershell region. The details of the stellar models 
become rather complicated at this stage. The situation is represented schematically in 
Fig. |5l It is sufficient to mention here that the hydrogen and helium burning shells, 
in fact, alternate as the main contributors to the overall luminosity. Most of the time, 
quiescent shell burning of hydrogen is the main nuclear energy source, giving rise to an 
increase of temperature and density in the barely active helium rich region above the 
carbon-oxygen core. The helium eventually ignites in a thermal pulse (instability). The 
sudden energy release produced by such a pulse drives convection within the helium-rich 
intershell and extinguishes the hydrogen burning shell. The star expands and cools until 
the thermal pulse quenches. At this point, the convective envelope reaches deep into 
the intershell region, carrying the products of helium burning (mainly carbon, but also 
elements heavier than iron; see below) to the stellar surface (third dredge- up). After 
a thermal pulse, the ashes of helium burning settle onto the core and the hydrogen 
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burning shell reignites. The cycle described above may repeat many times, with a 
period on the order of about 10^ y. This evolutionary stage is called the thermally 
pulsing asymptotic giant branch (TP-AGB). During the interpulse phase, in stars with 
masses of M > 4:Mq, the base of the convective envelope reaches down to the top of 
the hydrogen burning shell, with temperatures in excess of about 50 MK (hot bottom 
burning). Because the envelope is fully convective, it is completely cycled through this 
burning region and the products of hydrogen burning nucleosynthesis are enriched at the 
stellar surface. During their evolution, TP-AGB stars experience a significant mass loss 
via strong stellar winds, thereby enriching the interstellar medium with the products of 
their nucleosynthesis. The observational counterparts belong to the spectroscopic class 
M and are called MS, S, SC or C(N) stars, depending on their spectral properties. For 
example, C(N) stars are called carbon stars since they exhibit more carbon than oxygen, 
by number, in their atmospheres (since each thermal pulse and third dredge-up event 
increases the carbon abundance). Of great importance for nuclear astrophysics was the 
discovery of spectral lines from the element technetium in certain S stars [1]. All of the 
technetium isotopes are radioactive with half-lives of less than ~ 4 x 10^ yr. Such half- 
lives are very short compared to typical lifetimes of low- and intermediate-mass stars 
(> 10^ yr). Consequently, the discovery showed beyond doubt that the technetium must 
have been produced recently within the stars and that the products of nucleosynthesis 
could indeed reach the stellar surface via dredge-up processes. Asymptotic giant branch 
stars are believed to be the main sources of carbon and nitrogen in the Universe, and 
of many heavier nuclides, as will be seen below. For more information on AGB stars, 
see Refs. M, E 



5.2. Helium burning 

Helium burning in massive stars and AGB stars takes place in the temperature range of 
T = 0.1 — 0.4 GK and starts with the fusion of two a-particles. The composite nucleus 
^Be lives for only ~ 10"^^ s and quickly decays back into two a-particles. After a given 
time a tiny equilibrium abundance of ^Be builds up, sufficient to allow for capture of a 
third a-particle to form ^^C [95]. This process is referred to as triple-a reaction. It was 
pointed out by Fred Hoyle that the process would be too slow to account for the fusion of 
^^C, unless a resonance existed right above the ^Be + a threshold that is formed without 
inhibition by the centripetal barrier. The prediction and subsequent observation [961 [97] 
of the "Hoyle state" in ^^C represented a marvelous interplay of astrophysics and nuclear 
physics. The triple-a reaction is a (sequential) three-body interaction and thus has not 
been measured in the laboratory. From experimental knowledge of the nuclear masses 
and partial widths involved, the reaction rate can be estimated. Present uncertainties 
amount to about ±15% at typical hehum burning temperatures [HHl [99], a relatively 
precise value for a process that has not be measured directly. Recently, a quantum 
three-body calculation of the triple-a reaction [100] resulted in a rate that differs at 
T < 0.25 GK by large factors (for example, by a factor of ~ 10® at T = 0.1 GK) from 
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Figure 5. Schematic representation (mass versus time) of a low- or intermediate- 
mass [M < IOMq), thermally pulsing AGB star. The figure is not to scale. Shown 
are the convective H envelope (dark grey), the radiative ^He-^^C intershell (white), 
and the degenerate C-0 core (light grey). The thick blue line indicates the H-burning 
shell, which is active between thermal pulses, and the thin blue line indicates the 
weakly active He-burning shell. For regular short periods of time the barely active 
He shell ignites in a thermal pulse (TP), giving rise to a convective region (red) that 
extends over the entire intershell, extinguishing in the process the H-burning shell. 
When the thermal pulse comes to an end, the convective shell disappears and the 
quiescent He-burning resumes. Important mixing episodes occur at the end of each 
thermal pulse: (i) the convective envelope reaches into the intershell so that synthesized 
matter is transported to the stellar surface (third dredge- up; TDU); (ii) protons diffuse 
from the base of the envelope into the intershell, where they are captured by ^^C to 
produce (after /3-decay) ^'^C. Neutrons are then released by the ^'^C(a,n)^^0 reaction, 
producing in situ the main s-process component. During the subsequent thermal pulse, 
temperatures can be high enough to initiate the ^^Ne(a,n)^^Mg neutron source. See 
text. 



the established values. The proposed significant change of this important reaction rate, 
however, was found to be in disagreement with fundamental observational evidence of 
stellar evolution |101] and energetics of type I X-ray bursts |102] . 

Clearly, a resolution of this controversy is of utmost importance. 

Helium burning continues via the ^^C(a,7)^^0 reaction. There is no resonance 
near or above the a-particle threshold in ^^O and thus this process must proceed via 
broad-resonance tails (including subthreshold resonances) and direct mechanisms. These 
amplitudes may interfere, causing problems in the extrapolation of the S-factor from 
directly measured energies (> 1 MeV in the center of mass) to the astrophysically 
important energy range (~300 keV in the center of mass), which is at present not 
accessible experimentally. Knowledge of this reaction rate is of key importance. In 
massive stars, the ratio of the triple-a and ^^C(q:,7)^^0 reaction rates determines the 
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ratio of ^^C to ^^O at the end of core helium burning. This ratio, in turn, strongly 
influences the subsequent hydrostatic burning stages, thus affecting the presupernova 
stellar structure, the explosive nucleosynthesis, and the nature of the remnant (neutron 
star or black hole) left behind after the core collapse (see below). At present the 
^^C(a,7)^^0 reaction rate at typical helium burning temperatures is uncertain by ±40% 
|103j and a more accurate rate (uncertainty of <10%) is highly desirable. For more 
information on the status of this crucial reaction, see, for example, Ref. |104j . 

A number of other important reactions occur during the helium burning 
stage. Since the stellar plasma most likely contains ^"^N (from CNO cycle 
operation during the preceding hydrogen burning stage), the reaction sequence 
^^N(a,7)^8F(/3+z/)i80(a,7)22Ne will be initiated. The subsequent ^^Ne{a,nf^Mg 
reaction takes place in massive stars and during a thermal pulse in AGB stars with 
Ma_gb ^ 4:Mq near temperatures of T ~0.3-0.4 GK, giving rise to a flux of free 
neutrons. In AGB stars with Magb ^ 4M0, after a thermal pulse and third dredge- 
up event, some protons from the convective envelope mix into the radiative intershell, 
which consists mainly of ^He and ^^C. The nature of this mixing mechanism is not 
understood at present |105j and its magnitude is usually described by a free parameter 
in models of low-mass TP-AGB stars. The protons fuse with the abundant ^^C and 
form a thin ^^C-rich pocket via the sequence ^2C(p,7)^^N(/3"''z/)^'^C near the top of the 
intershell. The subsequent ^'^C(a,n)^^0 reaction occurs during the interpulse period 
under radiative conditions at temperatures near and in excess of T ~0.09 GK, giving 
rise to a flux of free neutrons. Recent experimental work on the ^'^C(a,n)^^0 reaction 
|106] reported a rate uncertainty of ±22% for a temperature near 0.1 GK. However, 
this reaction has been directly measured only for center of mass energies in excess of 
280 keV |107j . while the astrophysically important energy range is centered near 190 
keV. For the '^'^Ne{a,n)'^^Mg reaction, the lowest lying measured resonance is located at 
700 keV in the center of mass |108] . but typical temperatures of T ~ 0.3 GK translate 
to bombarding energies near 600 keV. Recent nuclear structure measurements of the 
2^Mg compound nucleus |109l IllOl 1111] resulted in a much improved reaction rate, with 
an estimated uncertainty of ±10% near T = 0.3 GK |112j . For both of these neutron 
producing reactions, future direct measurements will be important. 

5.3. s-process 

The neutrons released by the processes described above are of fundamental interest for 
the synthesis of about half of the elements heavier than iron. The heavier elements 
cannot be produced easily in stars via charged-particle reactions, since the rapidly 
decreasing transmission through the Coulomb barrier with increasing nuclear charge 
inhibits their formation. They are made instead by exposing lighter seed nuclei to 
a source of neutrons, such that neutron capture reactions can be initiated. The 
astrophysical environments discussed above give rise to relatively small neutron densities 
(~ 10^ — 10^^ neutrons per cm^). Under such conditions, successive neutron captures 
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by a chain of isotopes occur until a radioactive isotope is reached. In most cases, the /3- 
decay probabihty of this species is much larger than the probability for another neutron 
capture to occur. Thus the radioactive nuclide will /3-decay back to stability and another 
successive chain of neutron captures is initiated. This mechanism, referred to as the 
astrophysical s(low neutron capture)-process, gives rise to a nucleosynthesis path that 
runs close to the group of stable nuclides. The abundances synthesized by the s-process 
will depend on the magnitude of the neutron-capture cross sections and the total neutron 
flux. Nuclides with very small neutron-capture cross sections, for example, those with 
magic neutron numbers (A^ = 50, 82 and 126) are expected to pile up in abundance, 
while those with large cross sections will quickly be destroyed and achieve only small 
abundances. This is the reason for the narrow peaks at the neutron magic numbers 
(corresponding to mass numbers of A ~84, 138 and 208) in the solar system abundance 
distribution. At certain locations along the s-process path, depending on the conditions, 
the abundance flow encounters unstable nuclei with decay probabilities comparable in 
magnitude to the competing neutron capture probabilities. At these locations the s- 
process path splits into two branches. The analysis of the abundance pattern of such 
branchings provides crucial constraints, not only for the physical conditions during the 
s-process (temperature, mass density and neutron density) [113] . but also for convective 
mixing (branchings at ^^^I) and the time scale (branchings at ^^Kr or ^^Zr) of a thermal 
pulse |114j . Neutron captures at the other extreme, that is, at very high neutron 
densities, will be discussed in Sec. 16.91 

Detailed analyses of the solar system s-process abundance distribution indicated 
the existence of three distinct components, referred to as the main {A ~90-205), weak 
{A < 90) and strong {A ~ 208) s-process component, which are characterized by the 
capture of about 10, 3 and 140 neutrons per seed nucleus, respectively |115j . Modern 
theories of stellar evolution have revealed the origin of these components. Low-mass TP- 
AGB stars {Magb ^ 4:Mq) are believed to be the source of the main s-process component 
[116[ 1117 ]. Note that although in this case the ^^C(a,n)^^0 reaction represents the 
predominant neutron source, the ^^Ne(Q;,n)^^Mg neutron source is marginally activated 
during a thermal pulse, thereby altering sensitively certain abundance ratios of nuclides 
belonging to s-process branchings. Low-mass, early-generation (metal-poor) TP-AGB 
stars are thought to be the source of the strong s-process component, with ^^C(a,n)^^0 
providing most of the neutrons |118l I119J . Core helium burning and subsequent carbon 
shell burning in massive stars (M > IIM0) produces the weak s-process component 
|120[ 1121] . with ^^Ne(a,n)^^Mg as the most important neutron source. 

The signatures of the s-process are imprinted not only in the solar system abundance 
distribution and in spectroscopic observations of stellar atmospheres, but also in 
meteoritic Stardust (presolar grains). The latter observations are especially important 
for constraining models of AGB stars |122j . since the anomalous isotopic ratios in these 
grains can be measured with high precision (down to uncertainties of a few percent) |19j . 
In fact, low-mass AGB stars are predicted to be the most prolific sources of dust in the 
Galaxy [123] and the majority of discovered presolar grains indeed originate from such 
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stars |124] . Analysis of the isotopic composition of heavy elements in presolar grains 
from AGB stars can constrain, for example, the impact of the ^^Ne neutron source and 
the efficiency of ^^C pocket formation |125j . Furthermore, predictions of models of dust 
production by AGB stars |126j are in overall agreement with observations concerning 
low-mass stars. However, the models also predict a large production of oxygen-rich 
(from efficient hot bottom burning) dust in intermediate-mass AGB stars, which so 
far is not supported by observation. In fact, the complete lack of oxygen-rich presolar 
Stardust grains from intermediate-mass AGB stars represents a major puzzle |1271I128] . 
For accurate modeling of s-process scenarios, the neutron capture rates for a large 
number of nuclides in the A < 210 mass region must be known, in addition to the 
rates of the neutron source reactions (see above) and the rates for neutron-consuming 
reactions (neutron poisons ^^C, ^^N, ^^O, ^^Ne and ^^Mg). It must be emphasized that 
rather precise neutron capture rates are required for s-process models (uncertainties 
of less than a few percent) at energies ranging from kT ~8 keV in low-mass AGB 
stars to kT ~90 keV during shell carbon burning in massive stars. For many of the 
important neutron-capture reactions the required level of precision has been reached and 
the reliability of the current reaction rate data set is quite impressive |129] . Nevertheless, 
additional and more accurately measured cross sections are needed for many reactions, 
including those involving radioactive branching point nuclei. In addition, nuclear theory 
is indispensible for correcting the measured laboratory rates for effects of thermal target 
excitations in the hot stellar plasma. For more information on the s-process and related 
topics, see Ref. |130] . 

6. Core-collapse supernovae (type II, lb, Ic). The p-, z/-, up-, a-, and 
r-processes 

6.1. Classification 

The different types of supernovae are classified observationally according to their spectra. 
Type II supernova spectra contain hydrogen lines, while those of type I do not. Type 
I supernovae whose spectra show no absorption caused by the presence of silicon are 
referred to as type lb or Ic supernovae. (The distinction is based on a helium line feature 
in the spectrum). Type II, lb and Ic supernovae tend to occur in the arms of spiral 
galaxies, but not in early-type galaxies. Since the spiral arms contain many massive and 
thus young stars and early-type galaxies do not contain such objects, the observations 
suggest that massive stars are the progenitors of type II, lb and Ic supernovae. 

Stellar model simulations and observations support the conclusion that these three 
supernova types result from the collapse of an iron core. Stars with initial masses 
less than M ~2O-3OM0 explode as a type II supernova and form a neutron star as a 
remnant. Stars with masses above this range (Wolf-Rayet stars), or less massive stars 
in binaries, that have lost their hydrogen envelopes are thought to be the progenitors 
of type lb and Ic supernovae. It is not clear at present if the latter explosions leave 
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a neutron star or a black hole behind as a remnant, mainly because of our incomplete 
knowledge of post-main sequence mass loss and the details of fallback of matter onto the 
central object during the explosion. As will become clear in the following, core-collapse 
supernovae are of outstanding importance for three broad topics: (i) they are predicted 
to be the most prolific sources of the elements in the Galaxy; (ii) they are the place 
where neutron stars are born; and (iii) they are a plausible source of shock waves that 
are believed to accelerate most Galactic cosmic rays. 

6.2. Evolution beyond core silicon burning 

After silicon has been exhausted in the core, a massive star (M > IOMq) consists 
of several layers of different composition ( "onion- like" structure), separated by thin 
hydrostatic nuclear burning shells, surrounding a core mainly composed of ^®Fe (Sec. 
[5]). Considering the declining number of stars of increasing mass, weighted by the fact 
that more massive stars eject a larger yield of heavy elements, one arrives at a "typical" 
model star with a mass somewhere near ^25Mq. For clarity sake, we will discuss in 
the following the evolution and explosion of such a star. It must be kept in mind that 
stars of different masses exhibit significant differences in structure, shell composition, 
location of nucleosynthesis processes, and so on. The situation is shown schematically 
in Fig. [6] (left side) for a 25Mq star of initial solar composition, displaying the most 
abundant nuclides in each layer and, at the bottom, the nuclear burning stage that 
produced those particular ashes. (The subscripts C and S stand for core and shell 
burning, respectively.) For example, the most abundant nuclide in the core is ^^Fe, 
which is the product of core and shell silicon burning. In the next layer, on top of the 
core, the species ^^Si is the most abundant nuclide, which is the product of oxygen shell 
burning. The nuclear transformation of ^^Si to ^^Fe (silicon shell burning) continues at 
the intersection of the silicon layer and the iron core, the transformation of ^^0 to ^^Si 
(oxygen shell burning) takes place at the intersection of the oxygen and silicon layers, 
and so on. Note the location of the weak s- process component (Sec. [5l), both in the 
carbon-oxygen layer (resulting from helium core burning) and the oxygen-neon layer 
(resulting from carbon shell burning). The latter layer has been predicted [131] to be a 
main source of the important 7-ray emitter ^°Fe. 

At this point, the electron degenerate core has no further source of nuclear energy to 
its disposal. It grows in mass as a result of nuclear burning in the overlying shells. When 
the core mass reaches the Chandrasekhar limit (M ~1.4Mq), the electron degeneracy 
pressure is unable to counteract gravity and, consequently, the core collapses freely at 
about a quarter of the speed of light. When the density reaches values near nuclear 
density (p ~ 10^^ g cm~^), nuclei and free nucleons begin to feel the short-range nuclear 
force, which is repulsive at very small distances. The collapsing inner core reaches 
high inward velocities, overshoots nuclear density, and rebounds as a consequence [133] . 
The rebounding part of the core encounters infalling matter, giving rise to an outward 
moving prompt shock wave. The hot and dense inner core has become a proto-neutron 
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Figure 6. Structure and evolution of a 25Mq star of solar metallicity, as predicted 
by one-dimensional, spherically symmetric models [MJ 1132] , shortly before and after 
core-collapse (not to scale). Only the main constituents in each layer are shown. 
Minor constituents, among them important 7-ray emitters, are set in thin black boxes. 
Various nucleosynthesis processes are shown in green boxes: weak s-process (s); p- 
process (p); a-rich freeze-out (arf); //-process (v). (Left) Snapshot of presupernova 
structure. Nuclear burning takes place in thin regions (burning shells) at the interface 
of different compositional layers, where each burning shell migrated outward to the 
position indicated by the blue lines. The compositions result from burning stages 
indicated at the bottom (subscripts C and S stand for core and shell burning, 
respectively). The diagonally arranged numbers indicate the interior mass (in solar 
masses) for each burning shell. (Right) Explosive nucleosynthesis resulting from 
passage of the shock wave through overlying layers, giving rise to explosive burning 
of silicon (Si^;), oxygen (O^;) and neon-carbon (Nex/Cx)- Strictly speaking, this 
classification depends on the temperature range, not on the available fuel. Nevertheless, 
the names indicate approximately which compositional layers of the pre-supernova will 
usually be affected. Outside the outer dashed line, the composition is little altered by 
the shock. The inner dashed line indicates the approximate boundary of the part of 
the star that is ejected (mass cut). See text. 



star with a mass of M ~1.5M0. While the shock moves outward through the outer 
core, it loses energy, both by photodisintegrating iron peak nuclei and by emission of 
neutrinos. About 1 s after core collapse, when the prompt shock reaches the outer edge 
of the core, it has lost its kinetic energy and stalls. How exactly the shock is revived 
and how it will ultimately propagate through the stellar layers beyond the iron core 
and disrupt the star in a core-collapse supernova is one of the most elusive problems in 
nuclear astrophysics |134j . We will address some current ideas later. 

Since a self-consistent model of core-collapse supernovae is still lacking, current 
stellar models induce the shock wave artificially by depositing a significant amount of 
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energy somewhere near the iron core. The magnitude of energy deposition is highly 
constrained, since the total predicted kinetic energy of the explosion should agree with 
observation. For example, observations of SN 1987A (a type II supernova with an 
estimated mass near 20Mq) and other type II supernova light curves result in explosion 
energies of ~ (1 — 2) x 10^^ erg. The location of artificial energy deposition is also 
constrained by observation: it can neither be located inside the iron core or otherwise 
neutron-rich iron group nuclides are overproduced, nor can it be located beyond the 
oxygen burning shell or the resulting neutron star mass, after fall-back of matter, will 
be too large. In most simulations, the mass cut, that is, the boundary between ejected 
and fall-back matter, is located in the silicon layer (inner dashed line in top part of 
Fig. |6]). The shock wave moves through the star and heats matter to high temperatures 
for a time period of a few seconds. Subsequently, the hot and dense matter expands 
nearly adiabatically. As a result, the star experiences several episodes of explosive 
nucleosynthesis. The character of the explosive nuclear burning depends, among other 
things, on the location of the shock and the expansion time scale. During the explosion 
the nuclides that have been synthesized before and after the core collapse are ejected 
and are then mixed into the interstellar medium. In fact, core-collapse supernovae are 
predicted to be the source of the majority of nuclides {A > 12) in the Galaxy. The 
situation is shown schematically in Fig. |6] (top right) and will be discussed in more 
detail below. 

6.3. Explosive nucleosynthesis 

The first zone encountered by the shock, consisting mainly of ^^Si (Fig. |6]), is heated 
to temperatures of T > 5 GK, while the density exceeds p ~ 10*^ g cm~^. This 
matter undergoes complete explosive silicon burning, during which the rates for all 
forward and reverse strong or electromagnetic interactions achieve a global equilibrium. 
This condition is called nuclear statistical equilibrium (NSE) and has some interesting 
properties |135j . The abundance of any nuclide in NSE is independent of initial 
composition or thermonuclear reaction rates, but is instead given by its nuclear 
properties, such as binding energy, mass and spins. In fact, these properties are well- 
known for the nuclides participating in explosive silicon burning [27] . For known nuclear 
properties, the abundance of any nuclide in NSE is determined by three independent 
parameters: temperature, density, and neutron excess. The latter quantity is a measure 
for the number of excess neutrons per nucleon and can only change as a result of 
weak interactions. It is defined hj rj = ^j(A^i — Zi)Xi/Mi, where Ni, Zi, X, and M, 
denote the number of neutrons, number of protons, mass fraction and relative atomic 
mass (in atomic mass units), respectively, of nuclide i in the plasm4§| In a 25Mq pre- 
supernova star of initial solar composition the explosively burning silicon layer has a 

§ The neutron excess parameter is related to the electron mole fraction, Yg, by 77 = 1 — 2Ye. For 
example, a plasma consisting only of nuclei with equal numbers of protons and neutrons has 77 = or 

Fe = 0.5. 
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very small neutron excess of i] ~0.003. Under such conditions, NSE favors ^^Ni as the 
main constituent |136] . since it has the largest binding energy of any N = Z species 
among the iron peak nuclides. Other iron group species are synthesized as well, such 
as ^'^Ni, ^^Co and ^"^Fe. Note that exactly the same products would be obtained if the 
explosively burning nuclear fuel in this temperature regime would not consist mainly 
of ^^Si, but instead of any other species with t] ~0.003. The fate of matter in NSE 
depends on the subsequent expansion as the gas quickly cools and, consequently, nuclear 
reactions will start to fall out of equilibrium at a certain freeze-out temperature. If 
the density and expansion time scale are sufficiently large, NSE predicts only a very 
small light particle (p, n, a) abundances. As the temperature falls below the freeze- 
out temperature, there is not enough time to produce the light particles necessary to 
maintain NSE and the equilibrium will be terminated by a lack of light particles. The 
first reaction to drop out of equilibrium is the triple-a reaction. There are so few light 
particles that their subsequent capture during this normal freeze-out does not alter the 
NSE composition (mainly ^^Ni and other iron peak nuclides). On the other hand, if the 
density or expansion time scale are sufficiently small, NSE predicts a relatively large 
light particle abundance, especially for a-particles. When the temperature falls below 
the freeze-out temperature, the a-particles cannot be converted fast enough to iron 
peak nuclei in order to maintain NSE. The equilibrium is terminated by an excess of 
a-particles and their capture by nuclei during freeze-out will alter the NSE composition 
|137j . For a nuclear fuel with rj ~0, the main product of this a-rich freeze-out is again 
^^Ni, but additional species are produced, for example, the important 7-ray emitter "^^Ti. 

The next zone encountered by the shock, still composed mainly of ^^Si, is heated 
to temperatures in the range of T ~4-5 GK. Instead of global equilibrium among all 
species, equilibrium is established locally in two mass ranges: the region centered on 
silicon, and the iron peak. This local equilibrium in each of these regions is called quasi- 
equilibrium (QSE). The abundance of any nuclide in QSE, for example, with ^^Si, is 
specified by four parameters (apart from nuclear binding energies, masses and spins): 
temperature, density, neutron excess, and the abundance of ^^Si |138] . When the shock 
first encounters this zone, it starts to dissociate silicon and neighboring species, quickly 
forming a QSE cluster in equilibrium with the dominant species ^^Si, which has the 
largest binding energy of all nuclides in this region. Soon thereafter another QSE cluster 
forms in the iron peak region, which consists of nuclides of even higher binding energies, 
mediated by reactions in the A = 40— 44 mass range linking the two clusters. If sufficient 
temperatures and time would be available, all of the silicon group species would be 
transferred to the iron peak and NSE would be established. However, the expansion 
causes the freeze-out to occur before NSE can be established. Since a significant amount 
of ^^Si remains, the process is referred to as incomplete silicon burning |139j . The major 
products of this burning are iron peak species, ^^Si, and intermediate-mass nuclides. 

Subsequently, the shock reaches a zone composed mainly of ^^O (Fig. [6]). Explosive 
oxygen burning occurs mainly at temperatures in the range of T ~3-4 GK. The character 
of the burning is similar to incomplete silicon burning, in the sense that the fuel (^^O) 
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is dissociated, giving rise to two QSE clusters in the mass regions of silicon and the iron 
peak. However, since the peak temperature is lower, less matter is converted to the iron 
peak and much more material remains locked in the silicon region. After freeze-out the 
most abundant nuclides in this zone are ^^Si, ^^S, ^^Ar and ^°Ca. Finally, the shock 
encounters a zone mainly composed of ^^O, ^°Ne and ^^C, and heats it to temperatures 
of T ~2-3 GK. As a result, ^''Ne, and to a lesser extend ^^C, burns explosively. In 
this case, the temperature and expansion time scale are too small for establishing QSE 
and the forward and reverse nuclear reactions operate far from equilibrium. Hence 
the abundance of a given species does not depend on a few parameters only, such as 
temperature, density, and neutron excess, but in addition is sensitively influenced by the 
initial composition and the magnitude of the thermonuclear reaction rates. After freeze- 
out the most abundant species in this zone are ^^O, ^°Ne, ^^Mg and ^*Si. Stellar models 
|131] also predict that explosive Ne/C burning is the main source for synthesizing the 7- 
ray emitter ^^Al, which has been observed in the Galactic interstellar medium |140[ 1141] . 
This observation was of paramount importance since it provided unambiguous direct 
evidence for the theory of nucleosynthesis in stars. The half-life of ^^Al amounts to 
7.17 X 10^ yr [27], which is small compared to the time scale of Galactic chemical 
evolution (~ 10^° y), and thus nucleosynthesis must be occurring currently throughout 
the Galaxy. The outer zones of the star are heated to peak temperatures of less than 
T ~2 GK for very short periods of time and, consequently, do not experience significant 
nucleosynthesis. These layers, located beyond the outer dashed line in Fig. [6], are ejected 
with their composition resulting from various hydrostatic burning stages prior to the 
explosion (left side of figure). About one hour after core collapse, the shock, traveling at 
an average speed of several thousand kilometers per second, reaches the stellar surface. 
The sequence of events sketched above depends sensitively on all factors that 
influence the mass-radius relation (or density profile) of the pre-supernova star, since it 
determines the amount of matter exposed to the four explosive burning episodes. For 
example, the treatment of convection is important in this regard, because it impacts the 
size of the convective regions and the mixing efficiency in those zones. Also, as already 
mentioned in Sec. O the ^^C(a,7)^^0 reaction rate determines the ^^C/^^0 abundance 
ratio at the end of core helium burning and hence impacts the amount of fuel available 
for the subsequent advanced hydrostatic shell burnings and the location of those shells. 
Furthermore, the magnitude of the time delay between core collapse and the revival of 
the shock impacts the amount of matter located in each compositional layer. Besides 
the mass-radius relation, the neutron excess (or Yf.) profile of the pre-supernova star is 
crucial for the outcome of those explosive burning episodes that take place under the 
conditions of NSE and QSE. The neutron excess profile, in turn, is infiuenced by the 
treatment of convection, the time delay between core collapse and shock revival, and 
the previous hydrostatic evolution of the star. 
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6.4- Observations 

Direct evidence for the nucleosynthesis predicted by current core-collapse supernova 
models can be obtained from observations of radioactivity and neutrinos. For example, 
the long-lived radioisotopes ^^Al (Ti/a = 7.17 x 10^ yr [27]) and ^^Fe (T1/2 = 2.62 x 10^ 
yr [142] ) are both predicted to be synthesized mainly by core-collapse supernovae. 
Their half-lives are very long compared to a typical frequency of about two Galactic 
core-collapse supernovae per century and, therefore, these nuclides accumulate after 
ejection from thousands of supernovae in the Galactic interstellar medium. The diffuse 
7-ray emission from the decay of both species has been observed by detectors onboard 
satellites. The ^°Fe/^^Al 7-ray line flux ratios measured by both the RHESSI instrument 
|143] and the SPI spectrometer onboard INTEGRAL |144j are near ~0.15. Stellar model 
predictions |145l 1131] seem to be consistent with observation, when the theoretical yields 
from single massive stars are folded with a standard distribution of stellar masses (initial 
mass function). However, it must be emphasized that the theoretical yields are sensitive, 
among other things, to current uncertainties in nuclear reaction rates |146l 1147] . 

The radioisotope ^^Ti has a half-life (T1/2 = 60 yr |148] ) comparable to the Galactic 
supernova rate and thus probes individual supernovae. Gamma-rays from the decay of 
^^Ti have been detected ( [149J , and references therein) in the 330-year old supernova 
remnant Cassiopeia A, which may have been observed by the English astronomer 
Flamsteed in the year 1680, confirming the theoretical prediction that this nuclide is 
synthesized in the a-rich freeze-out during explosive (complete) silicon burning in core- 
collapse supernovae. However, it is not clear at present why '^^Ti has not been observed 
in any of the other supernova remnants. It has been speculated that perhaps Cassiopeia 
A resulted from an asymmetric explosion [150] or represents a peculiar event [151] . 

The best-studied supernova to date, SN 1987A, resulted from the core collapse of a 
progenitor star with a mass near 2OM0. A few months after the explosion, the evolution 
of the brightness (light curve) followed the radioactive decay of ^^Co (T1/2 = 77.2 d) 
[152] . the daughter nuclide of radioactive ^^Ni (T1/2 = 6.1 d) that is the main product 
of explosive silicon burning at low neutron excess. Subsequently, 7-rays from ^^Co 
decay were directly detected with the SMM satellite [153] . Together, both observations 
implied a total ^^Co mass of ~0.07Mq in the ejecta [154] . Finally, a few hours before 
the light from SN 1987A reached Earth, 11 electron antineutrinos were recorded by 
the KamiokaNDE-H detector [155] . 8 electron antineutrinos by the Irvine-Michigan- 
Brookhaven (1MB) detector [156j . and 5 electron antineutrinos by the Baksan Neutrino 
Observatory [157] . The number of detected neutrinos, their energies, and the measured 
burst duration are in agreement with theoretical predictions of the processes deep inside 
an exploding star (see below). 

6.5. p-process 

There are about 30 relatively neutron deficient nuclides with mass numbers of A > 74 
(between '^^Se and ^^^Hg) that cannot be synthesized by any of the known neutron 
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capture processes. These species are called p-nuclei, and they all contain more protons 
relative to other stable isotopes of the same element. The mechanism responsible for 
the synthesis of the p-nuclei is called the p-process. It turns out that the p-nuclei are the 
rarest among the stable nuclides. In fact, no single element has a p-process isotope as 
a dominant component. In principle, there are two kinds of reactions that could allow 
for the production of neutron-deficient nuclei starting from heavy seed nuclei: proton 
captures, (p,7), and photo-neutron emission, (7,n). The first mechanism requires high 
densities, temperatures, and relatively long time scales, which are unlikely to exist in 
any hydrogen-rich zones of common stars |158j . Thus it is generally accepted that the 
main mechanism of p-nuclei production involves the photodisintegration of heavy seed 
nuclei in a hot photon environment (peak temperatures near T ~2-3 GK), a hydrogen- 
exhausted stellar zone, and a short time scale. The abundance flow reaches from 
lead down to the iron peak, where further photodisintegrations become energetically 
unfavorable, all subject to the same hot photon environment. During this process, the 
seed nuclei are initially destroyed by several (7,n) reactions along a given isotopic chain. 
This continues until a neutron deficient nucleus is reached that prefers a (7,p) or (7,0;) 
over a (7,n) reaction, thus branching the abundance flow into an isotopic chain of a 
lighter element. Abundances tend to accumulate at these branching points, especially 
at nuclides with a closed neutron or proton shell structure (that is, at species with small 
photodisintegration rates). These waiting points become either directly p-nuclei (in 
the region of the lighter p-nuclei) or transmute via /3+-decays to p-nuclei after cooling, 
expansion and ejection of material (in the region of the heavier p-nuclei). It is important 
for any realistic site responsible for the synthesis of the p-nuclei that the combination of 
temperature and time scale guarantees the occurrence of some photodisintegrations, yet 
not so intense as to entirely reduce all nuclei to iron peak species. These arguments 
support the conclusion that the p-process occurs during stellar explosions with an 
associated rapid expansion and cooling of material. 

Although the site(s) of the p-process have not been identified unambiguously yet, 
many investigations have assumed that it occurs in core-collapse supernovae, when the 
shock wave passes through the oxygen-neon zone (Fig. [6]) and heats the matter to 
temperatures between ~2-3 GK [159j. During the explosion, different layers of this 
zone will undergo different thermodynamic histories and thus will achieve different 
peak temperatures. The weak s-process component operating during the preceding 
core helium and shell carbon burning stages (Secj5]) strongly enhances the p-process 
seed abundances in the A ~60-90 mass region. It has been demonstrated that p-nuclei 
with masses of A < 92, A ~92-144 and A > 144 are mainly produced in stellar zones 
with peak temperatures of T > 3 GK, T ~2. 7-3.0 GK and T < 2.5 GK, respectively. 
In fact, each p-nuclide is synthesized in a relatively narrow temperature range only 
|160] . The simulations reproduce the solar system abundance of most p-nuclei within a 
factor of 2-3. However, there is a notable underproduction of the light p-nuclei ^^Mo, 
^^Mo, ^^Ru and ^^Ru, as well as of the rare species ^^^In, ^^^In and ^^^La. Sites other 
than massive stars have also been considered and, interestingly, similar p-abundance 
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distributions are obtained in each case. Some of the underproduced p-nuchdes may be 
predominantly synthesized by a different mechanism (see below). For more information 
on the p-process, including a description of other proposed sites, see Ref. |161] . 

A number of current nuclear physics uncertainties impact p-process simulations. For 
example, the ^^C(a,7)^^0 rate influences the pre-supernova evolution of the massive 
star and hence the composition of the oxygen-neon zone prior to core collapse |162 



The ^^Ne(a,n)^^Mg reaction (Sec [5]) is important since it is responsible for the weak 
s-process component, which provides a fraction of the seed-nuclei for the p-process. 
With few exceptions, almost all of the many reactions (> 10, 000) in the region of the 
p-process, including the important (7,n), (7,p) and (7, a) photodisintegrations, have to 
be computed theoretically using the Hauser-Feshbach model of nuclear reactions. It 
has been shown that different prescriptions of the Hauser-Feshbach model (ct-nucleus 
and nucleon-nucleon optical potentials, and level densities) influence sensitively the 
final predicted p- nuclei abundances. Experimental (n,7), (p,7) and (0,7) rates on 
stable target nuclei in the mass range of A > 60 play an important role in this 
regard. The experimental cross sections are used for adjusting Hauser-Feshbach model 
parameters, thus improving theoretical rate predictions for a multitude of unmeasured 
reactions. Furthermore, the reverse photodisintegration rates can be calculated from 
the measured forward rates using the reciprocity theorem. For a sensitivity study of 
p-process nucleosynthesis to nuclear reaction rates, see Ref. jl63j . 

6.6. u-process 

The large energy release in the core collapse in form of neutrinos (see below) gives 
likely rise to the neutrino-induced synthesis of certain nuclides in the expanding mantle 
of the exploding star |164j . Specifically, neutrinos can interact with the nuclei via 
inelastic neutral-current neutrino scattering, (z/,z/'). Since /i and r neutrinos emerging 
from the proto-neutron star have larger predicted average energies compared to electron 
neutrinos (see below), and since the neutrino cross sections scale with the square of 
their energy, the interactions of the former neutrino species with nuclei predominate. 
On the other hand, electron neutrinos may interact with nuclei via charged-current 
interactions, (z/e,e~) or {i?e,e~^). All these neutrino interactions can populate excited 
nuclear levels that subsequently decay via emission of light particles (p, n, a, etc.). 
The released light particles, in turn, undergo reactions with other nuclei in the high- 
temperature environment and thus contribute to the synthesis of certain nuclides. In 
any given layer of the star, neutrino induced nucleosynthesis may occur before the 
arrival of the shock, during explosive burning, or after shock passage when the material 
expands and cools. This mechanism is referred to as the i^-process [165] • Simulations 
have shown [1651 11661 1167] that it may contribute appreciably to the solar abundance of 
the rare species ^^B (made in the carbon-oxygen layer; see Fig. [6]), ^^F and ^^^La (both 
made in the oxygen-neon layer). Recall that the latter nuclide is underproduced in the 
p-process. Yield predictions from //-process calculations are sensitively influenced by 
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current uncertainties in neutrino interaction cross sections, average neutrino energies of 
each flavor, total neutrino luminosity, as well as by the details of the adopted stellar 
evolution and explosion models. 

6. 7. Core collapse revisited. Role of neutrinos 

As discussed above, core-collapse supernovae are powered by the release of gravitational 
energy in the collapse of a massive star core to a proto-neutron star. We will now discuss 
the core collapse in more detail. At the onset of the collapse (t=0), temperatures 
and densities in the center of the iron core amount to T ~10 GK and p ~ 10^° 
g cm~^, respectively. Under these conditions, nuclear statistical equilibrium (NSE) 
is established. The silicon burning shell has continually increased the mass of the 
core, which is supported by electron degeneracy pressure. When the core exceeds the 
Chandrasekhar mass (~1.4Mq), it has no other thermonuclear energy source available 
to support the pressure and it becomes unstable to gravitational collapse. The collapse 
and the infall dynamics depend sensitively on two parameter^jj] : (i) the electron mole 
fraction. Ye, which is related to the neutron excess (see above); and (ii) the entropy per 
baryon, s. During the early stages, the core collapse is accelerated by two effects. First, 
as the (electron) density increases, electrons capture onto nuclei, (e~, Ue), hence removing 
electrons that were contributing to the pressure. As a consequence of a decreasing Ye, the 
neutron excess, 77 = 1 — 2Ye, increases in the core and, in addition, a burst of electron 
neutrinos is produced. Second, at the very high temperatures the thermal radiation 
becomes sufficiently energetic and intense that iron peak nuclei are photodisintegrated 
into lighter and less stable species, thus removing energy that could have provided 
pressure. Within a fraction of a second, the core with a size of several thousand 
kilometers collapses to a proto-neutron star of several tens of kilometer radius. The 
most important neutrino interactions during the collapse are (neutral current) elastic 
scattering on nuclei, {uey^e), electron-neutrino scattering, e~{h'e,^e)G~ , inverse /9-decay, 
{h'e,e~), and inelastic scattering on nuclei, (t'e, t'e) |170] . At t ~0.1 s, when the density 
reaches a value near ~ 10^^ g cm~'^, the neutrino diffusion time becomes larger than the 
collapse time and, consequently, the neutrinos become trapped |169j . Inside this region, 
called neutrino sphere, the neutrinos are coupled to the matter and are in thermal 
equihbrium. Outside, the neutrinos escape almost freely, with an average energy that 
is determined at the radius of the neutrino sphere, Ru. (Note that the location of the 
neutrino sphere depends both on neutrino type and flavor, and on neutrino energy.) 

II The electron mole fraction is equal to the electron-to-baryon ratio, Ye — N^/^^NiMi, and hence 
to the proton-to-baryon ratio, Yp] Ne, Ni and Mi denote the electron number density, nuclide number 
density and nuclide relative atomic mass, respectively, where the sum i is over all nuclides. For a 
radiation dominated environment, the entropy per baryon, s, is related to the photon-to-baryon ratio, 
(j), via s ~ 10(/) ~ T^ I p, where T and p denote the temperature and density, respectively. In simple 
terms, for a small value of the entropy per baryon, NSE favors a composition of iron peak nuclei, which 
are tightly bound. On the other hand, a large value of s implies that many photons are available per 
baryon, which favors the photodisintegration of heavier nuclei into free nucleons [168j . 
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At t ~0.11 s, the inner core (M ~O.5M0) reaches nuclear densities (~ lO^'^ g cm"^), 
bounces, and drives a shock wave into the infaUing matter. At the time of the bounce, 
the electron mole fraction in the inner core is relatively small {Ye ~0.3 for M < O.SM©), 
and increases gradually in the outer core {Yf. ~0.45 near M ^I.OMq). At t ~0.12 s, the 
prompt shock propagates outward, but loses severely energy (~9 MeV per nucleon) by 
dissociating iron peak nuclei into free nucleons. When the shock reaches the neutrino 
sphere, additional electron captures on free protons also remove energy from the shock, 
giving rise to a strong burst of electron neutrinos (prompt u^ burst, corresponding to 
~ 10^3 erg s"^ for ~10-20 ms). At t ~0.2 s, the shock stalls at a radius of ~100-200 
km in the outer core. 

During core collapse, a gravitational binding energy of several 10^'^ erg, representing 
a staggering ~10% of the iron core's rest mass, is released in form of neutrino radiation. 
Therefore, the stalled shock is thought to be revived by neutrinos and antineutrinos 
that emerge from the hot and dense proto- neutron star |171l 1172] . At the very 
high temperatures prevailing in the core, neutrinos and antineutrinos of all flavors 
(electron, muon and tau) are produced, mainly by electron-positron pair annihilation, 
e~ + e"*" — )■ Ue + i^e, neutrino- antineutrino pair annihilation, Ue + i^e ^ i^ti,T + ^'^.t, and 
nucleon bremsstrahlung, N + N — )■ A^ + A^ + z/^ ,- + 1/^ ,- [173] . While the neutrinos diffuse 
out of the core, the luminosities and average energies for each of the different neutrino 
types evolve with time. First, muon and tau neutrinos have smaller opacities compared 
to electron neutrinos since their energies are too low for charged-current interactions, 
such as u^+p ^ n + iJ,'^ (because of the large masses of the /x and r leptons). Thus, the 
muon and tau neutrinos decouple at smaller radii Ri, (at higher density and temperature) 
and hence emerge with higher average energies from the proto-neutron star compared 
to the electron neutrinos (see i/-process). Second, there are fewer protons than neutrons 
in the outer layers of the proto-neutron star and thus the charged-current interaction 
n + Ue ^ p+e~ occurs more frequently compared to p + i^e ^ n + e~^ . Hence, the electron 
neutrinos have a higher opacity than the electron antineutrinos and thus decouple at 
larger radii with smaller average energies. 

The situation at this stage is represented in Fig. [7] (left). The matter between 
the neutrino sphere and the stalling shock front consists mainly of free neutrons and 
protons. The gain radius divides this region into two parts: the first is located closer 
to the neutrino sphere and is characterized by the dominance oi p + e~ — )■ n + u^ 
and n + e~^ — )■ p + z/g over their reverse interactions, thus giving rise to effective 
neutrino cooling by neutrino emission; the second is located closer to the shock, where 
n + Ue — 7- p + e~ and p + Ue ^ n + e~^ dominate over their reverse interactions, hence 
giving rise to effective heating via neutrino absorption. The continuous neutrino energy 
deposition in the latter region keeps the pressure high and could then rejuvenate the 
shock and cause the supernova explosion (delayed shock model [1741 I175J ). Only a 
fraction (~1%) of the total gravitational binding energy, deposited by neutrinos as 
thermal energy of nucleons, leptons and photons in this region, would be required to 
initiate a powerful shock propagating through the stellar mantle and giving rise to an 
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Figure 7. Evolution of the region near the proto-neutron star (PNS) surface, as 
a result of neutrino-induced processes, expansion and cooling. The proto-neutron 
star surface can be defined by the radius of the energy-integrated electron neutrino 
sphere. The radii of the various zones differ from panel to panel and are not to scale. 
Approximate temperature ranges (in GK) are indicated. Only the main constituents 
are given in each zone. Symbols in green boxes indicate nucleosynthesis processes: 
i^p-process (t'p); a-process (a); r-process (r). (Left) Region between proto-neutron 
star surface and stalled shock, t '^0.2 s after core collapse. The main constituents 
are free protons and neutrons. Inside the gain radius neutrino cooling prevails, while 
outside matter is effectively heated by neutrinos. Radii: Ry (neutrino sphere); Rg 
(gain radius); Rs (shock radius). (Middle) The early neutrino-driven wind, t '^l.O, is 
proton-rich, Yg > 0.5, giving rise to the fp-process. (Right) The later neutrino-driven 
wind, t ^10, may become neutron-rich, Ye < 0.5, giving rise to the a- and r-processes. 
Note that the actual abundances of ^Be and ^^C are very small since these species are 
destroyed almost as quickly as they are produced. The label "Fe" denotes a distribution 
of seed nuclei in the A ~50-100 mass range. See text. 



explosion. The success of this model depends crucially, among other things, on the 
product of neutrino luminosity and neutrino interaction cross sections (i.e., the square 
of the average neutrino energy). However, almost none of the most advanced, self- 
consistent core-collapse models have produced an explosion and the exact mechanisms 
by which neutrinos give rise to the explosion remain elusive. Instead, many models 
change artificially the neutrino properties, such as the charged-current interaction rates, 
in order to increase the neutrino energy deposition behind the shock. The problem is 
highly complex, involving energy-dependent neutrino transport in three dimensions, a 
convectively unstable region near a compact hot and dense object, possible diffusive 
instabilities, magneto-rotational effects, and so on. It appears that a self-consistent 
3D simulation would require a computing power at least ten times as large as what is 
currently available. 

We will consider below nucleosynthesis processes that effect the deepest layers to be 
ejected in a core-collapse supernova. After the rejuvenated shock has been launched, the 
strong electron neutrino and antineutrino fluxes drive a continuous flow of protons and 
neutrons from the region near the proto-neutron star surface. This outflow of matter is 
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referred to as neutrino- driven wind |176] . Matter in this region expands and cools at high 
(possibly supersonic) velocity and eventually collides with the slower earlier supernova 
ejecta, resulting in a wind termination (or reverse) shock |177] . The nucleosynthesis in 
the early ejecta depends strongly on the entropy (photon-to-baryon ratio), expansion 
timescale, and electron mole fraction (electron-to-baryon ratio). The first two quantities 
depend on neutrino average energies and total luminosities, and the properties of the 
proto-neutron star (given by the equation of state). The latter quantity depends on the 
complicated interplay of the four (charged- current) electron neutrino and antineutrino 
interactions with free nucleons, n + u^ ^ p-\-e~ and p+ z/g ^ ^ + e"'", in the wind region. 
If the neutrino-driven wind is proton-rich, then conditions are favorable for the synthesis 
of neutron-deficient nuclides, while a neutron-rich wind may give rise to the production 
of neutron-rich species. The outcome of the nucleosynthesis is highly sensitive to the 
still uncertain physical conditions of the neutrino-driven wind. 

6.8. up-process 

Recent core-collapse simulations including energy-dependent neutrino transport have 
shown that the neutrino-driven wind is proton-rich {Y^ > 0.5), either at very early 
times |178j or for an extended time period up to 20 s |186j . The is shown schematically 
in Fig. [7] (middle). The wind is ejected at temperatures above 10 GK, consisting of 
free neutrons and protons in nuclear statistical equilibrium (NSE). Four distinct phases 
can be identified during the nucleosynthesis: (i) expansion and cooling through the 
temperatures range of T ~10-5 GK causes all neutrons to combine with protons, leaving 
a composition consisting of a-particles and an excess of protons; (ii) further cooling 
through temperatures of T ~5-3 GK allows the a-particles to combine to heavier nuclei, 
giving mainly rise to the production of the nuclides ^^Ni, ^°Zn and ^'^Ge, which consist 
of an equal number of neutrons and protons (A^ = Z). Since the slowest link in this 
sequence is the triple-a reaction, the number of heavy seed nuclei synthesized is equal 
to the number of ^^C nuclei produced. Abundance flows beyond ^^Ge are inhibited 
by strong reverse flows and by its long half-life (T1/2 = 64 s in the laboratory [27]): 
(iii) in the temperature range of T = 3 — 1.5 GK the (charged-current) interaction 
p + z/g — > ^ + e+ on the abundant protons produces free neutrons (typically 10^^ per cm^ 
for several seconds), which then participate in the nucleosynthesis. This aspect is very 
important, because fast (n,p) reactions on waiting points, such as ^^Ni, ^°Zn and ^^Ge, 
and subsequent (p,7) reactions allow for a continuation of the abundance flow toward 
heavier nuclides on the neutron-deficient side of the valley of stability. Note that a (n,p) 
reaction connects the same pair of nuclei as a corresponding /3+-decay; (iv) when the 
temperature falls below T ~1.5 GK, the (p,7) reactions freeze out and (n,p) reactions 
and /S^-decays convert the heavy nuclei to stable, neutron-deficient, daughters. This 
nucleosynthesis mechanism is called the i/p-process |187[ I188"| 1189] . 

The termination of the z/p-process is caused by the cooling of matter to temperatures 
below T ~1 GK, and not by the exhaustion of free protons. A few studies (for example. 
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Ref. |19U] ) have shown that the i/p-process may account for the solar abundance of 
the hght p-nuchdes up to ^°'^Cd, including ^^Mo, ^"^Mo, ^^Ru and ^^Ru, which are 
significantly underproduced in the p-process or other proposed scenarios (Sec. I6.5p . 
However, other studies |189[ 1191] could not reproduce the solar abundance of the hght 
p-nuclides and thus the issue remains controversial at present. It is also important 
to stress that, although many model simulations predict an early proton-rich neutrino- 
driven wind, the i^p-process is highly sensitive to the details of the explosion mechanism, 
the mass, and possibly the rotation rate of the proto-neutron star. Thus it is expected 
that the nucleosynthesis may vary considerably from event to event. Currently, the 
key uncertainties originate from the neutrino luminosities and average energies, which 
determine the electron mole fraction, Y^., in the neutrino-driven wind. A larger value 
of Ye will generally give rise to the synthesis of heavier p-nuclei. For example, models 
with Ye ~0.6 (i.e., at the onset of the z/p-process, when T ~3 GK) produce p-nuclei 
up to A = 152. For even larger, perhaps excessive, values of Ye, the density of free 
neutrons would increase as a result of the charged-current interaction between electron 
antineutrinos and the large proton abundance. Consequently, the (n,p) reactions would 
not only carry the abundance flow to the group of stable nuclides, but (n,7) reactions 
will, in fact, carry the flow to the neutron-rich side. The protons do not participate 
directly in the nucleosynthesis, but mainly act as a source of free neutrons |189j . 

The i^p-process is also sensitive to nuclear physics uncertainties. For example, the 
slow triple-a-reaction (Sec. 15. 2p represents a bottleneck for the synthesis of the heavy 
seed nuclei (before the start of the i/p-process, near T = 3 GK). Varying this rate will 
impact the production of the p-nuclei in the A = 100 — 110 mass region. Changing the 
rates of the (n,p) reactions on the waiting point nuclei that control the abundance flow to 
heavier nuclides, especially for ^^Ni(n,p)^^Co and ^°Zn(n,p)^*^Cu, strongly influences the 
production of the p-nuclei in the A > 100 mass region. In the most recent simulations 
[190] these (n,p) rates are based on theory (Hauser-Feshbach model), since experimental 
data are lacking at present. A study of the (modest) impact of nuclear mass uncertainties 
on the path of the z/p-process can be found in Ref. [192] . 

6.9. r-process 

About half of the nuclides beyond the iron peak, including gold, silver, platinum and 
uranium, are believed to be synthesized via neutron-induced reactions, under extreme 
conditions that are very different from those giving rise to the s-process (Sec. 15.31) . The 
abundances observed in the solar system and in (old) stars of very low metallicity hint at 
a mechanism that started to occur early in the history of the Galaxy, providing very high 
neutron densities at large temperatures for a very short period of time. This mechanism 
is called the astrophysical r(apid neutron capture)-process and its phenomenology has 
been understood for a long time (see, for example, Ref. |193j ). One of the most attractive 
sites, which has been studied intensively over the past two decades, in fact, is the 
neutrino-driven wind. We will first describe the phenomenological model and then 
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return to the connection between the r-process and the neutrino-driven wind. 

If one subtracts the main s-process component, estimated using modern models of 
AGB stars (Sec. [5]), from the measured solar system abundances in the A > 100 mass 
region, then the residual distribution exhibits two prominent peaks, at mass numbers 
near A = 130 and 195, about ten mass units removed from the s-process peaks near 
A = 138 and 208. The existence of the r-process peaks and of the long-lived radioisotopes 
232Th (Ti/2 = 1.4 X 10^° y), 235U (Ti/2 = 7.0 X 10« y), and ^^^V (T1/2 = 4.5 x 10^ y), 
which are located beyond the endpoint of the s-process, provide strong evidence for the 
existence of the r-process. The abundance peaks can be explained by the neutron magic 
numbers at N=82 and 126. Unlike the s-process, one may imagine a very strong neutron 
flux, such that neutron capture is more likely than /3~-decay. Thus the abundance flow is 
driven away from the stability valley toward the neutron-rich side. Here, the abundance 
accumulates at nuclides with neutron magic numbers [39]. After termination of the 
neutron flux in the r-process, these neutron-magic nuclei undergo series of /3~-decays 
along isobaric chains {A = const) until the most neutron-rich stable (or very long-lived) 
isobar is reached. Hence the r-process gives rise to abundance maxima in mass regions 
located below the corresponding s-process abundance peaks. 

Suppose that seed nuclei (for example, iron) are exposed to high neutron densities 
(A^„ > 10^° cm"^) and temperatures (T > 1 GK). Under such conditions, both 
(n,7) and (7,n) reactions are much faster than competing /3~-decays. Consequently, 
thermal equilibrium is quickly established along each chain of isotopes, implying that the 
abundance of a given isotope depends mainly on (n,7) reaction Q-values, temperature 
and neutron density, but is independent of neutron capture rates. Furthermore, the 
systematic behavior of the Q-values when moving from the stability valley toward the 
neutron drip line gives rise to an abundance maximum at one or at most two species in 
each isotopic chain. These nuclides represent waiting points for the abundance flows 
and, therefore, the (n,7)-f-)'(7,n) equilibrium condition is referred to as the waiting 
point approximation. The abundance flow must then continue with a slow /3~-decay, 
connecting one isotopic chain to the next, where again an independent equilibrium 
within the new isotopic chain is established. This repetitive sequence of events gives 
rise to the r-process path. The situation is analogous to a radioactive decay chain: in 
both cases the abundance flow will attempt to achieve a constant /3~-decay flow from one 
isotopic chain to the next. This condition is referred to as steady flow approximation. 
In this model, for given values of neutron density, temperature and time duration, 
the r-process path and r-process abundances can be calculated precisely, if the nuclear 
properties are known: neutron capture Q-values (or neutron separation energies), /3- 
decay half-lives, branching ratios for /3-delayed neutron emission, normalized partition 
functions, and so on. 

The phenomenological model discussed above is referred to as the classical r-process 
model. It is rather simple because it assumes: (i) a constant temperature, T, and neutron 
density, Nn, (ii) an instantaneous termination of the neutron flux after time r, and 
(iii) the waiting point and steady flow approximations. Nuclear masses play a central 
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role for the r-process since they determine directly or indirectly most of the nuclear 
properties listed above. The (n,7) Q-values (which are given by mass differences) enter 
exponentially in the determination of the equilibrium abundances of each isotopic chain 
and thus must be known rather accurately. It is obvious that the nuclear properties 
needed are for isotopes that are located far away from the valley of stability. With few 
exceptions, the required information is experimentally unknown since most of the nuclei 
on the r-process path cannot be produced yet in the laboratory. There is little alternative 
at present but to estimate the required nuclear properties by using nuclear models. In 
practice, attempts are made to derive semi-empirical formulas from the known properties 
of nuclei close to stability that can be extrapolated into the region covered by the r- 
process path. Such extrapolation procedures are subject to significant uncertainties 
even for the most sophisticated models (see, for example, Ref. |194j ). Obviously, any 
deficiencies in current nuclear models will have a direct impact on r-process predictions. 
The associated nuclear physics uncertainties affect most discussions of the r-process. 

Attempts to describe the entire observed distribution of solar system r-abundances 
by using the classical r-process model with a single set of T-N^-t conditions were 
unsuccessful |195j . It was found that at least three different sets of T-Nn-r conditions, 
each corresponding to a specific r-process path, are required in different mass regions 
in order to reproduce the observed solar-system abundance distribution. Temperatures, 
neutron densities and time scales for these components were in the ranges of T ~1.2- 
1.35 GK, A^„ ~ 3 X 10^° — 3 x 10^^ cm~^ and r ~1.5-2.5 s, respectively. Furthermore, 
it was demonstrated that the steady flow approximation applies locally in each of these 
mass ranges, but not globally over the entire mass region. The overall implication is that 
the solar system r-abundance distribution results from a superposition of components 
representing different r-process conditions. This may be caused, for example, by varying 
conditions in different zones of a given site. Clearly, the phenomenological model 
described above makes no assumption regarding the site of the r-process, but provides 
a general idea regarding some of its properties: very high neutron densities and short 
time scales. On the other hand, the temperature should not be too high or otherwise 
the heavy nuclei will be either destroyed by photodisintegration or the waiting point 
abundances shift too close to the stability valley, where the /3~-decay half-lives are too 
slow to allow for efficient r-processing. The classical r-process has also been applied 
for reproducing or predicting abundance ratios of neighboring nuclides, for example, for 
nuclear chronometers |196] or isotopic anomalies in primitive meteorites |197] . Such 
abundance ratios are most likely influenced by nuclear properties rather than by the 
details of the astrophysical r-process site. For an extensive discussion of the limitations 
of the classical r-process model, see Ref. |198j . 

It is worthwhile to visualize the results of a dynamic r-process calculation as opposed 
to the static model described above. Clearly, in reality the temperature and neutron 
density will change with time. Early during the r-process, the temperature and neutron 
density will be sufficiently large to ensure that (n,7)-H'(7,n) equilibria hold in all isotopic 
chains. Suppose now that the temperature and neutron density decrease with time. The 
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abundance flow will continuously adjust to the new conditions according to the waiting 
point approximation. Thus the r-process path, defined by T and Nn, must continuously 
move, starting from a location closer to the neutron drip line to one that is located closer 
to stability. For each location of the path, the f3~ -decay half-lives are different. When 
the neutron flux ceases, the r-process nuclei decay towards stability. It is apparent that 
only the r-process path just before freeze-out, and in particular the sections near neutron 
magic nuclei, matters for the observed flnal distribution of r-abundances. In other words, 
at freeze-out the r-process has mostly forgotten its earlier history. Furthermore, if the 
condition of an instantaneous termination of the neutron flux is relaxed, additional 
nuclear processes that may occur during the freeze-out will likely become important. 
For example, neutron capture reactions and /3-delayed neutron decays (apart from /3~- 
decays) may alter the flnal composition. 

In addition to studying the solar abundance distribution, important clues regarding 
the r-process are provided by stellar spectroscopy of metal-poor (and thus very old) 
Galactic halo stars |199[ 1200] . For example, observed elemental abundances for the 
remarkable halo giant star CS 22892-052 are shown in Fig. |8l It was found in this 
case, and for several other low-metallicity halo giant stars, that for elements of barium 
{Z = 56, A ~135) and beyond the observed heavy element abundances (data points) 
agree remarkably well with relative solar system r-abundances (solid line). This implies 
that most of the heavy elements observed in those stars were synthesized by the r- 
process early during the evolution of the Galaxy, with no apparent contribution from 
the s-process. Since the r-process elements could not have been synthesized in the halo 
stars themselves, they must have been produced by progenitors that evolved very rapidly 
and that ejected the matter into the interstellar medium before the formation of the 
currently observed halo stars. In fact, these stars formed so early in the history of the 
Galaxy that they may have received contributions from only one, or at most a few, 
r-process events. The most likely r-process sites seem to be associated with massive 
stars since low-mass or intermediate-mass stars evolve on much longer time scales. 
The agreement in abundances above A ~135 shows that the r-process mechanism is 
robust, in the sense that a similar abundance distribution pattern is produced in each 
r-process event. Interestingly, the good agreement between solar r-process and stellar 
abundances does not extend to the lighter elements between strontium (Z=38) and 
barium (Z=56). There have been many suggestions to explain their synthesis, including 
a "weak r-process" [I98], a "Light Element Primary Process" (LEPP) [2011 [202], and 
the a-process |203j . Clearly, more observations and studies will be needed to understand 
the synthesis of these elements. 

We will now describe the nucleosynthesis that may occur in the late neutrino-driven 
wind, assuming that conditions are favorable for the r-process. Four different phases 
can be identifled, as shown schematically in Fig. [7] (right): (i) above a temperature 
of T ~10 GK, the wind consists of free neutrons and protons in nuclear statistical 
equilibrium (NSE); (ii) as the wind expands and cools below T ~10 GK, the nucleons 
start to combine to a-particles. Near T ~7 GK, a-particles become the dominant 
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Figure 8. (Top) Elemental abundances of the low-metallicity halo giant star CS 
22892-052 (solid circles) compared to the solar system r-process abundances (solid 
line). The abundances are given in standard spectroscopic notation, where loge{A) = 
logiQ{NA/NH) + 12.0, and Ni denotes number abundance. The stellar data are from 
Ref. |199| . The solar system r-abundances (total solar system abundance minus s- 
abundances) are from Ref. [204] and have been normalized to the stellar Eu abundance. 
Note that the graph does not include the uncertainties introduced by the deconvolution 
of s- and r-abundances. (Bottom) Difference between total stellar and scaled r-process 
solar system abundances. The large and small scatter in the regions Z < 50 and 
Z > 50, respectively, is apparent. See text. 



constituent, leaving behind an excess of neutrons. Furthermore, some of the a-particles 
begin to assemble into nuclei via the strongly density- dependent process a+a+n -^^Be, 
followed by ^Be{a,ny^C; (iii) near T ~5 GK, nuclear statistical equilibrium starts to 
break down because the expansion time scale becomes shorter than the time required 
to maintain nuclear statistical equilibrium under conditions of high temperature, low 
density and large a-particle abundance. The slow process a + a + n -^'^Be is the first 
reaction to fall out of NSE. The situation is similar to the a-rich freeze-out (Sec. I6.3p . 
There is again an excess of a-particles and these cannot be consumed fast enough 
by the slow helium-induced reactions in the time available. The important difference, 
however, is that there is now a neutron excess. Both the abundant a-particles and the 
neutrons will participate in the buildup of heavier nuclei, starting with the sequence 
^'^C{n,'-)Y^C{a,nY^O{a,'^Y^Ne, and so on. Without the presence of neutrons, the a- 
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rich freeze-out produces mainly the N = Z isotope ^^Ni, while the flows beyond the 
iron peak would be negligible (Sec. I6.3p . The capture of neutrons has the important 
effect of shifting the abundance flow towards the stability valley, where the nuclei are 
more neutron-rich, thus extending the flow beyond the iron peak up to A ~110. The 
most important processes in the buildup of these heavy nuclei are (a,n) and (n,7) 
reactions. This neutron-rich, a- rich freeze-out is referred to as the a-process; (iv) at a 
temperature near T ~3 GK, the a-particle induced reactions become too slow to change 
the composition of the matter and the a-process ceases. At this time, the composition 
resides in a-particles (mostly), neutrons and seed nuclei in the A ~50-110 mass range. 
If the neutron-to-seed ratio is sufficiently large, an r-process can be launched while the 
matter further expands and cools. 

The r-process model described above has the advantage that the properties of the 
neutrino wind are determined by the proto-neutron star, not by the pre-supernova 
evolution. On the one hand, this implies that r-processing in this site should produce 
similar abundances for events involving proto-neutron stars of the same mass. On the 
other hand, the r-process will be sensitive to the proto-neutron star evolution and the 
neutrino interactions, and thus represents an important diagnostic of the event. It 
is crucial for a successful r-process that the preceding a-process is not too efficient. 
Otherwise, too many heavy seed nuclei are produced and too many neutrons are 
consumed, resulting in insufficient neutron-to-seed ratios for synthesizing nuclides up 
to A ~200 during the subsequent r-process. This requirement translates into relatively 
high entropies (or low densities) in the neutrino-driven wind so that the a + a + n -^^Be 
reaction is less efficient in converting a-particles to heavier nuclei. A fast expansion 
time scale limits the duration over which the freeze-out operates and is also helpful for 
reaching a high neutron-to-seed ratio. To summarize, an appropriate combination of 
low electron mole fraction, high entropy (i.e., high temperatures and low densities), and 
short expansion time scales is a necessary condition for a successful r-process in the 
neutrino-driven wind model. 

Some studies have attempted to reproduce the solar system r-process abundance 
distribution by assuming a superposition of many individual contributions in the 
neutrino-driven wind, each representing an appropriate combination of electron mole 
fraction, entropy and expansion time (see, for example, Ref. |205] ). It was found 
that the component with the heavier r-abundances {A > 120), including Th and U, 
could be reproduced for electron mole fractions below Y^ ~0.48 for modest entropies 
{s <300). Furthermore, it could be confirmed that the waiting point and steady flow 
approximations of the classical r-process apply locally (but not globally) over certain 
mass regions in the dynamic simulations. Recent studies |206j (see also Ref. [207] ) 
have found that many (though not all) nuclides between zinc {Z = 30, A ~65) and 
ruthenium [Z = 44, A ~100) can may be synthesized in the neutrino-driven wind, under 
suitable low- entropy conditions, via the a-process when the abundance of free neutrons 
at charged-particle freeze-out is negligible (i.e., no subsequent r-process occurs). The 
results are interesting because this scenario co-produces via a charged-particle process 
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a fraction of the A ~65-100 nuclides, which are commonly thought to originate from 
different processes, that is, the s-, r- and p-processes. Clearly, more abundances from 
stellar spectroscopy and more extensive nucleosynthesis simulations are called for. 

Calculations of r-process nucleosynthesis require an extensive input of nuclear 
physics information, such as nuclear masses, /9~-decay half-lives, branching ratios for /3- 
delayed neutron decay, partition functions, and so on. If the freeze-out from equilibrium 
is followed explicitly (that is, if the waiting point and steady flow approximations are not 
applied and the reaction network is solved numerically) then reaction rates for neutron 
captures and photodisintegrations are required as well. All of this information is needed 
for neutron-rich isotopes that are located far away from stability. If the a-process needs 
to be followed explicitly, then another large data set consisting of rates for charged- 
particle reactions and neutron- induced processes, such as (n,a) and (n,p), is required. 
Needless to say that almost all of this information must be obtained from global, 
semi-empirical, models for nuclear masses and /9~-decays, and from Hauser-Feshbach 
calculations. Information on directly measured properties, mainly life times and /3- 
delayed neutron emission probabilities, of nuclei located near or on the r-process path 
has been obtained in some cases (for recent work, see Refs. |179l 11801 11811 11821 11831 I184J ) . 
A notable measurement of the Qp (/3-decay end point) value of the waiting point nuclide 
^■^°Cd is described in Ref. |185j . There is reasonable hope that many crucial nuclear 
properties of very neutron-rich, unstable nuclides can be measured in the near future at 
radioactive-ion beam facilities. 

Although the neutrino-driven wind represents the most popular r-process site, 
recent long-term simulations of core-collapse supernovae predict unfavorable conditions 
(electron mole fractions, entropies, expansion time scales) for an r-process [1861 I208J . In 
particular, these studies obtained proton-rich {Yf. > 0.5) neutrino-driven winds for an 
extended time period up to 20 s. Clearly, this represents a severe difficulty for r-process 
scenarios in core-collapse supernovae. Thus the quest for the site(s) of the r-process 
remains an unsolved puzzle in nuclear astrophysics. An alternative r-process scenario 
will be discussed in Sec. [THl For additional information on topics related to the r-process, 
see Ref. [20511209] . 

6.10. The gamma-ray burst-supernova connection 

Gamma-ray bursts (GRBs) are brief, intense flashes of electromagnetic radiation with 
typical photon energies (~100 keV) in the 7-ray domain. Indeed, they constitute the 
most luminous electromagnetic events in the Universe. This initial 7 flash is usually 
followed by a longer-lived afterglow emitted at longer wavelengths (X-rays, ultraviolet, 
optical, infrared and radio). Gamma-ray bursts are isotropically distributed on the sky 
and hit the Earth from unpredictable locations several times per day [210] . To date, no 
repeated burst has been detected from the same source. 

Gamma-ray bursts were serendipitously discovered in 1967 by the military U.S. 
Vela satellites, originally designed to detect secret nuclear weapon tests (suspectedly 
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carried by the former USSR, after signing the Nuclear Test Ban Treaty in 1963). On 
July 2, 1967, the Vela 3 and 4 satellites detected a flash of 7 radiation unlike any known 
nuclear weapon imprint. As other Vela satellites were launched, additional 7-ray flashes 
were detected. The analysis of the different arrival times of the bursts, performed by 
researchers from Los Alamos National Laboratory led by Klebesadel, ruled out either a 
terrestrial or a Solar System origin. The work was later declassified and revealed to the 
scientific community |211] . 

In contrast with many other stellar variables and transients, which exhibit a 
characteristic time evolution, with a rapid rise to peak luminosity followed by an 
exponential-like decay, 7-ray burst light curves show very complex and diverse patterns. 
Indeed, no two 7-ray burst light curves look alike. Some 7-ray bursts are preceded 
by a weak burst (precursor) that anticipates the true bursting episode by seconds or 
minutes. The light curve spans for an overall duration ranging from milliseconds to tens 
of minutes, being single peaked or showing several individual minipulses. Moreover, the 
peaks themselves are sometimes symmetric or characterized by a fast rise and a very 
slow fading. In spite of the large differences found in the 7-ray burst light curves (also 
in the spectra), two broad categories have been established on the basis of the burst 
duration and its power spectrum: short (hard) 7-ray bursts, with a typical duration of 
0.3 s, and long (soft) bursts with a duration of about 20-30 s. Moreover, the former are 
characterized by a 50% larger average peak energy (360 keV) than the long 7-ray burst 
class (220 keV) pT2] . 

More than hundred theoretical models aimed at explaining the nature of 7-ray 
bursts have been proposed since their discovery (including collisions between comets 
and neutron stars |213] ). Extensive reviews on the different mechanisms proposed for 
the engines of 7-ray bursts can be found in Refs. |214[ I215[ 1216] . Current models require 
the acceleration of tiny amounts of matter to ultrarelativistic speeds |214] . as well as 
a beamed emission focalized into a very small fraction of the sky. Once this beaming 
effect is taken into account, the total energy emitted in 7-rays (at least, in the long- 
lasting 7-ray bursts) is inferred to be ~ 10^^ erg, roughly equal to the overall kinetic 
energy in core-collapse supernovae. Despite this similarity, no evident link between 
the two phenomena was established for decades. An exception can be found in Ref. 
|217] . perhaps the only work that anticipated the existence of 7-ray bursts before their 
discovery, in the context of relativistic shocks from supernovae. Actually, the most 
widely accepted scenario for 7-ray bursts in the 1990s involved the merging of degenerate 
binary systems, such as double neutron stars |218[ 1219] or binaries containing a neutron 
star and black hole |220t 1221] . At that time, nothing was pointing towards a possible 
connection between 7-ray bursts and supernovae. An important achievement in this 
direction was the detection of the first X-ray and optical afterglows in 1997, as well 
as the direct measurement of their redshifts using optical spectroscopy. This put the 
merger hypothesis aside, while clearly favoring an origin associated with the death of 
young, massive stars. The discovery of GRB 980425 in conjunction with one of the most 
unusual supernovae ever observed, SN 1998bw J222[I223] . represented a breakthrough in 
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the history of 7-ray bursts. This has been followed by the discovery of at least two other 
clearly connected GRB-SN pairs: GRB 030329 - SN 2003dh |22l] and GRB 031203 - SN 
20031w [225]. In such systems, the detected supernova seem to correspond to the type 
Ic class, that is, supernovae involving massive progenitors, showing no hydrogen in the 
spectra and also lacking strong lines of He I or Si II. Even though some models suggest 
that all long 7-ray bursts might be accompanied by supernova explosions, it is clear 
that not all supernovae (not even the type Ic class) give rise to 7-ray bursts. The reason 
why massive stars choose one path (a SN-GRB pair) or another (a SN explosion only) 
is not yet clear. It has been claimed, however, that 7-ray bursts may result only from 
the most rapidly rotating and massive stars. It is worth noting that in the cases of SNe 
2003dh and 1998bw, the explosion yielded 10^^ erg of kinetic energy into a solid angle 
of ~ 1 radian. This implies about ~ 10 times more energy than a regular supernova. 
Such energetic explosions have been unofficially coined as hypernovae |226] . 

About 30% of the sample of 7-ray bursts detected by the BATSE instrument on- 
board NASA's Compton Gamma-Ray Observatory (CGRO) correspond to the short 
GRB class. However, there are reasons to believe that because of observational bias, 
this may be the most frequent variety of 7-ray bursts in the cosmos. Indeed, short 7-ray 
bursts tend to release smaller amounts of energy than long bursts, and this may limit 
their detectability beyond a certain distance. Moreover, short 7-ray bursts for which 
well-localized counterparts have been determined, are associated with elliptical galaxies. 
This clearly points towards old progenitors, for example, the merging of binary neutron 
stars, as the likely sources powering short 7-ray bursts. It has also been claimed that 
a small subset of this GRB class is probably associated with giant flares from soft-7 
repeaters, that is, sources of short (~ 100 ms), repeating bursts of soft 7-radiation 
(< 100 keV), located in nearby galaxies. 

Some of the models aimed to explain long 7-ray bursts involve the rotation (near 
breakup) of a highly magnetized neutron star (the magnetar model |227l 12281 1229] ). 
or the collapse of the innermost layers of a massive star forming a rapidly rotating 
black hole surrounded by a disk (the collapsar model |230l 1231] ). It has been argued, 
however, that all magnetar models computed so far systematically disregarded the role 
played by mass accretion (at a rate ~ 0.1 M© yr~^) onto the proto-neutron star before 
its rotation is fully established. It has also been suggested that these models may require 
the initial operation of a successful neutrino-powered explosion before they can function. 
Provided that a disk and a black hole form, the greatest uncertainty associated with the 
latter models is the mechanism for turning disk binding energy (or black-hole rotation 
energy) into confined relativistic outfiows. Suggested mechanisms include neutrinos 
|230] . magnetic disk instabilities [232], or a magnetohydrodynamic extraction of the 
black hole rotational energy [233]. 

The nucleosynthesis contribution of 7-ray bursts has been explored by different 
groups. In the framework of the collapsar model, two possible explanations have been 
given for the synthesis of Ni that likely powers the light curves. In the first one, as in a 
regular supernova, it may be linked to a strong shock traversing the stellar mantle and 
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exploding the star. Nomoto and collaborators |234[ 1235] discussed the nucleosynthesis 
consequences of this scenario, showing that a very energetic shock, driven by bipolar 
jets in hypernova explosions, may yield a peculiar abundance pattern resembling that 
inferred from some extremely metal-poor stars (in particular, C-rich metal-poor stars, 
and hyper metal-poor stars). Notice, however, that a recent reanalysis performed in 3-D 
[236] revealed a lower ^^Ni production than in previous work. The second possibility 
involves Ni production in the wind blown off from the accretion disk surrounding the 
black hole |231l 12371 1238] . It has been found that matter ejected from the disk can 
exhibit a rather rich range of nucleosynthesis, depending on the specific mass accretion 
rate of the disk, the velocity and entropy of the outflow, and the radius at which 
matter is released |239] . Hence, the corresponding outflow can be either neutron- or 
proton-rich. This has been analyzed in detail in a number of papers |239t 12401 1241] , 
showing that the outflows from rapidly accreting disks (M ~ 10 M0 s~^) may produce 
r-process elements (Sec. 16. 9p . Moderate mass-accretion rates can yield ^^Ni and a suite 
of light elements, such as Li or B. For disks with low mass- accretion rates (M ~ 1 
Mq s~^), the main products are likely ^®Ni, ^^Ti, ^^Sc, ®^Zn and ^^Mo. The levels at 
which those elements are produced make 7-ray bursts likely contributors to the Galactic 
inventory. A recent reanalysis |242] has shown that disk models with mass accretion 
rates near 1 M© s~^ may result, for a certain combination of parameters, in i/p-process 
nucleosynthesis products (Sec. 16. 8p . Nevertheless, the study revealed that the time 
during which matter is subjected to neutrino interactions turns out to be critical, and 
may even halt i/p-process nucleosynthesis. In order to shed light into this issue, state- 
of-the-art hydrodynamic simulations will be needed. 

7. Classical and recurrent novae 

Novae have captivated the interest of astronomers since more than two Millennia. 
Indeed, their own name, from the latin stella nova (or "new star"), reflects the reasons 
that lay behind such interest: the sudden appearance of a luminous object in the sky 
(at a spot where nothing was clearly visible before), that fades away, back to darkness, 
in a matter of days to months. 

Until the astronomical scale distance was soundly established, (classical) novae 
and (type I/II) supernovae were both associated with the same celestial phenomenon, 
generically designated as "nova" . Shapley and Curtis were among the first astronomers 
to question the real distance to spiral "nebulae" , where novae had been serendipitously 
discovered, suggesting their extragalactic nature. A major step forward was achieved 
with the discovery and analysis of Nova S Andromedae (1885) by Hartwig, as well as by 
the new scale distance of Hubble that placed Andromeda outside the Milky Way. This 
pushed Nova S Andromedae far away and revealed that the underlying stellar object was 
characterized by a huge intrinsic luminosity. Thus the road was paved for the hypothesis 
that two distinct classes of stellae novae exist, which are in modern parlance referred to 
as novae and supernovae. 
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The understanding of the physics that gives rise to the sudden luminosity burst 
characterizing novae has motivated vivid discussions. Probably the first explanation 
that underscored the real mechanism powering nova explosions appears in Newton's 
Principia Mathematica: "So fixed stars, that have been gradually wasted by the light 
and vapors emitted from them for a long time, may be recruited by comets that fall upon 
them; and from this fresh supply of new fuel those old stars, acquiring new splendor, may 
pass for new stars. " Indeed, the concept of revitalization of old stars by fresh supply 
of new fuel (although not by comets) is at the base of the thermonuclear runaway 
model, in which mass accretion plays a central role. But the basic picture that define 
nova outbursts emerged from a number of observational breakthroughs in spectroscopy 
(the discovery of neon lines in Nova Persei 1901 |243i 1244] . revealing the existence of 
several nova types) and photometry (the interpretation of the dips in the optical light 
curves as caused by dust forming episodes |245] , or the discovery of the binary nature 
of nova systems [246[ 1247] ) . All in all, even though the observed features were soon 
interpreted as produced by the ejection of a shell from a star |248j . the explanation 
of the mechanism that causes the burst had to wait until mid 20th century. Indeed, 
its thermonuclear origin was first theorized by Schatzmann |249l 1250] . and Cameron 
[251] (see also Gurevitch & Lebedinsky [252] . and references therein), while the first 
hydrodynamic simulation of a nova outburst was performed by Sparks |253] . 

Novae have been observed in all wavelengths, except in 7-rayq5i- They constitute a 
very common phenomenon, representing the second most frequent type of thermonuclear 
explosion in the Galaxy after type I X-ray bursts. Although only a handful of novae (3 
to 5) are discovered every year mainly by amateur astronomers, a Galactic rate of about 
30 ± 10 yr~^ [255] has been predicted through extrapolation of the number of events 
observed in other galaxies. The reason for the scarcity of detections in our Galaxy is 
interstellar extinction by cosmic dust. 

Classical novae take place in close stellar binary systems, with orbital periods in 
the range 1-12 hr. They consist of a white dwarf star and a low-mass main sequence 
companion (a K-M dwarf, although evidence of more evolved companions also exists). 
Contrary to type la supernovae (Sec. [8]), in which the white dwarf is fully disrupted by 
the violence of the explosion, a nova outburst is basically confined to the outer layers. 
The star is not destroyed by the explosion and hence all classical novae are expected 
to recur, with periodicities of the order of 10^ — 10^ yr. Notice, however, that in the 
explosions that occur in very massive white dwarfs, the so-called recurrent novae, the 
expected periodicities range typically between 10-100 yr. Both novae and supernovae are 
characterized by a remarkable energy output, with peak luminosities reaching 10^ and 
10^*^ L0, respectively. These explosive phenomena also differ greatly in the amount of 
mass ejected (the whole star in a thermonuclear supernova compared to only 10~^ — 10~^ 
Mq in a nova) and in the mean ejection velocity (~ 10^ km s~^ in a supernova compared 

^ An unexpected emission at >100 MeV detected in the symbiotic binary V407 Cygni, attributed to 
shock acceleration in the ejected nova shell after interaction with the wind of the red giant, was reported 
by the Fermi-LAT Collaboration [254] . 
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to ~ 10'^ km s~^ in a classical nova). Interested readers are referred to some of the recent 
monographs and reviews on this subject |256[ I257[ I258[ 12591 I265i 1266] . 

Under nova conditions, the early evolution of the burst is dominated by 
the operation of both the proton-proton chains and the cold CNO cycle (mainly 
through ^^C(p,7)^^N(/3"'"t/)^^C(p,7)^^N). As the temperature increases, the characteristic 
timescale for proton captures onto ^^N becomes shorter than its /3-decay time, favoring 
a number of reactions of the hot CNO-cycle, such as ^'^N(p,7)^^0, together with 
^^N(p,7)^^0, or ^^0(p,7)^^F. Since such nuclear processes take place in an envelope 
that operates in degenerate conditions, with a pressure depending only on the density 
rather than on temperature, the star cannot react to the temperature increase with an 
expansion. This paves the road for a thermonuclear runaway (TNR), in which convection 
plays a critical role. Actually, convective transport sets in when temperatures exceed 
~20 MK and carries a substantial fraction of the short-lived, /3+-unstable nuclei ^^'^^O, 
i7p^ i3]\j^ synthesized in the CNO-cycle, to the outer and cooler layers of the envelope, 
where they escape destruction by proton captures. In fact, it is the sudden release 
of energy from these short-lived species, which become the most abundant elements 
in the envelope after H and He, that powers the expansion and ejection stages in a 
nova outburst |267j . As a result, the ejecta will likely be enriched in their daughter 
nuclei ^^N, ^''O, and ^^C, and hence these species are expected to be the fingerprints of 
nova explosions in the chemical abundance pattern |268l 12691 12701 12711 1272] . This basic 
picture emphasizes the key role played by the initial ^^C content as the main trigger of 
the thermonuclear runaway. Indeed, the reaction ^^C(p,7)^^N determines the amount 
of mass accreted in the envelope and the pressure at the envelope's base, which in turn 
determines the strength of the outburst (i.e., peak temperature, mass and velocity of 
the ejected shells). Note that the rates for this reaction seem to be well established at 
all temperatures [27311271] . 

From the nuclear physics viewpoint, novae are unique stellar explosions. Their 
limited nuclear activity, which involves about a hundred relevant species in the A < 40 
mass range, linked by a few hundred nuclear reactions, as well as the limited range of 
temperatures achieved in such explosions (10-400 MK), allows us to rely primarily on 
experimental nuclear physics input [259] . In contrast to related astrophysical explosive 
sites that operate at much higher temperatures (e.g., supernovae. X-ray bursts), during 
a classical nova outburst the main nuclear path runs close to the valley of stability 
and is driven by (p,7) and (p,a) reactions, as well as by /^"'"-decays. Simulations also 
show that contributions from any neutron- or a-capture reactions, including ^^0{a, 
7)^^Ne, become negligible. Different studies have analyzed the role played by nuclear 
reaction rate uncertainties affecting the overall nova nucleosynthesis pattern [275] . 
This sparked an extraordinary activity in many nuclear physics laboratories worldwide. 
Notable measurements were performed for the reactions ^^0(p,7)^^F, using stable beam 
facilities [2601 12611 [3821 [263], and for the reaction '^^Na{p,jf^Mg, using the ISAC 
radioactive ion beam facility at TRIUMF, Vancouver [264] . The former process strongly 
impacts production of the radioisotope ^^F and the galactic synthesis of ^^O, while the 
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latter influences the production of 7-rays from ^^Na radioactive decay. Actually, the 
list of reactions whose uncertainty still has a strong impact on nova yields has been 
dramatically reduced. The main interest is now focused on measuring the challenging 
reactions ^8F(p,a)i50, 25Al(p,7)26Si, and 30p(p,^)3is. 

Models of nova nucleosynthesis point toward Ca as the likely nucleosynthesis 
endpoint, in agreement with observations. Despite the current problems associated with 
the modeling of nova outbursts (mainly, the discrepancy in the amount of mass ejected 
between models and observations [26811276] ). no large differences are found between the 
overall abundance pattern inferred from observations and that derived from numerical 
calculations. It is worth noting however that abundances are not directly determined but 
inferred, and their associated uncertainty is often quite large. Moreover, observations 
reveal that the ejecta accompanying classical novae are significantly enriched in metals: 
on average, CO-rich novae show a mean mass fraction of about 0.25, while ONe-rich 
novae are characterized by values near 0.5. Because of the moderate peak temperatures 
achieved in these explosions, it is unlikely that such metallicity enhancements could 
result from thermonuclear processing of solar-like material. Instead, mixing at the 
core-envelope interface has been proposed as the likely explanation. Several mixing 
mechanisms have been proposed, including multidimensional processes. It has to be 
emphasized, however, that all multidimensional nova simulations performed so far rely 
on extremely limited nuclear reaction networks, containing only 13 isotopes to properly 
account for the energetics of the explosion. These studies are also restricted to some 
specific stages of the thermonuclear runaway, following the development of convection 
and the onset of the explosion for several hundred seconds only. Bearing in mind that 
the overall process, from accretion to ejection, lasts for about ~ 10^ yr, it is clear 
that such multidimensional models cannot provide presently nucleosynthesis constraints. 
The only two-dimensional studies of mixing in novae available until recent years 
[277t [2781 [279| 12801 [28T] reached totally different conclusions regarding the feasibility 
of mixing by hydrodynamic instabilities. To solve the existing controversy, other 
multidimensional simulations have been performed recently with the hydrodynamic 
code FLASH. The new 2-D studies |282[ 1283] showed large convective eddies naturally 
developing close to the core-envelope interface through Kelvin- Helmholtz instabilities 
(instabilities resulting from a velocity shear between two fluids), with a size comparable 
to the height of the envelope (Fig. [9]), which dredge- up CO-rich material from the 
outermost layers of the underlying white dwarf into the accreted envelope. The mean 
metallicity achieved in the envelope, Z~0.20-0.30, is in agreement with observations 
of CO nova ejecta. Similar results have been found in a 3-D framework by the same 
group (in contrast to the results reported in Ref. [284] ). despite of the different way 
in which convection develops and extends throughout the envelope in the 2-D and 3-D 
simulations. 

Infrared and ultraviolet observations have significantly contributed to the 
understanding of the nova phenomenon, frequently revealing dust forming episodes in 
the ejecta |285] . Since the pioneering studies of dust formation in novae by Clayton & 
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Figure 9. Snapshots of the development of Kelvin- Helmhohz instabilities at the 
core-envelope interface in 2-D simulations of nova outbursts, at t = 215 s (upper left 
panel), 235 s (upper right), 279 s (lower left), and 498 s (lower right), shown in terms of 
^^C mass fraction along the computational domain. Notice how the Kelvin- Helmholtz 
instabilities dredge-up and transport unburnt CO-rich material from the outmost layers 
of the white dwarf core and inject it into the envelope. The mean CNO abundance in 
the envelope at the end of the simulations reaches ^0.30, by mass |282j . 



Hoyle |286] (a concept already suggested by Cameron in 1973 |287j ). all attempts devoted 
to the identification of possible nova grains relied mainly on the search for low ^''Ne/^^Ne 
ratios. Since noble gases, such as Ne, do not condense into grains, ^^Ne was attributed 
to in situ ^^Na decay, a clear imprint of a classical nova explosion. Indeed, Clayton and 
Hoyle pointed out several isotopic signatures (large overproduction of ^^'^''C, ^^O, ^^Na, 
^^Al or '^°Si), that may help to identify such putative nova grains. About fourty years 
later most, but not all, of these signatures still hold in view of our current understanding 
of such explosions (except for ^'^C, which is bypassed by the main nuclear path in novae, 
and ^^O, which is slightly overproduced in novae although grains nucleated in such 
environments are expected to be more anomalous in ^^O.). A major step forward in 
the identification of presolar nova grains was achieved by Amari et al. |288[ 1289] , who 
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reported on several SiC and graphite grains, isolated from the Murchison and Acfer 
094 meteorites, with abundance patterns qualitatively similar to those obtained from 
nova models: low ^^C/i^C and i^N/^^N ratios, high ^oSi/^^Si, and close-to-solar ^sSi/^^Si 
(plus high ^^Al/^'^Al and ^^Ne/^'^Ne ratios for some of the grains). However, in order 
to quantitatively match the grain data, one had to assume mixing between material 
synthesized in the outburst with more than ten times as much unprocessed, isotopically 
close-to-solar material before grain formation. One possible source of dilution might 
be mixing driven by the collision of the ejecta either with the accretion disk or the 
outer layers of the stellar companion |290] . Concerns about the likely nova paternity 
of the Amari et al. grains have been raised [291] . after three additional /im-sized SiC 
grains were isolated from the Murchison meteorite with similar trends, in particular, low 
i2q^i3q ^^^ ^^N/^^N ratios, but also non-solar Ti ratios. A supernova origin cannot 
therefore be excluded for at least some of the grains, although additional work on this 
subject seems mandatory to disentangle the controversy. Recently, other nova candidate 
(oxide) grains have been identified [292]. 

Among the species synthesized in classical nova outbursts, several radioactive nuclei 
have deserved particular attention as potential 7-ray emitters |293l I294[ I295[ 1296] . The 
species ^^N (T1/2 = 598 s) and ^^F (T1/2 = 1-83 hr) may provide a prompt 7-ray signature 
at 511 keV (through electron-positron annihilation), plus a continuum at lower energies 
due to Comptonized photons (with a cut-off at 20-30 keV because of photoelectric 
absorption) |297] . Unfortunately, this emission is short-lived, lasting only a few hours 
and takes place before the nova optical discovery. Hence, this signature must be detected 
a posteriori, through data mining by analyzing the data collected by any 7-ray satellite 
that was looking at the right spot on the sky at the right time. Even though the 511 keV 
line (and the lower energy continuum) has the largest predicted flux of all 7-ray signals, 
better chances of observing 7-rays from novae may be achieved with two other (long 
lasting) signatures: a 478-keV line from ^Be decay (T1/2 = 53.3 days), and a 1275-keV 
line from ^^Na decay (T1/2 = 2.6 yr). The species ^^Al is another important radioactive 
isotope that may be synthesized during nova outbursts, although only its cumulative 
emission can be observed because of its slow decay (T1/2 = 7.16 x 10^ yr). More details 
on the predicted 7-ray emission from novae can be found in Ref. |298] . 

8. Thermonuclear (type la) supernovae 

Type la supernovae (SN la) have become valuable cosmological tools. Through 
calibrated light curve analysis, they have been used as probes to outline the geometrical 
structure of the Universe, unraveling its unexpected acceleration stage |299t 1300] . 
Moreover, their leading role in the chemical and dynamical evolution of the Galaxy 
has motivated titanic efforts in the modeling of such stellar beacons. 
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8.1. Observational constraints 

Type la supernovae are observed in all types of galaxies, in sharp contrast with other 
supernovae types (SN lb, SN Ic or SN II), which are only found in spiral and irregular 
galaxies. This suggests that SN la are likely associated with old stellar populations and 
that the progenitors of SN Ib/c and SN II are much younger. 

Photometrically, type la supernovae exhibit a sudden rise in luminosity, up to a 
maximum visual absolute magnitudcJ of My ~ —19.3"^"^ (Lpeafe ~ 10^° L©) in about 20 
days. This phase is followed by a steep decline in brightness by about 3 magnitudes in 
~ 30 days, and later by a second, smoother decline over a period of ~70 days |30H 1302] . 

The energy required to power a supernova explosion can be inferred from 
determinations of the kinetic energy of the expanding ejecta (at typical velocities ranging 
between 5,000-10,000 km s~^), Ekin ~ 10^^ erg, as well as from the energy integrated 
over the light curve, Erad ~ 10"^^ erg. This energy budget obviously requires a powerful 
source. Early models attributed this to thermonuclear processing of CO-rich material 
into Fe-peak elements inside a white dwarf [303] . Indeed, the incineration of ~ 1 M© of 
a C-0 mixture, which releases ~ 10^^ erg g~^, accounts for the required 10^^ erg. 

During the late stages the light curve is powered by the radioactive decay chain 
^^m{/3+uf^Co{P+iyf^Fe ^M, [305], with the two different slopes attributed to the 
different half-lives of the corresponding nuclear decays, Ti/2{^^Ni) = 6.1 days, and 
^1/2(^^^*0) = 77.2 days. Estimates of the ^^Ni mass produced in type la supernovae 
can be determined empirically from modeling the late-time (nebular) spectrum and 
through application of Arnett's rule |306] : the peak luminosity is equal to the energy 
inputs from radioactive decays within the expanding ejecta. Assuming a rise time to 
peak luminosity of 19 days, Arnett's rule gives for the peak luminosity Lpeafc = 2 x 10^^ 
(Mtvj/Mq) erg s~^ [307] . Estimates for the ^^Ni masses from 17 type la supernovae, 
inferred from these methods, are in the range of 0.1-1 Mq [308] . It is worth noting that 
a large fraction of the energy released in the decays from ^^Ni and ^^Co is in the form 
of 7-rays, with energies of 158, 812, 750 and 480 keV for ^^m{/3+uf^Co, and 847 and 
1238 keV for ^^Co(/3'^z/)^®Fe. Hence, supernovae have become priority targets for many 
7-ray missions, like INTEGRAL, or in mission proposals, such as DUAL or GRIPS. 

Type la supernovae are spectroscopically identified by the absence of hydrogen 
emission Balmer lines and the presence of a prominent blue-shifted silicon II (P Cygni) 
emission feature near maximum light |302] . Whereas the first observational constraint 
poses limits on the maximum amount of hydrogen that can be present in the expanding 

"*" The apparent magnitude of a star, m, indicates its brightness as seen from Earth. It is based on a 
logarithmic measure of luminosity, such that a decrease in magnitude, Am, translates into an increase in 
Imninosity (from L to L') as L'/L = io^™/2-5. The absolute magnitude, M, is defined as the apparent 
magnitude that a star would have if it were located at a distance of 10 parsecs (32.62 light-years) 
from Earth. The connection between apparent and absolute magnitudes of a star can be expressed 
as m - M = 5 log (d/10), where d is the distance to the star in parsecs. The absolute magnitude 
and the luminosity of a star are related through M = Ms„„ - 2.5 log(L/LQ), where Ms„„ — 4.76 and 
L0 = 3.83 X lO'^'^ erg s^^ are the absolute magnitude and the luminosity of the Sun. 
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atmosphere of the star (i.e., Mh < 0.03 — O.IMq), the second feature reveals that nuclear 
processing must be taking place during the event. Spectroscopic snapshots at different 
epochs provide valuable clues on the physical mechanism powering such events: the 
spectra at maximum light are dominated by the presence of neutral or singly ionized 
lines of Si, Ca, Mg, S, and O, moving at high speeds (in the range of f ~ 8,000-30,000 
km s~^) |309] . Around two weeks after peak luminosity prominent, permitted Fe II lines 
are clearly seen. Finally, a month after maximum, the spectrum is dominated by the 
presence of forbidden emission lines from Fe II, Fe III and Co III. The overall picture 
clearly stresses key nuclear physics issues: the presence of intermediate-mass elements in 
the early spectra, when only emission from the outermost ejected layers is seen, reveals 
that the thermonuclear explosion did not incinerate the whole star (incomplete burning). 
This, in turn, provides clues on the flame propagation regime (a yet unsolved issue in 
supernova theory), since the presence of intermediate-mass elements rules out a pure 
detonation (Sec. 8.3). In sharp contrast, when inner regions of the star become later 
on accessible, the prominence of the Fe lines clearly points towards complete nuclear 
processing of matter to Fe-peak elements. Finally, the presence of Co lines at late stages 
strongly supports the hypothesis of a light curve tail powered by ^^ Co-decay. 

8.2. Progenitors 

The most widely accepted scenario for type la supernova explosions involves a carbon 
deflagration or detonation of a CO white dwarf that reaches its maximum mass (the so- 
called Chandrasekhar limit, Men ~ 1-4 M©) by accretion from a companion star. Helium 
white dwarfs have been ruled out as potential progenitors, since these objects would 
undergo extremely violent and too energetic detonations |310[ I311J , and no intermediate- 
mass elements would be produced. In turn, ONe(Mg) white dwarfs are expected to 
collapse, because of the efficiency of electron captures on ^"^Mg [3121 13131 13141 1315] . 

Two different progenitor systems have been proposed so far: a single degenerate 
one, dealing with a main sequence or giant companion that transfers matter onto a 
white dwarf star [316] : and a double degenerate system composed by two white dwarfs 
[3171 1318] . In the latter scenario, the explosion takes place when the two white dwarfs 
merge as a result of angular momentum loss driven by the emission of gravitational 
waves. During the last decades, both scenarios were regarded as competing, since SN 
la were considered a quite homogeneous class of objects, both at maximum light [319] 
and at late stages [320] . Hence, a unique progenitor was assumed. Arguments in favor 
of double-degenerates include the fact that those systems account in a natural way for 
the lack of hydrogen in the spectra. A major problem, however, involves the statistics 
derived from the search of systems with double-degenerates: only KPD 0422-1-5421 
[321] has a mass that could drive a Chandrasekhar-mass explosion after a merging 
episode lasting less than the age of the Universe. Although theoretical estimates based 
on population synthesis codes provide an appropriate frequency of massive mergers to 
account for the observed type la supernova rate [322], numerical simulations do not 
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necessarily support an explosive fate for such systems (see Refs. |323[ I324[ 1325] , and 
Sect. 10). These problems have revitalized the interest in single degenerate scenarios, 
with potential candidates like cataclysmic variables (formed by a white dwarf and a low- 
mass main sequence star) or symbiotic variables (a white dwarf plus a red giant star). 
In the first case, the accreted matter is H-rich. This poses problems since, even though 
some hydrogen should have been stripped from the secondary star during the explosion, 
this element is not seen in the spectra of type la supernovae |326l 1327] . But its major 
drawback is related to its path toward the Chandrasekhar mass: the theory of stellar 
evolution predicts that the maximum mass for a CO- white dwarf is about ~ 1.1 Mq 
(see Ref. |328] and references therein). If a CO-white dwarf is required for a successful 
explosion, such objects would need to accumulate roughly 0.3 Mq. This is difficult to 
achieve and, in particular, mass-accretion rates leading to nova-like explosions, implying 
loss rather than accumulation of mass, must be avoided (Sect. 7). In symbiotic variables, 
the accreted matter is instead He-rich. Several works |329[ 1330] have estimated a very 
high frequency of such systems in the Galaxy. Actually, to account for the observed 
Galactic SN la rate it would be sufficient if 1 — 4% of all symbiotic variables end their 
evolution with the accreting white dwarf reaching the Chandrasekhar limit. However, it 
has been pointed out |331] that white dwarfs in symbiotic binary systems scarcely reach 
the Chandrasekhar mass. Calculations yield about 0.1 Mq accreted during an assumed 
lifetime of ~ 5 x 10^ yr, suggesting that some white dwarfs may even explode before 
reaching the Chandrasekhar limit. 

Although the canonical scenario is based on the explosion of a Chandrasekhar 
mass CO white dwarf. He-rich accretion onto low-mass white dwarfs has also received a 
particular interest, mostly in connection with subluminous supernovae. In these models, 
helium ignites at the core-envelope interface, driving a double shock front (Fig. [TOj) : a 
detonation wave that moves outward, reaching the stellar surface after incinerating most 
of the envelope to ^^Ni, and an ingoing compressional wave that propagates toward the 
center, where it provokes carbon ignition |333] . This scenario has been suggested as 
a likely site for a number of nuclear processes, including ap-process, pn-process (i.e., a 
variant of the rp-process; Sec. |9]) and, more importantly, p-process (Sec. 16. 5p . capable of 
overproducing almost all p-nuclei in solar proportions, including the elusive Mo and 
Ru species (Refs. |334[ 1335] : see also Ref. |336] for recent work on the synthesis 
of p-nuclei in type la supernovae). The crucial issue turns out to be the precise 
He ignition density that determines the strength of the explosion, and therefore, the 
range of temperatures expected in the envelope. The overall amount of ^^Ni produced 
in these sub- Chandrasekhar models is often found to be smaller than in standard C- 
deflagration or detonation models, leading to fainter events and hence providing a 
possible origin for dim supernovae, like the sub-luminous SN 1991bg |333[ 1337] . From 
the nucleosynthesis viewpoint, these sub- Chandrasekhar models are characterized by 
a significant overproduction of a number of nuclear species, including ^^Ca (from the 
/3-decay of ^^Ti) and ^^Ti, which are systematically underproduced in type II and type 
la models. The main problem associated with this scenario is that only a very narrow 
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Figure 10. Temperature, density and velocity profiles resulting from He-rich accretion 
onto a 0.8 Mq CO white dwarf, at a rate of M = 3.5 x 10"^ Mq yr"\ following the 
progress of the detonation front originating at the envelope's base. The horizontal axis 
corresponds to the (internal) mass coordinate, ranging from (center of the star) to 1 
(surface) . Notice how the ingoing compressional wave that propagates inward provokes 
carbon ignition near the center (last panel). The simulation has been performed with 
the 1-D hydrodynamic code SHIVA 332j . 
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range of mass-accretion rates (~ 10~^ < Macci^e yr~^) < 5 x 10~^) gives rise to a 
central C-detonation and thus the models require some fine tuning. It is also worth 
noting that most of these studies have been performed in spherical symmetry and 
claims have been made that the quasi-central C-detonation is in fact imposed by the 
adopted geometry. Nevertheless, 3-D simulations of sub-Chandrasekhar models with 
asynchronous ignition (i.e., ignition takes place in a suite of different envelope points) 
produced successful supernova-hke explosions |338j . Nevertheless, the main arguments 
against sub-Chandrasekhar models (at least, as the likely scenario for the majority of 
type la supernova explosions) result from a number of spectroscopic anomalies, namely 
the prediction of a fast expanding layer consisting of ^^Ni and ^He, which is not observed 
in the spectra of SN la, and shorter than observed rise times to peak luminosity |339j . 
More recently, other supernova-like explosions, driven by He-mass transfer in AM 
CVn stellar binariecl have received some attention. Bildsten et al. |340j conclude that 
He-detonations producing primarily ^^Cr, ^^Fe, and ^^Ni are likely expected in some 
of these systems. They also suggest that these faint ".la" supernovae may represent 
between 1.6% and 4% of the local type la supernova rate. 

8.3. Type la supernovae as standard candles 

The increasing number of supernovae discovered has revealed a diversity among SN la, 
raising serious doubts on the uniqueness of the progenitor. About 85% of the observed 
SN la are believed to belong to a remarkably homogeneous class of events |341] , with a 
dispersion of only AM < 0.3™'^^ |342[ 1343] when normalized to peak luminosity. These 
SN la are known as "Branch-normals", with canonical examples such as SN 1972E, 
1981B, 1989B or 1994D. But the number of SN la deviating from this homogeneous 
class may be as high as ~30% |344[ 1345] . Already in the 1970s j319] . it was evident 
that two speed classes existed, a fast and a slow one, characterized by different rates of 
decline in their light curves. This diversity was soon interpreted in terms of explosion 
strength and envelope opacity: weaker explosions result in less luminous (and redder) 
events, showing a faster declining light curve and somewhat smaller velocities in the 
ejecta compared to more energetic supernovae |346] . Whereas SN 1991T is considered 
one of the most energetic supernovae ever observed, SN 1991bg and 1992K are the 
reddest, fastest and most subluminous SN la known to date |309[ 1347] . The difference 
in brightness between both groups amounts to ~ 2™"^. The existence of these different 
SN la classes was firmly established by Phillips |348] on the basis of an empirical relation 
between maximum brightness and rate of decline during the first 15 days after peak. 
This correlation is actually used to renormalize the observed peak magnitudes, thereby 
substantially reducing the dispersion in absolute magnitudes, and thus making type la 
supernova one of the most precise cosmological distance indicators. 

* AM CVn stars are a class of cataclysmic variables named after AM Canum Venaticorum. In these 
variables, a white dwarf accretes H-deficient matter from a compact stellar companion. They have 
extremely short orbital periods (<1 hr), and are predicted to be strong sources of gravitational waves. 



Nuclear astrophysics 56 

The variations in light curve (shape, peak luminosity, characteristic timescales, 
and so on) are interpreted as driven by a different amount of ^^Ni synthesized in the 
explosion, for the same overall ejected mass. Since peak luminosity is proportional to the 
^^Ni mass, brighter events are found when more ^^Ni is produced. This is accompanied 
by larger expansion velocities and broader peaks because of the increased opacity (deeply 
influenced by the presence of Fe-peak elements). Models based on light curve analysis of 
the dim SN 1991bg suggest that the amount of ^®Ni present in the ejecta was only 0.07 
Mq |349j ■ much smaller than the canonical value of ~ 0.6 M© associated with normal 
SN la [339] . 

8.4- Ignition regimes 

Carbon ignition in the degenerate interior of a white dwarf close to the Chandrasekhar 
limit can result in the propagation of a supersonic wave (a detonation driving nuclear 
burning by a shock wave) or a subsonic wave (a deflagration inducing burning by thermal 
conduction of the degenerate electrons). The actual outcome depends on the specific 
density, temperature, chemical composition and velocity profiles at the time ignition 
sets in (with low densities favoring the onset of detonations). In this section, we will 
briefly summarize the main characteristics of the different regimes proposed so far. More 
details can be found in Refs. [2661 I3501I351J . 

8.4-1 ■ The prompt detonation model. It assumes a pure detonation at the white dwarf 
center that propagates outward [352] . Although energetically feasible, it fails to account 
for the distribution of intermediate-mass elements inferred from the spectra of type la 
supernovae. Indeed, a detonation wave prevents the expansion of the layers ahead of 
the burning front and hence the fuel is almost totally incinerated into Fe-peak elements. 
The tiny amount of intermediate-mass elements produced is insufficient to provide the 
characteristic strong Si II feature, thus ruling out this model as the likely mechanism 
powering SN la. 

8.4.2. The deflagration model. When assuming a subsonic rather than a supersonic 
flame propagation, fuel layers can react to the advancing front with an expansion [353j . 
Thus, a fraction of the star's matter can in principle be processed into intermediate- 
mass elements. The feasibility of a deflagration model depends on its ability to avoid 
unburnt pockets of carbon and oxygen behind the front |354j . This has been explored in 
the framework of high-resolution, multipoint ignition simulations [355] , which show the 
inability to reach ^^Ni masses of >0.7 M© and kinetic energies exceeding 0.7 x 10^^ erg, 
and thus are unable to account for the bulk of Branch-normal type la supernovae. For 
the typical conditions that characterize white dwarf interiors, such deflagration fronts are 
highly subsonic [356] and hence allow for an expansion of the outer white dwarf layers. 
The expansion ultimately quenches the burning, at a time when the white dwarf is still 
gravitationally bound. Moreover, large amounts of unburnt C and O (> 0.57 Mq) are 
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left behind |357] . Two additional problems associated with deflagration models include 
the lack of chemical stratiflcation in the ejecta as well as the presence of big clumps of 
radioactive ^^Ni in the photosphere at peak luminosity [266]. 

Many multidimensional turbulent deflagration models have been published in recent 
years |358l 13591 13601 13611 136211363] . examining the effects of turbulence and instabilities 
on the propagation of the flame. 

8.4-3. The delayed detonation models. In this model, a deflagration front propagates 
and pre-expands the star. Subsequently, the deflagration wave switches into a detonation 
{deflagration-detonation transition). Several versions of this model exist, with [364] 
or without [3651 1366] assuming a pulsation. One-dimensional simulations of delayed 
detonations proved successful in reproducing many observational features of type la 
supernovae. Indeed, the combination of the two ignition regimes allows for the synthesis 
of intermediate-mass elements, while providing the required energy budget. Moreover, 
the model accounts for the expected light curves and photospheric expansion velocities 
[339] . The main uncertainty affecting these models deals with the mechanism driving 
the deflagration-detonation transition. Recent studies suggest that this naturally takes 
places once the flame reaches the so-called distributed burning regime [367] , at typical 
densities ~ 10^ g cm~^. This corresponds to the time when the laminar flame width 
equals the Gibson scale, that is, the scale at which turbulent velocity fluctuations equal 
the laminar flame speed. When the Gibson scale becomes smaller than the laminar 
flame width, turbulence affects the laminar flame structure. 

A somewhat related model, the gravitationally confined detonation model, has been 
recently proposed in the framework of multidimensional simulations 13681 . It relies on 
asymmetric deflagration flame ignitions, pushing burnt fuel towards the white dwarf 
surface. The collision of these ashes on the far side may generate a detonation wave 
propagating inward and burning the stellar core, potentially triggering very energetic 
supernovae 13691 370^ . This model also produces substantial amounts of intermediate- 



mass elements. Notice, however, that the feasibility of this scenario has been questioned 
by other studies {371^ . 



Alternative models based on pulsational delayed detonations have also been 
proposed to account for the bulk of type la supernovae. As in standard delayed 
detonation models, the initial deflagration wave pre-expands the star. Because of 
the slow (subsonic) flame velocity, the burning front is quenched and fails to unbind 
the star. But during recontraction, compressional heating at the interface between 
burnt and unburnt material ultimately triggers a detonation. Recent multidimensional 
deflagration models [358i 1359] . with a more accurate treatment of flame instabilities, 
suggest that the star may get unbound rather than recontracting. Finally, other models 
propose that when no detonation is triggered by the collision of ashes, a pulsational 
contraction may trigger a detonation wave that revitalizes the explosion [3721 13731 1374] . 
The reader is also referred to Refs. [3751 13761 13771 13781 1379] for additional information 
on multidimensional studies of delayed detonations. 
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Figure 11. Yields from 2-D models of type la supernova explosions: (left panel) 
carbon deflagration; (right panel) delayed detonation. Mass fractions are integrated 
over the whole ejecta, are shown relative to solar values, and are normalized to ^^Fe. 
Since delayed detonation models are initiated by a deflagration, the white dwarf 
can expand prior to the onset of the detonation. Therefore, the detonation moves 
through relatively low density layers producing only intermediate-mass elements (Fe- 
peak elements are only synthesized during the previous deflagration stage). As a result, 
the distribution of Fe-peak elements is very similar in both C deflagration and delayed 
detonation models. Figure from Maeda et al. [379], reproduced by permission of the 
American Astronomical Society. 



8.5. Nucleosynthesis: type la supernovae as Galactic Fe factories 

About once per century in our Galaxy, but probably once per second in the observable 
Universe, a white dwarf is totally disrupted by a titanic supernova explosion. The 
specific abundance pattern of the nuclear species freshly synthesized during such 
explosions depends critically on the peak temperature achieved, as well as on the neutron 
excess (Sec. 16. 3p . The latter depends in turn on the metallicity of the white dwarf 
progenitor, in particular the distribution of neutron-rich species (such as ^^Ne), as well 
as on the density at which the thermonuclear runaway occurs. In broad terms, the 
abundance pattern of the ejecta is the result of four different burning regimes: nuclear 
statistical equilibrium in the inner regions, incomplete Si-burning, 0-burning, and C/Ne- 
burning in the outermost layers (see Sec. 16.31 and Ref. |380] ). 

The impact of nuclear uncertainties on the nucleosynthesis produced in type la 
supernovae has not been analyzed in the same detail as for other astrophysical scenarios, 
such as classical novae or X-ray bursts. Indeed, no sensitivity study aimed at identifying 
the key nuclear reactions whose uncertainties have the largest effect on the yield has 
been performed to date. It is however worth mentioning the efforts made to accurately 
determine the cross section of the ^^C+^^C reaction, since it triggers the explosion of 
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these supernovae when temperature exceeds ~ 700 MK. A new low-energy resonance, 
with a strong impact on the overall rate, has been reported in recent years |381] . but 
follow-up studies have not confirmed its presence. Further studies aimed at improving 
our knowledge of this rate are clearly required because of its relevant role as a trigger 
of type la supernova explosions, and hence, on the ignition density. Other properties 
of the explosion can depend as well on the specific composition (more precisely on the 
^^C and ^^Ne content). This has been analyzed in detail in a number of papers (see 
Refs. [38 2 [ 13831 1384] . and references therein) that identified systematic influences of 
composition on supernova explosions, including the ignition density, the specific density 
at which the initial deflagration transforms into a detonation, the release of energy, or the 
flame speed, among others. Hence, reactions involved in the synthesis and destruction of 
such nuclear species are particularly relevant for our understanding of type la supernova 
explosions. 

Estimates of the frequency of type la supernovae in our Galaxy yield a value of 
~0.2-0.3 events per century [385]. This is similar to the expected occurrence of SN Ib/c 
(~0.65 per century), but smaller than the inferred rate for core-collapse supernovae 
(~3.32 per century) |386] . Since the amount of ^^Ni synthesized and ejected by each 
of the different supernova types is about 0.6, 0.3 and 0.07 Mq, respectively, it turns 
out that type la supernovae produce about half of the iron content in the Milky Way. 
Therefore, any model of type la supernovae must be able to reproduce the isotopic 
abundances of the Fe-peak elements determined in the Solar System (a task requiring 
Galactic chemical evolution models that take into account the different contributions 
from all possible Fe-peak producers). For decades, models systematically overproduced 
neutron-rich species, such as ^^Cr or ^°Ti, with respect to the Solar System values 
)387] 138811389] . This has greatly improved (but not yet solved) after the revision of the 
stellar weak interaction rates |390] , a key ingredient of all SN la nucleosynthesis models. 

For the different type la supernova ignition mechanisms discussed above, somewhat 
different thermal histories and thus different nucleosynthesis imprints are expected. This 
has been investigated in detail in the framework of spherically symmetric, hydrodynamic 
models, but only partially through multidimensional simulations. An early attempt, 
the famous carbon deflagration model W7 of Nomoto et al. |391] . was successful in 
providing the abundance pattern required by Galactic chemical evolution models, as 
well as in reproducing the basic features of the observed spectra and light curves of 
some "Branch-normal" SN la. The extension of this model to the delayed detonation 
regime improved the agreement with observations. 

Recent efforts to model the nucleosynthesis accompanying the explosions in a 
truly multidimensional framework have been carried out by a handful of groups 
[367] 1375] 1378] 1379] . An example of the expected nucleosynthesis has been recently 
reported in the framework of 2-D simulations by Maeda et al. j379] (Fig. [Till . We will 
briefly describe the differences reported between a pure turbulent deflagration model and 
a delayed detonation model. At the first stages of the deflagration, nuclear statistical 
equilibrium ensues. Electron captures become very important at the densities achieved. 



Nuclear astrophysics 60 

The most abundant species in this fully incinerated region are ^^'^^Fe, and ^^Ni (which 
later decays into ^^Fe). Because of the subsonic nature of deflagrations, the white 
dwarf reacts by expanding moderately, and hence the deflagration wave propagates 
through progressively lower density layers. Hence, the front reaches layers where nuclear 
statistical equilibrium still holds, but where electron captures are no longer relevant. 
The most abundant species in these layers are 56,58]\j-j_ rjj-^g front continues to advance 
through successively lower temperature and density regions, where matter undergoes 
0-burning or incomplete Si-burning, producing intermediate-mass elements, such as 
^*Si or ^^S. Carbon- or Ne-burning regions are encountered subsequently, resulting in 
large amounts of ^^O and ^^Mg. At some point, the low temperatures prevent further 
thermonuclear reactions. 

Similarly, in the delayed detonation framework, the detonation is ignited after 
significant expansion of the white dwarf has taken place, thus reducing the impact 
of electron captures. As the wave propagates outward, where the density decreases, the 
resulting nucleosynthesis is driven by oxygen burning, carbon burning, and eventually 
burning ceases. 

In terms of chemical stratification, the C deflagration model exhibits a broad region 
in which electron capture products (^^'^^Fe, ^^Ni), complete-Si burning products (^^Ni), 
intermediate-mass elements ('^^S, ^^Si, ^'^Mg), and unburnt elements (^^O, ^^C) co-exist. 
This is surrounded by a shell with unburnt material in regions that are never reached 
by the flame. The delayed detonation model results in a similar structure, but in the 
innermost regions unburnt C and O are efficiently processed into intermediate-mass 
elements by the delayed detonation wave. In the particular model reported by Maeda 
et al. [379], C is almost completely burned, but O remains with a mean mass fraction 
of ~ 0.1. Surrounding the deflagration region, there is a layer containing products of O 
burning, namely ^^Si, ^^S, and ^^Mg. This, in turn, is surrounded by a shell containing 
a mixture of C-burning products and unburnt material. To summarize, the delayed 
detonation model reproduces the characteristic layered structure, in contrast to what is 
found in the C deflagration model. 

It is finally worth noting that a layered structure, with a surface dominated by 
intermediate-mass elements in agreement with observations, has also been obtained in 
other delayed detonation models, such as the 3-D simulations of Ropke & Niemeyer 
[367] . In contrast. Bravo & Garcia-Senz |378j (Fig. [12]) found large amounts of Fe- 
peak elements near the surface regions and a structure that is not clearly layered. 
This abundance stratification has also been obtained in the 3-D delayed detonation 
models of Gamezo et al. |375j . The reasons for the reported differences between 
Bravo & Garcia-Senz and the other groups are yet unclear, although they have 
been attributed to differences in the modeling of the detonation, which is poorly 
understood in multidimensional frameworks. Extensive multidimensional simulations 
are needed to clarify these issues. The interested reader is referred to some of 
the existing monographs and reviews on type la supernovae, for further information 
[26611350113511 [3921I3931I39I1 1395] . 
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Figure 12. Snapshots of the evolution of a type la supernova explosion, in the 
framework of a 3-D delayed detonation model. Panels show a temperature map 
during the detonation regime, at times t — 0.8, 0.9, 1.0, and 1.1 s. Regions in red 
are characterized by temperatures of about 5 GK, while blue regions achieve at most 
0.5 GK. Credit: E. Bravo and D. Garcia-Senz, A&A, vol. 478, p. 849, 2008 [378] . 
reproduced with permission © ESO. 



9. X-ray bursts. The ap— and rp— processes 



Since the launch of NASA's X-ray sateUite Uhuru, in December 1970, astronomers have 
opened a new window to observe violent phenomena that occur in the cosmos. Indeed, a 
suite of different astrophysical sources have revealed an intense X-ray emission, including 
AGNs (active galactic nuclei powered by black holes), supernova remnants, and binary 
stellar systems hosting compact objects, such as white dwarfs (cataclysmic variables; 
super soft X-ray sources), neutron stars (X-ray bursters) or black holes (X-ray novae). 
Among these objects. X-ray bursters (XRBs), serendipitously discovered in the 1970s, 
have deserved particular attention. 
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XRBs exhibit brief, recurrent bursts. Their hght curves show a large variety of 
shapes (with single, double or triple-peaked bursts) and are characterized by a fast rise 
(~ 1 — 10 s) toward a peak luminosity of ~ 10^^ erg s~^, followed by a slower (sometimes 
exponential-like) decline. Recurrence periods range from hours to days |396ll397t 1398] . 

Even though the bursting episodes detected in the X-ray source 3U 1820-30 by 
the Astronomical Netherlands Satellite (1974-1976) |399j are often referred to as the 
first X-ray bursts ever observed, similar events were previously recorded by the Vela- 
5b satellite (1969-1970) |400[ 1401] during a survey in the constellation Norma. These 
early discoveries were soon followed by the identification of additional bursting sources, 
one of which, the enigmatic Rapid Burster (MXB 1730-335 |402j ). is characterized by 
puzzling recurrence times that are as short as ~10 s. To date, 92 Galactic X-ray bursting 
sources have been identified [403] . The first extragalactic X-ray burst sources have been 
discovered in two globular cluster candidate sources of the Andromeda galaxy [404] . 
Plans to extend the survey to other galaxies within the Local Group are currently 
underway. 

The thermonuclear nature of X-ray bursts was suggested shortly after their 
discovery. Maraschi & Cavaliere |405] . and independently Woosley & Taam |406] . were 
the first to propose that X-ray bursts are driven by thermonuclear runaways that occur 
on the surface of accreting neutron stars. Since the fast succession of bursts exhibited by 
the Rapid Burster did not match the general behavior of most X-ray bursts, two different 
categories, type I and type II bursts, were established. The former are accretion-powered 
thermonuclear runaways, while the latter are most likely associated with instabilities in 
the accretion disk. In the following, we will restrict the discussion to type I X-ray bursts, 
the most frequent kind of thermonuclear stellar explosion in the Galaxy, and the third 
most energetic (after supernovae and classical novae). 

The spatial distribution of type I X-ray bursts shows a clear concentration toward 
the Galactic center, and the sources are frequently located inside globular clusters. 
This observational pattern points toward old progenitors |396] . Indeed, the donors that 
transfer material onto the neutron star surface in X-ray bursts are known to be faint, 
low-mass stars (with M < 1 Mq), corresponding either to main sequence or red giant 
stars. The nature of the compact object hosting the explosion has been particularly 
controversial. Grindlay & Gursky [407] proposed a pioneering model of X-ray bursts that 
assumed accretion onto massive black holes (> 100 Mq). But shortly after, observations 
of bursting sources belonging to globular clusters |408] . from which accurate distance 
determinations could be inferred, suggested instead a much smaller object, either a 
neutron star or a stellar mass black hole. A number of spectroscopic features, together 
with estimates of masses inferred for those systems, pointed clearly toward a neutron star 
[409] . It has been suggested that time- resolved spectroscopy may allow measurements of 
the surface gravitational redshift [4101 1411] . from which constraints on the mass-radius 
relationship could be derived. 

The orbital periods determined in X-ray burst sources range typically between 1 - 
15 hours |412] . Hence, these binary systems are so close that mass transfer episodes. 
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driven by Roche-lobe overflows, ensue. This leads to the formation of an accretion 
disk surrounding the compact star. The maximum mass-accretion rate is set by the 
Eddington limit (M^^^ ~ 2 x 10~^ Mq yr~^, for H-rich accretion onto a 1.4 Mq neutron 
star). It was soon realized that the ratio of the gravitational potential energy released 
by the matter falling onto the neutron star during the accretion stage (G MNsf^NS ~ 
200 MeV per nucleon) and the nuclear energy released during the burst (~ 5 MeV per 
nucleon, for a fuel of solar composition), matches the observed ratio of time- integrated 
persistent and burst fluxes (~ 40-100). 

Progress in the understanding of the physics of X-ray bursts has been achieved 
through multi-wavelength observations beyond the X-ray domain. In 1978, the flrst 
simultaneous optical/X-ray emission was detected from the X-ray burst source 1735-444 
|413j . The flux in the optical was interpreted as the low-energy tail of the blackbody 
X-ray burst emission [396J. But more importantly, the optical burst was delayed by ~3 
seconds with respect to the burst observed in X-rays [414] . A similar delay has been 
observed in many other X-ray burst sources. It has been suggested that this optical 
emission is caused by reprocessing of X-rays in material within a few light-seconds 
from the source, either at the accretion disk that surrounds the neutron star or at the 
hemisphere of the secondary, which is directly illuminated by the source. In other words, 
the delay in the optical emission could be explained by travel-time differences between 
the X-rays directly emitted by the source, and those that hit the disk or the secondary 
and subsequently are scattered down in energy and become optical photons. Although 
efforts have been invested at other wavelengths, neither infrared nor radio emission has 
been unambiguously detected simultaneously with X-rays from a bursting source. 

9.1. Models of type I X-ray bursts 

Van Horn & Hansen |415l 1416] flrst pointed out that nuclear burning on the surface of 
neutron stars can be unstable. The different regimes for unstable burning have been 
extensively discussed elsewhere |398[ I417[ 1418] . For a chemical mixture with a CNO mass 
fraction of ~ 0.01, for instance, mixed H/He-burning is expected for mass accretion rates 
of M< 2 X 10^^° Mq yr~^, which is triggered by thermally unstable H ignition. Pure 
He-shell ignition results for mass accretion rates in the range from 2 x 10^^*^ M© yr^^ 
to 4.4 X 10~^° Mq yr~^, following completion of H-burning. Mixed H/He-burning is 
expected for M> 4.4 x 10^^*^ Mq yr^^, this time triggered by thermally unstable He 
ignition. 

Since the pioneering work of Maraschi & Cavaliere |405j . Woosley & Taam |406j . 
and Joss |419] , many groups have performed simulations of the dynamics of X-ray bursts 
and their associated nucleosynthesis. With a neutron star as the underlying compact 
object hosting the explosion, temperatures and densities in the accreted envelope 
reach rather high values (Tpeak ^ 10^ K, p^ax ~ 10^ g cm~^). As a result, detailed 
nucleosynthesis studies require the use of hundreds of nuclides, up to the SnSbTe-mass 
region [420] or beyond [421] , and thousands of nuclear interactions. Indeed, the endpoint 
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of the nucleosynthesis in X-ray bursts is still a matter of debate. For example, recent 
experimental work suggests that it will be more difficult to reach the SnSbTe-mass region 
[422] . Because of computational limitations, studies of X-ray burst nucleosynthesis have 
usually been performed with limited nuclear reaction networks that are truncated near 
Ni 11231 MM 1125], Kr |l26l \m\, Cd |12H], or Y |129]. More recently, Schatz et al. 
|420[ 1430] carried out detailed nucleosynthesis calculations with a network containing 
more than 600 nuclides (up to Xe), but assuming a single burning zone only. Koike et 
al. |421] also performed one-zone nucleosynthesis calculations, with temperature and 
density profiles obtained from a spherically symmetric evolutionary code, linked to a 
1270-nuclide network extending up to ^^^Bi. Extensive 1-D hydrodynamic calculations 
with detailed nuclear reaction networks (up to Te) have been published in recent years 
[43 It I432t 1433] . To date, no multidimensional calculation for realistic X-ray burst 
conditions has been successfully performed, although several attempts are currently 
underway with hydro codes such as MAESTRO |434] or FLASH. Notice that the 2-D 
study reported in Zingale et al. |435] analyzes He-detonations on neutron stars, whereas 
thermonuclear runaways are expected to propagate defiagratively in X-ray bursts. 

Detailed accounts of the main nuclear paths during type I X-ray bursts, both in 
the framework of one-zone models and hydrodynamic (1-D) simulations, have been 
given elsewhere [391 ESIl SSSl gSB |l32l gSSl HSZ]. To summarize. X-ray bursts are 
powered by a suite of different nuclear processes, including the rp-process (rapid proton- 
captures and /S^-decays), the 3a-reaction (Sec 15. 2p . and the ap-process (a sequence of 
(«,p) and (p,7) reactions). The main nuclear fiow proceeds far away from the valley of 
stability (Fig. [T3|), reaching the proton drip-hne beyond A = 38 [430]. Contrary to the 
thermonuclear runaways characteristic of other scenarios (e.g., classical novae). X-ray 
bursts are halted by fuel consumption instead of expansion when degeneracy is lifted. 
In fact, breakout from the CNO region via ^^0(a,7)^^Ne leads the path toward the 
synthesis of intermediate- mass and heavier species. During the final stages of the burst, 
the main nuclear fiow is dominated by a multitude of successive /3-decays that power the 
decline of the light curve. The final composition of the envelope depends on a number of 
parameters in the simulations (mainly, the mass-accretion rate and the adopted initial 
metallicity), but tends to favor species in the mass range A=60-70 (mainly ^^Zn and 
^^Ge [433] ). Other studies [432] find production of additional species with masses near 
A~105. 

A major inconvenience in the modeling of X-ray bursts comes from the lack of 
clear observational nucleosynthesis constraints. The potential impact of X-ray burst 
nucleosynthesis on the Galactic abundances is indeed a matter of debate. Although 
ejection of matter from a neutron star is very unlikely because of the extremely large 
escape velocities required, it has been proposed that late radiation-driven winds during 
photospheric radius expansion may lead to the ejection of a tiny fraction of envelope 
matter that contains processed material [438t 1439] . Although it has been claimed 
[437] that X-ray bursts may help to explain the Galactic abundances of the elusive 
light p-nuclei (Sec. 16. 5p . recent hydrodynamic calculations j432] have shown that the 
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Figure 13. Main nuclear activity shown in terms of the most abundant species at the 
innermost envelope shell, for two different XRB models. Color legends indicate different 
ranges of mass fractions (in logarithmic scale). Calculations have been performed with 
the code SHIVA [332] . and rely on 1.4 Mq neutron stars {Uni = 1.6 x 10^** erg s"^ 
= 4.14 L©), accreting mass a rate Mace = 1-8 x 10"^ M© yr"! (0.08 MEdd)- The 
composition of the accreted material is assumed to be solar-like (X = 0.7048, Y = 
0.2752, Z = 0.02) [upper panel], or metal-deficient (X = 0.759, Y = 0.240, and Z = 
0.001 ~ Z0/2O) [lower panel]. Both panels display the nuclear activity at Tpeak, for 
the first computed burst. 
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concentration of these p-nuclei in the outermost layers of a neutron star envelope is 
orders of magnitude below the required levels. It has also been claimed that the 
nuclear ashes produced by X-ray bursts may provide characteristic signatures, such 
as gravitationally redshifted atomic absorption lines that could be identified through 
high-resolution X-ray spectra, representing a valuable tool to constrain X-ray burst 
models |438i I440[ 144 ![ 1442] . However, this issue is controversial at present. Cottam et 
al. |443j reported features in the burst spectra of 28 X-ray bursts detected from the 
source EXO 0748-676, which were interpreted as gravitationally redshifted absorption 
lines of Fe XXVI (during the early phase of the bursts), Fe XXV, and perhaps O VIII 
(at late stages). But no evidence for such spectral features was found either during the 
16 X-ray bursts observed from GS 1826-24 |444j or from the bursting episodes detected 
from the original source |445t 1446] . 

Reliable estimates of the neutron star surface composition are nevertheless 
important since they determine a suite of different thermal |430l 1447] , radiative |448] . 
electrical |430l 1449] , and mechanical properties |450l 1451] of the neutron star. Also, the 
diversity of shapes of X-ray burst light curves is probably driven by different nuclear 
histories [i26lH52] . 

Finally, it is important to stress that most of the reaction rates adopted in these 
extensive nucleosynthesis calculations rely on theoretical estimates using statistical 
models of nuclear reactions and, therefore, may be subject to significant uncertainties. 
Efforts to quantify the impact of such nuclear physics uncertainties on the overall 
abundance pattern for X-ray burst nucleosynthesis have been performed by a number of 
groups Hai 11531 11291 HSU HSU 0331 0551 HSU SSS]. The most extensive works to date 
were performed by Parikh et al. |457l 1458] . who identified the key nuclear processes 
(about 50 nuclear reactions out of 3500, including ^^As(p,7)^^Se, ^^Ga(p,7)^^Ge, 
"'^^C(a;,7)^^0, and ^^Ag(p,7)^'^Cd), as well as key nuclear masses (in particular of ^^Ge, 
65 As, 66Se, 69Br, ^^Kr, ^^Nb, ^^Mo, ^6,87^^ 96 Ag, 97q^^ losg^^ ^^^ loeg^)^ „hose 

uncertainties deeply influence the final yields. 

9.2. Superbursts: the next frontier 

A few extremely energetic bursts have been reported in recent years as a result of better 
performances in monitoring achieved with X-ray satellites, such as Chandra or XMM- 
Newton. These rare and rather violent events are known as superbursts |459[I460] . They 
are characterized by long durations with typical exponential-like decay time ranging 
from 1 to 3 hours (including an extreme case, KS 1731-260, that lasted for more than 
10 hours [461] ). With a released energy of about ~1000 times that of a typical X-ray 
burst (~ 10^^ erg), superbursts exhibit much longer recurrence periods. For example, 
the period amounts to 4.7 yr for the system 4U 1636-53, for which two superbursts 
have been observed to date |462j . Although superb urst sources also exhibit type I X-ray 
bursts, their occurrence is quenched for about a month after each superb urst. The first 
superburst was reported by Cornelisse et al. [463] in the framework of a "common" 
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type I bursting source, the BeppoSAX source 4U1735-44. To date, 15 superbursts from 
10 different bursting sources fiave been discovered, including GX 17+2, for which 4 
superbursts have been identified |464] . 

The duration and energetics of superbursts suggest that they result from 
thermonuclear flashes occurring in fuel layers at greater depths than what is typically 
encountered for X-ray bursts (corresponding to densities exceeding 10^ g cm~^ |465j ) . 
Thus they are more likely to involve C-rich ashes resulting from prior type I X-ray bursts 
[ii9l 11651 M6\ \Wn MSi SM] , as first proposed by Woosley & Taam [106] . Unfortunately, 
it is not clear whether enough carbon is left after a type I burst |420[ I430[ I433[ 1432] . 
Nevertheless, Gumming & Bildsten |465] suggested that low concentrations of carbon 
are sufficient to power a superburst, especially in neutron star "oceans" enriched in 
heavy nuclei. Moreover, recent studies indicate that both stable and unstable burning 
of the accreted H/He mixture may be required to power a superburst |470] . 

10. Stellar mergers 

A small fraction of all stars, perhaps 10%, are expected to merge. Such merging episodes 
have been invoked to explain the peculiar properties of a suite of different objects, 
such as blue straggler stars (presumably products of two old, red stars merging in a 
dense stellar environment), R stars (suggested to result from the merger of a He white 
dwarf with a H-burning red giant branch star |471] : but see also Ref. |472] ) . HdG, 
R Gr B, and AM GVn stars (presumably produced in double white dwarf mergers 
|473l 14741 14751 I476[ 1477] ) . V838 Mon-type stars (mergers of a main sequence and a 
pre-main sequence star |478] ). and short 7-ray bursts (from double neutron star or 
neutron star-black hole mergers; Sect. 6.10)1^. Because of their importance as potential 
sources of gravitational waves, particular attention has been focused on the merging of 
compact objects |479[ 1480(1481] . Galactic compact binaries, such as double neutron stars, 
neutron star-white dwarf binaries, or close double white dwarfs, are prime targets for 
the planned space-borne gravitational wave interferometer LISA. In fact, the emission 
from Galactic (close) white dwarf binary systems is expected to be the dominant source 
of the background noise in the low frequency range. Simulations of the merging of a 
white dwarf and a neutron star |482] . and of double white dwarfs |483[ I484[ 1485] . have 
received little attention so far compared to the coalescence of double neutron stars (Ref. 
|487] . and references therein). 

tt Blue stragglers are main sequence stars that belong to open or globular clusters, and are more 
luminous and bluer than stars at the main sequence turnoff point for the cluster; R stars are K-type 
giants, enriched in 12,13q g^j^j i4-]\j. ^^l of these are single stars and standard stellar evolution theory 
cannot explain why they are C-rich; HdC stars are a type of H-deficient, supergiant C-rich star; R Cr 
B stars are H-poor, C- and He-rich, high-luminosity objects resembling simultaneously erupting and 
pulsating variables. Finally, V838 Monocerotis is a red variable star that experienced a major outburst 
in 2002. Originally believed to be a classical nova, V838 Mon exhibited an unusual light curve. The 
nature of the outburst is still not well understood. 
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10.1. Double white dwarf mergers 

The merging of two white dwarfs by emission of gravitational radiation can be considered 
the final fate of a significant fraction of this type of systems, with an estimated frequency 
of ~ 8.3 X 10-3 jT~\ 

Several hydrodynamic simulations of the coalescence of double white dwarf systems 
have been performed [483[ I485[ I486[ 1488] for different white dwarf masses and 



composition affl . Mass transfer ensues as soon as the less massive white dwarf fills its 
Roche lobe (Fig. [T^ . For certain combinations of masses of the double white dwarf 
system, mass transfer is stable |490j . In those cases, since the white dwarf radius scales 
as M~^/3, when the secondary loses mass its radius increases, and hence, its mass- loss 
rate would increase, reinforcing the overall mass transfer episode. An accretion "arm" 
forms, extending all the way from the secondary to the surface of the more massive 
white dwarf. In fact, the arm becomes entangled because of the orbital motion of the 
coalescing white dwarfs, giving rise to a spiral shape. Finally, the secondary becomes 
fully disrupted and a Keplerian disk surrounding the primary star forms after some 
orbital periods, depending on the initial conditions. 

Models involving a He white dwarf undergo partial thermonuclear fusion as the 
temperatures achieved during the merger exceed 10^ K. This translates into noticeable 
abundances of Ca, Mg, S, Si and Fe. In fact, the abundances of some of these species 
are even larger in the hot corona surrounding the primary star. This has been claimed 
as a feasible explanation for the origin of the recently discovered metal-rich DAZd white 
dwarfs, a subclass with strong H lines and dusty disks |482] . Indeed, the reason for 
such a metallicity enhancement could be related to the accretion of a minor planet that 
condensed in the metal-rich disks resulting from a merger episode. 

The expected gravitational wave emission from these systems would be 
characterized by an initial, almost sinusoidal, pattern of increasing frequency (the so- 
called chirping phase), followed by a sudden disappearance of the gravitational wave 
signal, on a timescale of the order of the orbital period as the merger proceeds. All such 
merging episodes could be detectable by LISA with different signal-to-noise ratios. 

10.2. Neutron star - white dwarf mergers 

Another interesting case of compact star mergers involves systems composed of a neutron 
star and a white dwarf. About 40 binaries of this type have been discovered to date 
[491] . although the number of systems expected to merge in less than the age of the 
Universe within a distance of < 3 kpc from the Sun could be as high as ~850 |492j . This 
results from the relatively high expected number of these binaries, with an estimated 
neutron star-white dwarf formation rate that is about 10-20 times larger than what is 
predicted for double neutron stars [493]. 

ft The nucleosynthesis in double white dwarf systems undergoing collisions or close encounters has also 
been modeled recently [489] . 
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Figure 14. Snapshots of the merging episode of two white dwarfs of masses 1.2 and 0.4 
Mq. Calculations have been performed with a smoothed particle hydrodynamics (SPH) 
code. Each panel shows the equatorial temperature contours as well as the location of 
the SPH particles. When the first particles of the secondary star hit the surface of the 
primary (top right panel), the shocked regions reach temperatures up to 5 x 10* K. 
The colliding particles are reflected almost instantaneously by the hard boundary of 
the primary. As mass transfer proceeds (central panels) the accreted matter is shocked 
and heated on top of the surface of the primary. This matter is actually ejected from 
the surface of the primary and a hot toroidal structure forms (right middle panel). 
Temperatures exceed 10^ K and nuclear reactions transform the accreted He into C 
and O. Credit: J. Guerrero, E. Garcia-Berro, and J. Isern, A&A, vol. 413, p. 268, 
2004 |483j ■ reproduced with permission © ESO. 
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A neutron star-white dwarf binary is expected |494] to form after an explosion 
of a massive primary star in a core-collapse supernova (Sec. [6]). For an appropriate 
range of primary star masses, the explosion would be accompanied by the formation 
of a neutron star. In its evolution toward a neutron star-white dwarf binary, the 
system must first avoid disruption by the supernova blast. This depends critically on a 
number of parameters, including the mass ratio between the stellar components. The 
final evolution is governed hereafter by the secondary, which undergoes severe mass-loss 
episodes during its red giant phase, and ends up as a white dwarf. 

Two classes of neutron star- white dwarf binaries can be distinguished |494[ 1495] : 
intermediate-mass binary pulsars, characterized by very eccentric orbits and masses of 
the secondary stars ranging between 0.5-1.1 Mq, and low-mass binary pulsars, with 
almost circular orbits and masses of the secondaries between 0.15 and 0.45 Mq. Both 
classes are expected to merge because of the emission of gravitational waves with an 
estimated frequency of 1.4 x 10~^ yr~^ |496j . Neutron star-white dwarf binaries expected 
to merge in less than the age of the Universe include the sources PSR 30751-1-1807, PSR 
31757-5322, and PSR 31141-6545. 

Simulations of these compact star mergers are scarce and despite their interest 
as potential nucleosynthesis sites, most of the efforts have focused on their role as 
future laboratories for gravitational wave studies. After coalescence, theses systems 
evolve into a configuration clearly resembling that discussed for a double white dwarf 
merger, with a a central compact object (a neutron star, in this case) surrounded by a 
heavy Keplerian, self-supported accretion disk |497] . The expected gravitational wave 
emission from those systems bears also a clear resemblance to that of a double white 
dwarf merger, with a short chirping phase detectable by the future LISA interferometer 
up to an estimated distance of about 250 kpc. Even though the simulations performed 
to date face severe limitations (e.g., no post-Newtonian corrections, and very limited 
nuclear reaction networks), some interesting results have already been obtained. In 
particular, it has been estimated that the mass ejected from these systems amounts to 
about 0.01% of the overall system mass assuming a merger of a 1.4 M0 neutron star and 
a 0.6 Mq white dwarf |497j . Further studies aimed at determining the nucleosynthesis 
impact of these mergers on the Galactic abundances are clearly needed. 

10.3. Double neutron star and neutron star-black hole mergers. Another r-process site? 

Mergers of stellar binaries containing double neutron stars, or a neutron star and a 
black hole, have been extensively studied as r-process nucleosynthesis sites, as potential 
sources of gravitational waves, and as the central engines for powering short 7-ray bursts. 
To date about ten stellar binaries that are likely containing double neutron stars 
have been identified |498] . Of these, about half are expected to merge over the age of the 
Universe. This is roughly in agreement with the predicted frequency, between 40 and 
700 Myr~^ |499l I500J . for a galaxy like ours. The frequency of neutron star-black hole 
mergers is highly uncertain, with estimates varying by orders of magnitude [500[ 1501] . 
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Nevertheless, with one of the stars being a black hole, these systems are expected to 
play a minor role in the chemical evolution of the Galaxy. 

Pioneering simulations of double neutron star and neutron star-black hole mergers 
were performed by Oohara & Nakamura [502] . assuming Newtonian gravity and 
polytropic equations of state (see also Refs. |503l 15041 1505] ). Several improvements, 
aimed at providing a more realistic description compatible with general relativity 
include the use of pseudo-relativistic potentials |506[ 1507] . as well as post-Newtonian 
corrections J5081 1509] . Another approximation adopted, halfway between Newtonian 
gravity and general relativity, is the conformal flatness approximation (CFA), in which 
the spatial part of the metric is assumed to be conformally flat |510l 15111 15121 1513] . The 
first fully general relativistic calculations of such mergers were performed by Shibata 
and coworkers |514] . It is worth noting that two complementary general relativistic 
formulations are presently used, the Baumg arte- Shapiro -Shibata- Nakamura (BSSN) 
formulation, and the generalized harmonic formalism (see Refs. |515[ I516[ 1517] . and 
references therein). More details on the different methods adopted in the modeling of 
double neutron star and neutron star-black hole mergers can be found in Ref. |487] . 

It has been suggested that the contribution of these mergers to the Galactic 
inventory of chemical species can proceed through three different processes: (i) direct 
ejection, (ii) neutrino-driven winds, and (iii) disk disintegration. Lattimer & Schramm 
[518] first provided estimates of the fraction of a neutron star that becomes unbound 
after a neutron star-black hole merger. Their calculations yield a mean ejected mass per 
merger of about 0.05 Mats, or equivalently, an accumulated value of 2 x 10""^ Mq pc~^, 
which is approximately equal to the mass fraction of r-process elements in the Solar 
System. More realistic simulations |507l 15191 1520] yield a mean ejected mass of 3 x 10^^ 
Mq, with a mean electron mole fraction in the range of 0.01 < Y^. < 0.5. Neutrino- 
driven winds emitted by the merger remnant are considered a promising alternative to 
the more traditional core-collapse supernovae as a production site for r-process nuclei 
(see Refs. |242[ I52H 1522] . and Sect. I6.10p . Simulations suggest that such neutrino- 
driven winds drive ejection of about ~ 10~^ Mq of matter with Y^ ~ 0.1 — 0.2 per event 
|523] . Finally, the disk resulting after the merging episode may give rise to the ejection 
of neutron- rich material, either through viscous heating in an advective disk or through 
the energy produced during recombination of nucleons into nuclei |524[ 1525] . It has been 
claimed that about 0.03 Mq of moderately neutron-rich material (0.1 <Y^ < 0.5) may 
be ejected from the disk |526] . 

First estimates of the impact of such double neutron star mergers on the Galactic 
nucleosynthesis |527] yielded a promising agreement with the observed r-process pattern 
beyond Ba. However, this was questioned by detailed inhomogeneous chemical evolution 
studies |528] on the basis of the low expected frequency of such mergers. The reasons 
included a lack of consistency between the theoretical predictions and observations at 
very low metallicities, and the fact that injection of r-process material produced by a 
single merger would lead to a scatter in the ratios of r-process and Fe abundances at 
times too long compared to observations. While core-collapse supernovae seem to be 
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favored as the main r-process site (Sec. 16. 9p . a contribution from neutron star mergers 
and the existence of multiple r-process sites still remain open questions. 

11. Non-stellar processes. Cosmic ray nucleosynthesis 

11.1. Quest for the origin of the light nuclides 

We have addressed so far the origin of most nuclides, except of the light species ^Li, 
^Be, ^°B and ^^B. Their solar abundances are smaller by about 6 orders of magnitude 
compared to other light nuclides (although their abundances are much larger compared 
to almost all s-, r-, or p-nuclides). As was the case for ^H and ^Li, their cross sections 
for proton-induced reactions are so large, because of the small Coulomb barriers, that 
they are destroyed during the hydrogen burning phase in stellar interiors already at 
temperatures below a few million kelvin. In fact, under such conditions their mean 
lifetimes amount to less than a few Gyr. On the other hand, standard big bang 
nucleosynthesis (Sec. [3]) produces only negligible amounts of ^Li, ^Be, ^''B and ^^B, 
with predicted number abundances less than 4 orders of magnitude compared to ^Li 
|529j . Their origin seemed so obscure that they were attributed in the seminal work of 
Burbidge, Burbidge, Fowler and Hoyle [5] to some unknown "x process". Among the 
suggestions made in that work was the production by spallation reactions involving 
high-energy protons, neutrons or a-particles of energies in excess of 100 MeV per 
nucleon on abundant CNO nuclei in stellar atmospheres. While it was convincingly 
demonstrated that the energy available during the T Tauri phase of young stars is 
insufficient to produce the light nuclides in stellar atmospheres |530] . the association of 
lithium, beryllium and boron synthesis with spallation reactions proved correct. 

An important piece of evidence was uncovered around 1970, when it was pointed out 
that in the Solar System the ratio of number abundances for Li, Be and B compared to 
C, N and O amounts to ~10~^, while the value in Galactic cosmic rays is near ~0.2 |531j . 
Consequently, the idea was born that the source of Li, Be and B are spallation reactions 
between Galactic cosmic rays and nuclei in the interstellar medium. The conjecture 
is supported by the observation that the spallation cross sections of protons and a- 
particles on C, N and O nuclei at energies in excess of 100 MeV per nucleon favor the 
production, in decreasing order, of B, Li and Be |532] . The observed Galactic cosmic ray 
abundances of these light species exhibit the same exact ordering. More information 
on abundances and associated spallation cross sections is provided in Fig. [151 First 
quantitative models |534] . based on measured spallation cross sections, showed that this 
scenario could reasonably account for the abundances of ^Li, ^Be and ^°B observed in 
the Solar System and in stars after 10 Gyr of Galactic chemical evolution. 

On the other hand, ^Li and ^^B were considerably underproduced in the standard 
cosmic ray spallation scenario. It is now believed that additional mechanisms contribute 
to their production. During the first billion years. Big Bang nucleosynthesis dominates 
the synthesis of ^Li (Sec. [3]), while later a stellar source, perhaps classical novae or 
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Figure 15. (Left panel) Number abundances of light elements in Galactic cosmic Rays 
(solid circles) and the Solar System (open squares), normalized to silicon. The cosmic 
ray values for H and He are adopted from measurements by the balloon-borne BESS 
instrument and the GSFC instrument on the IMP-8 spacecraft, while those for the 
other elements shown are from CRIS measurements. The cosmic ray abundances were 
measured during solar minimum at 170 MeV per nucleon |533j . The Solar System 
abundances are adopted from Lodders |539j . Note that in cosmic rays, B is more 
abundant than Li, which in turn is more abundant than Be. However, in the Solar 
System, Li is more abundant than B, which is more abundant than Be, supporting the 
conjecture that a significant fraction of the Solar System ''Li must have been produced 
by some stellar source (classical novae or AGB stars). (Right panel) Cross sections 
for the production of Li, Be and B by the spallation of protons incident on ^^C as a 
function of energy. Data are adopted from Ramaty et al. |532j . Note the decreasing 
sequence of B, Li and Be. 



AGB stars, contributed significantly to its production (Fig. [T5|) . Thus ^Li is exceptional 
among all naturally occurring nuclides in that its origin can be attributed to three 
distinct sites: the Big Bang, Galactic cosmic rays, and stars. In the case of ^^B it was 
found that neutrino- induced interactions with ^^C during core-collapse supernovae may 
also produce significant amounts of this species (Sec. E]). However, the predicted yields 
for the stellar production of ^Li and for the neutrino-induced synthesis of ^^B presently 
have large associated uncertainties. 



11.2. Galactic cosmic rays and spallation 

Galactic cosmic rays consist predominantly of energetic protons, a-particles and heavier 
nuclei, with kinetic energies up to the several 100 TeV per nucleon range and beyond. 
Their energy density amounts to about 2 eV cm~^ |535] and thus exceeds the energy 
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density of star light (~ 0.3 eV cm~^). Near an energy of 10 GeV per nucleon, for 
example, the number fractions of hydrogen, helium, CNO and LiBeB nuclei amount 
to 95%, 4.5%, 0.4% and 0.07%, respectively |536[ 1537] . An energetic cosmic ray 
particle propagating through the Galaxy faces a number of possibilities. First, it may 
undergo a high-energy nuclear collision (spallation) with interstellar nuclei. Second, an 
energetic cosmic ray particle may slow down significantly because of many collisions with 
interstellar electrons. Such ionization losses are responsible for the eventual deposition 
of cosmic ray matter into the interstellar medium. Third, the cosmic ray particle may 
escape into intergalactic space, depending on the value of the escape length, that is, 
the amount of matter the cosmic ray particle traverses between its source and the 
Galactic boundary. The escape length, which can be inferred from the abundance ratio 
of spallation products (for example, Li, Be and B) and target nuclei (mostly C, N and 
O), amounts to ~10 g cm~^ for Galactic cosmic rays near an energy of 1 GeV per 
nucleon ^538j . During propagation an energetic charged particle is strongly deflected by 
Galactic magnetic fields, such that most information regarding its original direction of 
motion is lost. 

Spallation may involve a proton or a-particle with energies in excess of several MeV 
per nucleon impinging on a heavy nucleus in the interstellar medium, of which C, N and 
O are most abundant. Alternatively, an energetic heavy nucleus, again mainly C, N or 
O, may interact with interstellar H and He. A third possibility, thought to be especially 
important for the production of the lithium isotopes in the early Galaxy when the CNO 
abundance was small, is the spallation of a-particles with energies in the tens of MeV per 
nucleon range with other ^He nuclei of the interstellar medium. The production rate of 
a given LiBeB species, i, expressed in terms of the number abundance ratio with respect 
to hydrogen, Fj = Ni/N^, is given by an expression of the form dYi/dt = T,j^k F^jajkY^Pi, 
with F the average cosmic ray flux (number of protons per cm^ and s), a the average 
spallation cross section (in cm^), and P the probability that after production species 
i will be thermalized and retained in the interstellar medium; the indices j and k 
denote the projectile and target, respectively, that is, protons, a-particles or CNO 
nuclei. It is instructive to consider a simple example: a typical cosmic ray proton flux 
of Fp ~10 cm~^ s~^, an average spallation cross section of cr^^QQ ~ 5 x 10~^^ cm^ [532] 
for the production of ^Be, assuming energetic protons impinging on C or O nuclei near 
an energy of 200 MeV per nucleon (i.e., the maximum of the spectrum), an abundance 
of YcNO ~ 10^^ in the interstellar medium, a probability of Pbs ~1, and integrating 
over a time period of ~ 10^° yr, yields a value of Y^e ~ 1-6 x 10~^^. The measured 
Solar System ^Be abundance from solar photosphere and meteoritic data amounts to 
'^Be — (2-1 ±0.1) X 10^^^ [539]. in good overall agreement with the predicted result. Of 
course, this simple estimate does not account for the change of the ^Be abundance over 
time. 
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11.3. Evolution of Li, Be, B abundances. Origin of Galactic cosmic rays 

The study of the Li, Be and B abundances, and in particular their evolution over 
the past several billion years, yields important clues for testing ideas regarding the 
nucleosynthesis in stars, for the physics of stellar outer layers, and for models of the 
early Galaxy. Moreover, the Li, Be and B abundance evolutions may hint at the origin 
of Galactic cosmic rays. Although it can be safely assumed that the sources of the 
cosmic ray particles, except perhaps those of the very highest energies, are located 
within the Galaxy, their origin and acceleration mechanism have not found a satisfactory 
explanation yet. In fact, the mystery of cosmic ray origin represents a long-standing 
problem of foremost importance. Energy arguments severely constrain the kind of 
candidate sources: for an assumed energy density of e = 2 eV cm~^, maintained on 
average over the volume of our Galaxy, V ^ 7 x 10^^ cm^, and a cosmic ray lifetime 
of T ~20 Myr near 1 GeV per nucleon, as derived from the cosmic ray abundances 
of the radioactive "clock" isotopes ^°Be, ^^Al, ^^Cl and ^^Mn [540] . one finds a total 
power of PcB. = ^V/r ~ 4 x 10"^" erg s~^. Thus, Galactic cosmic rays are expected 
to be associated with the most energetic astronomical objects. For example, a typical 
core-collapse supernova releases a kinetic energy of (1 — 2) x 10^^ erg (Sec. [6]). Assuming 
a frequency of about two such events per century gives a power of Psn ~ 10"^^ erg 
s~^, that is, significantly larger than what is needed for accelerating Galactic cosmic 
rays. Indeed, there is now substantial observational evidence that the origin of Galactic 
cosmic rays is related to supernova remnants [54111542] . 

Consider now the abundance evolution of ^Be and suppose that supernovae are the 
source of Galactic cosmic rays. One possibility for producing this species is energetic 
protons or a-particles colliding with interstellar C, N and O nuclei. The production rate 
of ^Be will then depend on the cosmic ray flux and the CNO abundance. The former 
factor is proportional to the supernova rate, dnsN/dt, while the latter is proportional 
to the total number of supernovae up to that time, usn- Consequently, the observed 
^Be abundance at any given instant should be proportional to n|^ or, in other words, 
to the square of the metallicity (as measured, for example, by the Fe abundance). 
On the other hand, ^Be may be produced by energetic CNO nuclei colliding with 
interstellar protons or a-particles. In this case, assuming a constant abundance of 
CNO nuclei in Galactic cosmic rays, the cosmic ray flux is again proportional to the 
supernova rate, but the H or He abundance in the interstellar medium is independent 
of nsN- Thus one would expect that the ^Be abundance at any given time is directly 
proportional to the metallicity |543[ 1544] . Alternatively, a linear relationship between 
the ^Be abundance and metallicity could naturally arise if supernova-produced protons 
or a-particles undergo spallation reactions with CNO nuclei in the same supernova's 
ejecta |545] . Similar arguments hold for the boron abundance evolution. 

Observations in the 1990s of beryllium and boron in metal-poor stars were of 
paramount importance in this regard |543l I546[ 1547] . Unlike lithium, which exhibits 
the Spite plateau at low metallicities (Sec. |3]), the beryllium and boron abundances 
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showed no such plateau and, furthermore, revealed a linear dependence on metallicity. 
Therefore, it is unlikely that Galactic cosmic rays originate from supernova-induced 
energetic protons or a-particles colliding with CNO nuclei of the interstellar medium 
(which would give rise to a beryllium or boron abundance proportional to the metallicity 
squared). Neither can the observations be explained by a CNO abundance in Galactic 
cosmic rays that, similar to the CNO abundance in the interstellar medium, would 
increase with time (since then the relationship between the beryllium or boron 
abundance and metallicity could not be linear). In addition, the CRIS instrument 
onboard the Advanced Composition Explorer (ACE) detected significant amounts of 
^^Co in Galactic cosmic rays, but only very small amounts of the radioactive precursor 
^^Ni |540j . The latter isotope decays by electron-capture (with a half-life of T1/2 = 
7.6 X 10^ yr in the laboratory), but is stable once fully stripped of its electrons and 
accelerated to high energies. From the observed abundance ratio of ^^Ni and ^^Co it can 
be estimated that the time delay between explosive nucleosynthesis and acceleration 
of radioactive ^^Ni must exceed ~ 10^ yr, when the ejecta of a particular supernova 
are presumably diluted and mixed in the interstellar medium. Consequently, it seems 
unlikely that individual supernovae accelerate their own ejecta. 

Several ideas have been proposed. A popular model is the production of Galactic 
cosmic rays in superbubbles, which result from the evolution of massive star clusters 
[5481 15491 I550J . First, through their strong winds, the massive stars create cavities of 
hot, low-density, gas in the interstellar medium. The cavities increase in size with each 
supernova explosion, and eventually merge to create superbubbles. Assuming that, in 
the simplest case, the composition of the metal-rich, low-density, material within the 
superbubble does not change with time. Galactic cosmic rays could be launched from this 
reservoir by subsequent supernova shocks. The superbubble origin has been criticized 
recently and other ideas have been proposed |551j . Future developments in this field at 
the intersection of nucleosynthesis, astroparticle physics and Galactic chemical evolution 
promise to be exciting. 

12. Conclusions and future perspectives 

The last fifty years have witnessed extraordinary progress in the modeling of 
stellar evolution, as well as in our understanding of the origin of the elements. 
First, supercomputers have allowed astrophysicists to perform complex hydrodynamic 
simulations coupled with extensive nuclear reaction networks at astounding levels 
of resolution. Indeed, last decades have seen the emergence of multidimensional 
simulations in many areas of stellar astrophysics research. Second, better observational 
constraints have been obtained through combined efforts in spectroscopy, the dawn of 
high-energy astrophysics with space-borne observatories, and laboratory measurements 
of presolar meteoritic grains that condensed in a suite of stellar environments. And 
third, nuclear physics has witnessed the birth of radioactive ion beam facilities, as well 
as direct measurements at astrophysically important energies using stable ion beam 
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facilities. 

Further progress eagerly awaits future developments aimed at solving key remaining 
questions in nuclear astrophysics, chief among them: 

(i) Why do predictions of helioseismology disagree with those of the standard solar 
model? 
What is the solution to the lithium problem in Big Bang nucleosynthesis? 

What do the observed light-nuclide and s-process abundances tell us about 
convection and dredge-up in massive stars and AGB stars? 

What are the production sites of the 7-ray emitting radioisotopes ^^Al, '^^Ti and 

60pg7 

What is the origin of about 30 rare and neutron deficient nuclides beyond the iron 
peak (p-nuclides)? 

What causes core-collapse supernovae to explode? 

What is the extend of neutrino-induced nucleosynthesis (//-process)? 

What is the extend of the nucleosynthesis in proton-rich outflows in the early ejecta 
of core-collapse supernovae (z/p-process)? 

What are the sites of the r-process? 

What causes the discrepancy between models and observations regarding the mass 
ejected during classical nova outbursts? 

Which are the physical mechanisms driving convective mixing in novae? 

What are the progenitors of type la supernovae? 

What is the nucleosynthesis endpoint in type I X-ray bursts? Is there any matter 

ejected from those systems? 
(XIV) What is the impact of stellar mergers on Galactic chemical abundances? 
[XV) What are the production and acceleration sites of Galactic cosmic rays? 

These questions will challenge young and courageous researchers for years to come, 
thereby keeping nuclear astrophysics in the spot light of science. 
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